Missing MQ tag from bwa mem alignments
1
0
Entering edit mode
22 months ago

I am using BWA MEM 0.7.17 for all my paired-end alignments.

Can someone clue me in on why some alignments have a MQ tag and others do not?

The commands are identical in some cases. Does BQSR apply this tag, there could be a difference in the options at this step perhaps? And I think some alignments use GATK 3.7 or earlier and others with the MQ tag has GATK 3.8

Thank you.

bwa mem wgs alignment • 694 views
ADD COMMENT
0
Entering edit mode
22 months ago
benformatics ★ 2.6k

Older versions of GATK (i.e. pre-4.x), and thus the bugs they contain, are no longer officially supported.

The missing MQ likely involves an incompatibility when combining VCFs. There doesn’t seem to be a retroactive solution.

https://gatkforums.broadinstitute.org/gatk/discussion/4618/mq-annotation-missing

ADD COMMENT
0
Entering edit mode

I was referring to SAM/BAM file MQ tag. Unfortunately sometimes I have to use legacy code.

ADD REPLY

Login before adding your answer.

Traffic: 2162 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6