Fasta.open(infile, mode="w") do |fasta_output|
Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
input.each_record do |record|
- if record.has_key? :SEQ
+ if record[:SEQ]
seq = Seq.new_bp(record)
total += record[:SEQ].length