We focused on identifying heteroplasmy levels in nanopore sequencing data using Mutserve2 with default settings. We simulated various levels of heteroplasmy for one mitochondrial variant using bioinformatics. What we observed is that the proportion of called variants starts to decline at a heteroplasmy ratio between 0.93 and 0.97, with values mainly ranging between 0.93 and 0.94.
This suggests that the variant, which is present in more than 90% of reads, wasn't called in about 70 out of 999 attempts. Furthermore, this variant is in the .txt file but flagged as strand bias, so it's filtered out and doesn't appear in the final VCF file.
Do you have any ideas to explain this anomaly? Are there specific parameters that can be adjusted to address this issue?
We focused on identifying heteroplasmy levels in nanopore sequencing data using Mutserve2 with default settings. We simulated various levels of heteroplasmy for one mitochondrial variant using bioinformatics. What we observed is that the proportion of called variants starts to decline at a heteroplasmy ratio between 0.93 and 0.97, with values mainly ranging between 0.93 and 0.94.
This suggests that the variant, which is present in more than 90% of reads, wasn't called in about 70 out of 999 attempts. Furthermore, this variant is in the .txt file but flagged as strand bias, so it's filtered out and doesn't appear in the final VCF file.
Do you have any ideas to explain this anomaly? Are there specific parameters that can be adjusted to address this issue?