What does RAW_MQ annotation mean ?
REQUIRED for all errors and issues:
a) GATK version used: 4.3
b) Exact command used: I didn't create the file
c) Entire program log:
I have a VCF where the records have no RMS mapping quality (MQ) annotation but a RAW_MQ annotation (see example below). I've search the web and cannot find out what these values mean? I'm trying to determine what the mapping quality is like at these positions. Can I use these values at all ? Thank You - Robert
LKNJ01018100.1 585 . A G 9674.13 PASS AC=22;AF=0.611;AN=36;BaseQRankSum=1.03;ClippingRankSum=0;DP=500;ExcessHet=3.01;FS=53.516;InbreedingCoeff=-0.6723;MLEAC=21;MLEAF=0.583;MQRankSum=0;QD=27.17;RAW_MQ=95735;ReadPosRankSum=-0.042;SOR=2.007 GT:AD:DP:GQ:PL 0/1:2,15:17:26:496,0,26 0/1:4,15:19:71:515,0,71 0/1:10,12:22:99:346,0,264 0/1:4,11:15:96:396,0,96 0/1:8,9:17:99:312,0,276 1/1:1,17:18:13:616,13,0 0/1:4,16:20:39:519,0,39 1/1:1,13:14:1:485,1,0 0/1:4,13:17:94:429,0,94 1/1:0,27:27:81:1065,81,0 0/1:5,14:19:99:510,0,123 1/1:1,32:33:63:1215,63,0 0/1:6,10:16:99:296,0,106 0/1:7,12:19:99:396,0,249 0/1:2,18:20:18:602,0,18 0/1:4,18:22:58:572,0,58 0/1:7,13:20:99:430,0,251 0/1:6,15:21:99:548,0,162
-
Thank you for your post, Robert Baldwin! 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.
-
Hi Robert Baldwin,
Thank you for writing to the GATK forum! I hope that we can help you sort this out.
I brought your inquiry to our developers and received some feedback to share with you.
Firstly, would you be able to clarify how exactly you created your VCF? It would be helpful if you try checking inside the header of the VCF for all of the commands you have run on the file so far. Once you locate this, please include it in your response to this post, and we can move forward with diagnosing the problem.
It might look something like this:
##GATKCommandLine.SelectVariants=<ID=SelectVariants,Version=3.4-46-gbc02625,Date="Fri Sep 25 10:45:27 EDT 2015",Epoch=1443192327688,CommandLineOptions="analysis
I hope to hear back from you soon! Please include any further questions that may arise in your response.
Best,
Anthony
-
Hi Robert Baldwin,
We haven't heard from you in a while so we're going to close out this ticket. If you still require assistance, simply respond to this email and we'll be happy to pick up where we left off!
Kind regards,
Anthony
Please sign in to leave a comment.
3 comments