bio-polyploid-tools 0.4.4 → 0.4.5

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: c65f3649faacf8bbf2fdb2269a6ed125841ba74c
-  data.tar.gz: a3ef8e79ccc3feeaece52dca1dafdd1bf74888cd
+  metadata.gz: 9e11509fcfcbe8aa33caccd377dd0eca756713ef
+  data.tar.gz: ee00b6ac0ca58b233e3154caf2cc3e00b909a869
 SHA512:
-  metadata.gz: 0fe43f7748842d15ac6e08feab09c7fba77b7fe9b7879a091df5c622cca453f5d7878c6ae39447dc4bb26ca4239d4a80700d69a84fb20348d7891d7be68d689a
-  data.tar.gz: 257ec9bc355fa195148b67600d0507edaff4b5144c1a05db4b95ac468bede2ca09dab8b459dd2131d0f98a23018bf777d6fe8e819084bc173049c8d4e0ab8878
+  metadata.gz: 50afb04d54bd58599d730cfe44bcf7f8e51ec2aea62af870ebca9af0e685fc7a68e06da8a82b82973f1f8dcd6cb3f3414b1b93753dd64e37b19d12473b46adc1
+  data.tar.gz: 3b2313544ce3954ca7e6f0c6ca4e4de32255329468511a662ea15abb87e5cbcc038328371d93d3213c6daba56b453bb2f0bac5660e3009010414f4284fdac2a9

data/VERSION

@@ -1 +1 @@
-0.4.4
+0.4.5

data/bio-polyploid-tools.gemspec

@@ -2,11 +2,11 @@
 # DO NOT EDIT THIS FILE DIRECTLY
 # Instead, edit Jeweler::Tasks in Rakefile, and run 'rake gemspec'
 # -*- encoding: utf-8 -*-
-# stub: bio-polyploid-tools 0.4.4 ruby lib
+# stub: bio-polyploid-tools 0.4.5 ruby lib
 
 Gem::Specification.new do |s|
   s.name = "bio-polyploid-tools"
-  s.version = "0.4.4"
+  s.version = "0.4.5"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.require_paths = ["lib"]

data/lib/bio/PolyploidTools/SNP.rb

@@ -265,7 +265,7 @@ module Bio::PolyploidTools
       str << "=\n"
 
 
-      #In case that we don't have a right primer, wi do both orientation
+      #In case that we don't have a right primer, we do both orientations
       unless opts[:right_pos]
         sequence =  opts[:sequence]    
         left = sequence.size - left - 1
@@ -290,6 +290,7 @@ module Bio::PolyploidTools
       seq_original = String.new(pr.sequence)
       #puts seq_original.size.to_s << "-" << primer_3_min_seq_length.to_s
       return primer_3_propertes if seq_original.size < primer_3_min_seq_length
+      
       seq_original[pr.snp_pos] = self.original
       seq_original_reverse = reverse_complement_string(seq_original)
 

data/lib/bio/db/primer3.rb

@@ -141,11 +141,11 @@ module Bio::DB::Primer3
           values << primer3_line_1.right_primer 
           values << primer3_line_1.type.to_s 
           values << primer3_line_1.orientation.to_s 
-          values << primer3_line_1.shortest_pair.left.tm 
+          values << primer3_line_1.best_pair.left.tm 
           values << primer_2_tm
-          values << primer3_line_1.shortest_pair.right.tm
+          values << primer3_line_1.best_pair.right.tm
           values << "first" 
-          values << primer3_line_1.shortest_pair.product_size
+          values << primer3_line_1.best_pair.product_size
         elsif  primer_1_tm != "NA"
           values << primer_1
           values << primer3_line_2.left_primer
@@ -153,60 +153,62 @@ module Bio::DB::Primer3
           values << primer3_line_2.type.to_s
           values << primer3_line_2.orientation.to_s
           values << primer_1_tm
-          values << primer3_line_2.shortest_pair.left.tm
-          values << primer3_line_2.shortest_pair.right.tm
+          values << primer3_line_2.best_pair.left.tm
+          values << primer3_line_2.best_pair.right.tm
           values << "second"
-          values << primer3_line_2.shortest_pair.product_size
+          values << primer3_line_2.best_pair.product_size
         else
 
           first_candidate = find_primer_pair_first
           second_candidate = find_primer_pair_second
 
           if first_candidate
-
-            #puts "input to search #{first_candidate.left_coordinates}"
             primer_2 = primer3_line_2.left_primer_with_coordinates(first_candidate.left_coordinates, first_candidate.orientation)
             primer_2_tm = find_left_primer_temp(primer_2)
-            #puts "In the funky if #{primer_2}"
           end
           if second_candidate
+            #puts "input to search #{first_candidate.left_coordinates}"
             primer_1 = primer3_line_1.left_primer_with_coordinates(second_candidate.left_coordinates, second_candidate.orientation) 
             primer_1_tm = find_left_primer_temp(primer_1)
+            #puts "In the other funky if #{primer_2}"
           end
 
           if first_candidate and second_candidate and first_candidate < second_candidate 
+            #puts "A"
             values << first_candidate.left_primer
             values << primer_2
             values << first_candidate.right_primer 
             values << first_candidate.type.to_s 
             values << first_candidate.orientation.to_s 
-            values << first_candidate.shortest_pair.left.tm 
+            values << first_candidate.best_pair.left.tm 
             values << primer_2_tm
-            values << first_candidate.shortest_pair.right.tm
+            values << first_candidate.best_pair.right.tm
             values << "first" 
-            values << first_candidate.shortest_pair.product_size
+            values << first_candidate.best_pair.product_size
           elsif  second_candidate 
+            #puts "B"
             values << primer_1
             values << second_candidate.left_primer
             values << second_candidate.right_primer
             values << second_candidate.type.to_s
             values << second_candidate.orientation.to_s
             values << primer_1_tm
-            values << second_candidate.shortest_pair.left.tm
-            values << second_candidate.shortest_pair.right.tm
+            values << second_candidate.best_pair.left.tm
+            values << second_candidate.best_pair.right.tm
             values << "second"
-            values << second_candidate.shortest_pair.product_size
+            values << second_candidate.best_pair.product_size
           elsif  first_candidate 
-            values << primer_2
+            #puts "C"
             values << first_candidate.left_primer
+            values << primer_2
             values << first_candidate.right_primer
             values << first_candidate.type.to_s
             values << first_candidate.orientation.to_s
             values << primer_2_tm
-            values << first_candidate.shortest_pair.left.tm
-            values << first_candidate.shortest_pair.right.tm
+            values << first_candidate.best_pair.left.tm
+            values << first_candidate.best_pair.right.tm
             values << "first"
-            values << first_candidate.shortest_pair.product_size
+            values << first_candidate.best_pair.product_size
 #          else
 #            values << "" 
           end
@@ -220,12 +222,12 @@ module Bio::DB::Primer3
         values << primer3_line_1.right_primer 
         values << primer3_line_1.type.to_s 
         values << primer3_line_1.orientation.to_s      
-        values << primer3_line_1.shortest_pair.left.tm 
+        values << primer3_line_1.best_pair.left.tm 
         values << "NA"
-        values << primer3_line_1.shortest_pair.right.tm
+        values << primer3_line_1.best_pair.right.tm
 
         values << "first+"
-        values << primer3_line_1.shortest_pair.product_size
+        values << primer3_line_1.best_pair.product_size
       elsif primer3_line_2 
         values << primer3_line_2.polymorphism
         values << primer3_line_2.left_primer_snp(self) 
@@ -234,10 +236,10 @@ module Bio::DB::Primer3
         values << primer3_line_2.type.to_s
         values << primer3_line_2.orientation.to_s
         values << "NA"
-        values << primer3_line_2.shortest_pair.left.tm
-        values << primer3_line_2.shortest_pair.right.tm
+        values << primer3_line_2.best_pair.left.tm
+        values << primer3_line_2.best_pair.right.tm
         values << "second+"
-        values << primer3_line_2.shortest_pair.product_size
+        values << primer3_line_2.best_pair.product_size
 
       end 
       values.join(",")
@@ -274,7 +276,9 @@ module Bio::DB::Primer3
         end
       case
       when primer3record.line == @line_1
+
         @line_1_template = primer3record.sequence_template
+
       when primer3record.line == @line_2
         @line_2_template = primer3record.sequence_template
       else
@@ -283,10 +287,11 @@ module Bio::DB::Primer3
 
       if  primer3record.primer_left_num_returned.to_i > 0 
         case
-        when primer3record.line == @line_1
+        when primer3record.line == @line_2
           primers_line_1 << primer3record
+          #puts primer3record.inspect
           @primer3_line_1 = primer3record if not @primer3_line_1  or @primer3_line_1 > primer3record
-        when primer3record.line == @line_2
+        when primer3record.line == @line_1
           primers_line_2 << primer3record
           @primer3_line_2 = primer3record if not @primer3_line_2 or @primer3_line_2 > primer3record
         else
@@ -299,15 +304,18 @@ module Bio::DB::Primer3
   class Primer3Record
     include Comparable
     attr_accessor :properties, :polymorphism
+    attr_accessor :socres
 
-    def shortest_pair
-      return @shortest_pair if @shortest_pair
-      @shortest_pair = nil
-      @primerPairs.each do | primer |
-        @shortest_pair = primer if @shortest_pair == nil
-        @shortest_pair = primer if primer.size < @shortest_pair.size
-      end
-      @shortest_pair
+
+    def best_pair
+      return @best_pair if @best_pair
+      #@best_pair = nil
+      #@primerPairs.each do | primer |
+      #  @best_pair = primer if @best_pair == nil
+      #  @best_pair = primer if primer.size < @best_pair.size
+      #end
+      @best_pair = @primerPairs.first
+      @best_pair
     end
 
     def primer_error
@@ -335,44 +343,16 @@ module Bio::DB::Primer3
       return nil
     end
 
-    def <=>(anOther)
-      ret = snp <=> anOther.snp
-      return ret if ret != 0
-
-
-      #Sorting by the types. 
-      if type == :chromosome_specific 
-        if anOther.type != :chromosome_specific
-          return -1
-        end
-      elsif type == :chromosome_semispecific
-        if anOther.type == :chromosome_specific
-          return 1
-        else anOther.type == :chromosome_nonspecific
-          return -1
-        end
-      elsif type == :chromosome_nonspecific
-        if anOther.type != :chromosome_nonspecific
-          return 1
-        end
-      end
-
-      #Sorting if it is in intron or not This will give priority 
-      #to the cases when we know for sure the sequence from the line
-      #and reduce the chances of getting messed with a short indel
-      if self.exon?
-        unless anOther.exon? 
-          return -1
-        end
-      else
-        if anOther.exon?
-          return 1
-        end
-      end
-
-      #Sorting for how long the product is, the shorter, the better 
-      return  product_length <=> anOther.product_length
+    def score
+      ret = 0
+      ret += @scores[type]
+      ret += @scores[:exon] if exon?
+      ret -= product_length
+      ret
+    end
 
+    def <=>(anOther)
+      return  anOther.score <=> score
     end
 
     def parse_coordinates(str)
@@ -385,13 +365,13 @@ module Bio::DB::Primer3
 
     def left_coordinates
       #@left_coordinates = parse_coordinates(self.primer_left_0) unless @left_coordinates 
-      @left_coordinates = shortest_pair.left.coordinates
+      @left_coordinates = best_pair.left.coordinates
       @left_coordinates 
     end
 
     def right_coordinates
       unless @right_coordinates 
-        @right_coordinates = shortest_pair.right.coordinates
+        @right_coordinates = best_pair.right.coordinates
         @right_coordinates[0] = @right_coordinates[0] - @right_coordinates[1] + 1
       end
       @right_coordinates 
@@ -399,7 +379,7 @@ module Bio::DB::Primer3
 
     def left_primer
       #@left_primer = self.sequence_template[left_coordinates[0],left_coordinates[1]] unless @left_primer
-      @left_primer = shortest_pair.left.sequence
+      @left_primer = best_pair.left.sequence
       @left_primer
     end
 
@@ -448,15 +428,21 @@ module Bio::DB::Primer3
     end
     
     def right_primer
-      return shortest_pair.right.sequence
+      return best_pair.right.sequence
     end
 
     def product_length
-      return shortest_pair.size
+      return best_pair.size
     end
 
     def initialize
       @properties = Hash.new
+      @scores = Hash.new
+      @scores[:chromosome_specific] = 1000
+      @scores[:chromosome_semispecific] = 100
+      @scores[:chromosome_nonspecific] = 0
+      @scores[:exon] = 50
+    
     end
 
     def snp
@@ -587,7 +573,7 @@ module Bio::DB::Primer3
   end
 
   class PrimerPair
-
+    include Comparable
     attr_reader :record
     attr_reader :left, :right
     
@@ -602,6 +588,10 @@ module Bio::DB::Primer3
       return product_size.to_i
     end
 
+    def <=>(anOther)
+      penalty.to_f <=> anOther.penalty.to_f
+    end
+
     def initialize(record, index)
       raise Primer3Exception.new(), "Index #{index} is greater than the number of records" unless index < record.size
       @record = record

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: bio-polyploid-tools
 version: !ruby/object:Gem::Version
-  version: 0.4.4
+  version: 0.4.5
 platform: ruby
 authors:
 - Ricardo H.  Ramirez-Gonzalez