gatk
gatk copied to clipboard
Merging phased indels
We are performing somatic variant calling in tumour-only mode where we call variants with both Mutect2 (GATK version 4.1.2.0) and Vardict-java. We have noticed that Mutect2 calls adjacent phased indels as two separate variants, whereas Vardict calls these as a single variant. We have tried increasing the value of the --max-mnp-distance option, and this merges SNP-SNP combinations into an MNP but does not merge phased combinations of INDEL-SNP or INDEL-INDEL. Is there any way this could be implemented?