fixing broken biopieces

[biopieces.git] / bp_bin / patscan_seq

diff --git a/bp_bin/patscan_seq b/bp_bin/patscan_seq
index 194aaae7c812c42f6790d57272cf858bd8918d4f..1b533a9ebee71f22ed6b769371bcc7b607d4bcfa 100755 (executable)
--- a/bp_bin/patscan_seq
+++ b/bp_bin/patscan_seq
@@ -146,14 +146,16 @@ Biopieces.open(options[:stream_in], tmp_file) do |input, output|
         record[:SEQ_NAME] = seq_name_count
         seq_name_count += 1
 
+        seq = Seq.new_bp(record)
+
         if seq_type.nil?
-          seq      = Seq.new("", record[:SEQ])
           seq_type = seq.type_guess
         end
+
+        fasta_io.puts seq.to_fasta
       end
 
       output.puts record
-      fasta_io.puts record
     end
   end
 end