This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
Combine phased calls from Mutect2
Hello, is there a way to have Mutect2 emit multi-nucleotide variants instead of multiple adjacent SNVs?
For example, consider this variant:
Mutect will call the G/T SNP in position 2 as one line, and the G/C SNP at position 3 as another line. Then the fact that they are part of the same haplotype is indicated by the phasing information in the info column of the vcf.
I would prefer to have it call a multiple nucleotide variant: REF GG and ALT TC.
Can I get Mutect to do this? Or is there any post-processing tool you can recommend?
The reason for preferring MNVs instead of SNVs is that I am using ensembl-VEP to predict the protein consequences of the variants. In that case it's quite important to represent the actual haplotypes instead of stepping through variant sites one-by-one.