bp[:SEQ] = sam[:SEQ].seq
unless sam[:SEQ].qual.nil?
- bp[:SCORES] = sam[:SEQ].convert_phred2illumina!.qual
+ bp[:SCORES] = sam[:SEQ].qual_convert!(:base_33, :base_64).qual
end
if sam[:NM] and sam[:NM].to_i > 0
entry[:RNEXT] = rnext
entry[:PNEXT] = pnext
entry[:TLEN] = tlen
- entry[:SEQ] = (qual == '*') ? Seq.new(qname, seq) : Seq.new(qname, seq, qual)
+ entry[:SEQ] = (qual == '*') ? Seq.new(seq_name: qname, seq: seq) : Seq.new(seq_name: qname, seq: seq, qual: qual)
entry[:QUAL] = qual
# Optional fields - where some are really important! HATE HATE HATE SAM!!!