The frontline support team will be slow on the forum because we are occupied with the GATK Workshop on March 21st and 22nd 2019. We will be back and more available to answer questions on the forum on March 25th 2019.
Missing variant while useing HaplotypeCaller
We are analysing monozygotic twins while using the HaplotypeCaller but we see some weird things. Not all SNPs are seen by the HaplotypeCaller while the bam file as input has clearly the SNP in there, this while in the twin it is called.
I did let the HaplotypeCaller output a bamfile and it seems that the haplotypecaller does not seem to use any reads on this SNP. In the igv screenshot there are the gvcf files of both samples and the pre-haplotypecaller bam and the post-haplotypecaller bam.
The pre-bams have used the indelrealignment and baserecalibration like in the best practice.
There multiple of this locations so this screenshot is just an example.
Someone know why this is caused and maybe how to solve this?