diff --git a/src/alignment_io.cpp b/src/alignment_io.cpp index d4dd3da..4436d36 100644 --- a/src/alignment_io.cpp +++ b/src/alignment_io.cpp @@ -913,6 +913,12 @@ gafkluge::GafRecord alignment_to_gaf(function node_to_length, string cs_cigar_str = "+" + aln.sequence(); gaf.opt_fields["cs"] = make_pair("Z", std::move(cs_cigar_str)); } + + // emit base qualities for unaligned sequences + // optional base qualities + if (base_quals && !aln.quality().empty()) { + gaf.opt_fields["bq"] = make_pair("Z", string_quality_short_to_char(aln.quality())); + } } // optional frag_next/prev names