Mutect2 -could not detect expected Indels
REQUIRED for all errors and issues:
a) GATK version used:
b) Exact command used:
c) Entire program log:
Hi GATK team,
I am calling somatic variants using Mutect2 for 10 samples. I could see lots of missense mutations but could not find expected small Indels there in all 10 samples. What may be the reason ? I think Mutect2 is able to detect both SNVs and Indels well. Please clarify and let me know if I need to tweak some parameters.
a) GATK version used: 4.2.2.0
b) Exact command: python3 gatk mutect2 -R hg19.fasta -I sample1.bam --max-mnp-distance 0 -o sample1.vcf.gz.
-
Mutect2 and HaplotypeCaller of GATK utilizes local reassembly for genotyping of variants and the default parameters of both tools are mostly optimized for whole genome and exome sequencing samples and therefore may not be enough for other special NGS data types such as amplicon sequencing. Depth of data may also be an issue for some variants however these tools have bunch of parameters for tuning for those cases and some of the most useful parameters are summarized in the following article
Can you try some of the recommended tunable parameters and let us know whether that works for you?
Regards.
-
Okay, thanks. I am looking into it and will get back to you soon.
Thanks,
Shivangi
Please sign in to leave a comment.
2 comments