fixed broken pcr_seq

[biopieces.git] / bp_bin / pcr_seq

diff --git a/bp_bin/pcr_seq b/bp_bin/pcr_seq
index dec3ea0e4f1b4c01c6841a9499a5c47394248ec5..f54577373e75f7f32f24752c3a54ff31eaedafa8 100755 (executable)
--- a/bp_bin/pcr_seq
+++ b/bp_bin/pcr_seq
@@ -164,7 +164,11 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
   Fasta.open(infile, mode="w") do |ios|
     input.each_record do |record|
       output.puts record
-      ios.puts record
+
+      if record.has_key? :SEQ
+        entry = Seq.new_bp(record)
+        ios.puts entry.to_fasta
+      end
     end
   end