Fasta.open(infile, mode="w") do |fasta_io|
input.each_record do |record|
output.puts 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
prodigal.run
prodigal.each do |record|
- record[:DEFINITION].match /seqhdr="([^"]+)/ do |m|
- record[:S_ID] = $1
- end
-
output.puts record
end
end