Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,498,176 | A→G | 8.0% | I117V (ATT→GTT) | cysG → | fused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,498,176 | 0 | A | G | 8.0% | 160.3 / 4.1 | 50 | I117V (ATT→GTT) | cysG | fused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (25/21); new base G (2/2); total (27/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 |
CAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGT > NC_000913/3498086‑3498321 | cAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt < 4:441081/110‑1 (MQ=255) cAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt > 3:441081/1‑110 (MQ=255) cAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGtc > 1:423267/1‑120 (MQ=255) cAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGtc < 2:423267/120‑1 (MQ=255) cgcgTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt > 4:667809/1‑107 (MQ=255) cgcgTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt < 3:667809/107‑1 (MQ=255) cgCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACt > 1:101679/1‑140 (MQ=255) cgCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACt < 2:101679/140‑1 (MQ=255) tctGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCa > 1:570505/1‑140 (MQ=255) tctGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCa < 2:570505/140‑1 (MQ=255) tGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGcc < 1:663032/115‑1 (MQ=255) tGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGcc > 2:663032/1‑115 (MQ=255) gTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTTCTGCCGtt > 2:920426/1‑148 (MQ=255) tAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCActgct < 4:123138/141‑1 (MQ=255) tAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCActgct > 3:123138/1‑141 (MQ=255) tCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt > 4:1130879/1‑61 (MQ=255) tCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCAt < 3:1130879/61‑1 (MQ=255) cgcCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGt < 4:140088/132‑1 (MQ=255) cgcCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGt > 3:140088/1‑132 (MQ=255) aaaGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTg > 3:816429/1‑67 (MQ=255) aaaGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTg < 4:816429/67‑1 (MQ=255) cAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAAt > 1:1004141/1‑107 (MQ=255) cAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAAt < 2:1004141/107‑1 (MQ=255) tttATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCg > 1:328395/1‑135 (MQ=255) tttATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCg < 2:328395/135‑1 (MQ=255) ttattaTGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTg > 4:499241/1‑74 (MQ=255) ttattaTGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTg < 3:499241/74‑1 (MQ=255) ttatGCCGTCGATTGTTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAAtt > 3:99990/1‑148 (MQ=255) ttatGCCGTCGATTGTTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAAtt > 3:98838/1‑148 (MQ=255) tGCCGTCGATTGTTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTAcg < 4:98838/148‑1 (MQ=255) tGCCGTCGATTGTTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTAcg < 4:99990/148‑1 (MQ=255) tccGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCAct > 2:901752/2‑99 (MQ=255) gCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCAct < 1:901752/99‑1 (MQ=255) cGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTg > 3:319875/1‑113 (MQ=255) cGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTg < 4:319875/113‑1 (MQ=255) cGTCGATTATTGACCGCTCACCGCTCATAGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCt > 3:932003/1‑112 (MQ=255) cGTCGATTATAGACCGCTCACCGCTCATAGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGACACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCt < 4:932003/112‑1 (MQ=255) gTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGc < 3:42453/66‑1 (MQ=255) gTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGc > 4:42453/1‑66 (MQ=255) gTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCggg > 4:126220/1‑148 (MQ=255) tCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGa > 1:690600/1‑89 (MQ=255) tCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGa < 2:690600/89‑1 (MQ=255) tCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGc > 3:955687/1‑148 (MQ=255) tCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGc > 2:412086/1‑148 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCtctc < 4:686786/54‑1 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCtctc > 3:686786/1‑54 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTgcg < 4:353090/76‑1 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTgcg > 3:353090/1‑76 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACg > 4:624815/1‑129 (MQ=255) cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACg < 3:624815/129‑1 (MQ=255) ttattGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGt < 1:920426/148‑1 (MQ=255) | CAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGT > NC_000913/3498086‑3498321 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |