fgbio icon indicating copy to clipboard operation
fgbio copied to clipboard

SetMateInformation should add the MQ tag to any supplementary alignments

Open msto opened this issue 1 year ago • 1 comments

Currently, SetMateInformation does not add the MQ tag to supplementary alignments.

I'm not sure if this is intended behavior - the docs suggest that supplementary alignments would have this tag added.

Adds and/or fixes mate information on paired-end reads. Sets the MQ (mate mapping quality), 'MC' (mate cigar string), ensures all mate-related flag fields are set correctly, and that the mate reference and mate start position are correct.

Supplementary records are handled correctly (updated with their mate's non-supplemental attributes). Secondary alignments are passed through but are not updated.

GroupReadsByUmi requires that all alignments, including supplementary alignments, have the MQ tag set, so it would be helpful if SetMateInformation produced a compatible BAM 🙂

msto avatar Jan 26 '24 17:01 msto