Missing mutations
I am using the latest version of Mutect2 (4.1.9.0) and trying to call mutations from Mutect. I noticed that Mutect is not calling muations with clear read support in IGV. I tried some suggestions found on this forum like this (https://gatk.broadinstitute.org/hc/en-us/community/posts/360059696811-Mutect2-not-calling-a-4-bp-deletion-in-BRCA1-with-50-AF) but they didn't help. You can see the mutation here:
This is a call for a BAM file that doesn't have an associated normal file. I run Mutect2 and this is the command through Cromwell:
```
${gatk_path} \
Mutect2 \
-R ${ref_fasta} \
-I ${input_bam} \
--panel-of-normals ${panel} \
--germline-resource ${gnomad} \
--f1r2-tar-gz f1r2.tar.gz \
-L ${regions_list} \
-O ${output_filename}
```
The regions list that I use is simply a file that splits the reference into different regions to call. I do not specify particular regions to call, simply to break up the calling into chunks.
There should be a mutation somewhere between chr13:32,914,430-32,914,440 in the photo, I think it is at 32,914,440 based on the photo. Is there anything else I can try?
-
I read that using `force-active` in Mutect2 may help. What exactly is this doing, and why is this turned off by default?
-
Hi vctrymao, please see this troubleshooting document and work through these steps: When HaplotypeCaller and Mutect2 do not call an expected variant
In terms of the --force-active argument, it is turning on all regions as active regions. Please see this article about active regions for more information.
Also, is this region in IGV near any of the ends of your regions-list intervals?
-
I read this, but I did not see an exact description of what an "active region" is defined as mathematically. Regardless, this seems to have fixed the problem.
-
Thanks for posting the update vctrymao!
-
Could you explain how exactly GATK defines active regions? It isn't intuitive to me that Mutect2 would omit a region with high read level support.
-
vctrymao please see this paper: https://github.com/broadinstitute/gatk/blob/master/docs/mutect/mutect.pdf There is a section on the active regions and how they are determined mathematically.
-
vctrymao -force-active makes Mutect2 run very slowly. It is possible this variant is not output because it overlaps the germline resource or the panel of normals.
Please sign in to leave a comment.
7 comments