Fig. 1

Heteroplasmy detection workflow. The raw read files obtained from the ChIP-seq experiments were first aligned to the respective reference genomes. The aligned reads were then pre-processed, filtering out duplicate reads and extracting reads mapping with a high quality score to the mtDNA. The heteroplasmy detection algorithm was then used across the samples. Finally, we analyzed the genomic properties of heteroplasmic positions across vertebrates