Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,915,756 | A→G | 12.5% | L111P (CTG→CCG) | ccmC ← | protoheme IX reservoir complex subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,915,756 | 0 | A | G | 12.5% | 77.8 / 3.6 | 32 | L111P (CTG→CCG) | ccmC | protoheme IX reservoir complex subunit |
Reads supporting (aligned to +/- strand): ref base A (17/11); new base G (2/2); total (19/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01 |
GCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGA > NC_002947/4915612‑4915903 | gCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGTATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAgcg > 1:80882/1‑149 (MQ=255) cGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAAAGCCGAACTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGAAAACCAGCGCCa > 1:92169/1‑149 (MQ=255) tCACTATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAATGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCAcggcgg > 5:212986/1‑148 (MQ=255) aCGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCa > 5:253930/1‑149 (MQ=255) cTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgc < 6:75433/149‑1 (MQ=255) aTCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcgc > 6:41032/1‑148 (MQ=255) gCCTGGCCCAGCGGAATTATGCCGTAGTTCAGGAACAGCAGG‑ATGAGCATGGAGGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCTGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCAGGGCGGTCATCCAGGCGCCaa < 2:92169/148‑1 (MQ=255) aataatGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACt < 5:41032/149‑1 (MQ=255) atGCCGCAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAg < 7:168901/149‑1 (MQ=255) aCAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGTCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCg > 4:100322/1‑149 (MQ=255) ggAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTg > 4:167332/1‑133 (MQ=255) ggAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTg < 3:167332/133‑1 (MQ=255) gAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGc > 6:21848/1‑147 (MQ=255) cagcagG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGtt < 2:80882/149‑1 (MQ=255) tGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTc > 6:121161/1‑149 (MQ=255) cGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGCCatc < 6:25612/76‑3 (MQ=38) cGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGCCatc > 5:25612/1‑74 (MQ=38) cATCCCACACCCACCAGCTTCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCAGGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCg < 6:212986/149‑1 (MQ=255) cATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCg > 8:75812/1‑149 (MQ=255) cATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCg < 5:257561/149‑1 (MQ=255) caccaGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCcag > 8:26374/1‑149 (MQ=255) caGCTGCCCCAGGTCGGCTTGCCCCAGAAGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCcagcag > 1:151941/1‑149 (MQ=255) cccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg > 3:245078/1‑149 (MQ=255) ccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg > 3:201356/1‑148 (MQ=255) ccAGGTCGACTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg > 3:51330/1‑148 (MQ=255) gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcgc < 7:26374/147‑1 (MQ=255) ccccAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAg > 7:172831/1‑149 (MQ=255) cccAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGGCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTGATCTTCCACACCAGCCCCACCACCCCGGCCACAGCCAGCAGCACAAAGCAGGACTGCGCCAg > 1:215403/1‑148 (MQ=255) ggTTACCGGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGcc > 7:36483/1‑133 (MQ=255) ggTTACCGGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGcc < 8:36483/133‑1 (MQ=255) ggTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGca < 4:201356/149‑1 (MQ=255) ccAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATgg < 4:51330/149‑1 (MQ=255) aGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGa > 6:208346/1‑148 (MQ=255) | GCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGA > NC_002947/4915612‑4915903 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |