fixed bug in plot_scores

[biopieces.git] / bp_bin / pcr_seq

diff --git a/bp_bin/pcr_seq b/bp_bin/pcr_seq
index 5cee30b152e49d2d0a6c7013c6d05952b80bc57a..f6c874e7e9555e021c83f5213f6277e2f8f591ff 100755 (executable)
--- a/bp_bin/pcr_seq
+++ b/bp_bin/pcr_seq
@@ -146,7 +146,7 @@ class Pattern
     seq      = Seq.new
     seq.seq  = primer
     seq.type = 'dna'
-    seq.revcomp
+    seq.reverse!.complement!
 
     descriptor ? seq.seq + descriptor : seq.seq
   end
@@ -164,11 +164,11 @@ tmpdir  = Biopieces.mktmpdir
 infile  = File.join(tmpdir, "in.fna")
 
 if options[:forward_rc]
-  options[:forward] = Seq.new("test", options[:forward_rc], 'dna').revcomp.seq
+  options[:forward] = Seq.new("test", options[:forward_rc], 'dna').reverse.complement.seq
 end
 
 if options[:reverse_rc]
-  options[:reverse] = Seq.new("test", options[:reverse_rc], 'dna').revcomp.seq
+  options[:reverse] = Seq.new("test", options[:reverse_rc], 'dna').reverse.complement.seq
 end
 
 raise ArgumentError, "no adaptor specified" unless options[:forward] or options[:reverse]
@@ -177,7 +177,7 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
     input.each_record do |record|
       output.puts record
 
-      if record.has_key? :SEQ
+      if record[:SEQ]
         entry = Seq.new_bp(record)
         ios.puts entry.to_fasta
       end