Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
Fasta.open(file_fasta, mode="w") do |fasta_io|
input.each_record do |record|
- fasta_io.puts record
+ if record[:SEQ_NAME] and record[:SEQ]
+ seq = Seq.new_bp(record)
+
+ fasta_io.puts seq.to_fasta
+ end
end
end