+ Fasta.open(@outfile, 'r') do |ios|
+ ios.each do |entry|
+ record = {}
+
+ fields = entry.seq_name.split(" # ")
+
+ record[:REC_TYPE] = "GENE"
+ record[:S_ID] = fields[0]
+ record[:S_BEG] = fields[1].to_i - 1
+ record[:S_END] = fields[2].to_i - 1
+ record[:S_LEN] = record[:S_END] - record[:S_BEG] + 1
+ record[:STRAND] = fields[3] == '1' ? '+' : '-'
+ record[:SEQ] = entry.seq
+ record[:SEQ_LEN] = entry.length
+
+ yield record