Mutect2 & FilterMutectCalls: Kind request for exact derivation of MBQ values in INFO column of VCFs

Follow

Franziska

• December 06, 2022 13:45

Dear gatk community,

for QC I was extracting the MBQ values of my somatic variantc calls from Mutect2 (gatk/4.2.0.0, tumor-normal mode + FilterMutectCalls, from INFO column of VCF). I assumed MBQ column is giving median base quality of ref and median base quality of alt allele (seperated by comma), meaning it gives me the median of all phred scaled sequencing base qualities at the variant position of all reads considered for the ref or alt allele. Is this right?

I am asking, because some of my data shows significant differences between MBQ distributions of ref and alternative alleles, which I cannot reproduce if I extract variant and reference reads with gatk PrintReads and extract the median base qualities manually.

Thank you very much for clarifying the exact derivation of MBQ values given in the INFO column of Mutect2 VCFs.

Best regards,

Franziska

 

 

 

 

0

• 

  Genevieve Brandt (she/her)

  • December 07, 2022 22:08

  Thank you for your post, Franziska! I want to let you know we have received your question. We'll get back to you if we have any updates or follow up questions. 

  Please see our Support Policy for more details about how we prioritize responding to questions. 

  0

  Comment actions Permalink
• 

  Franziska

  • January 25, 2023 10:40

  Are there any news regarding my question?

  0

  Comment actions Permalink

Please sign in to leave a comment.