Hi All, I am analysing mitochondrial genomic sequence data sets to perform mutation rate analysis among my samples. I performed variant calling by using bcftools and have vcf files. At the time of filtration I need to consider the heteroplasmy state of the mitochondria among the samples. Will keeping the variants supported by low allele depth address the heteroplasmy state?
I am really confused as mitochondria are haploid and heteroplasmic what will be the combination of parameters for variant filteration?
Also can anyone explain me the basic of how genotpye and allele depth related in vcf file?