"read_sff -i #{@sff_file} |
progress_meter |
analyze_seq |
- analyze_vals -k GC%,HARD_MASK%,SOFT_MASK% |
+ analyze_vals -k SEQ,GC%,HARD_MASK%,SOFT_MASK% |
write_tab -o #{@anal_file} -x"
)
STDERR.puts "done.\n"
fields = line.split("\t")
case fields.first
- when "SEQ" then @count, @min, @max, @sum, @mean = fields[2], fields[3], fields[4], fields[5], fields[6]
+ when "SEQ" then @count, @min, @max, @sum, @mean = fields[2 .. 6]
when "GC%" then @gc = fields[6]
when "HARD_MASK%" then @hard = fields[6]
when "SOFT_MASK%" then @soft = fields[6]