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?
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 mergesSNP-SNP
combinations into an MNP but does not merge phased combinations ofINDEL-SNP
orINDEL-INDEL
. Is there any way this could be implemented?