Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,915,752 | A→T | 12.1% | V112V (GTT→GTA) | 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,752 | 0 | A | T | 12.1% | 80.2 / 2.1 | 31 | V112V (GTT→GTA) | ccmC | protoheme IX reservoir complex subunit |
Reads supporting (aligned to +/- strand): ref base A (17/10); new base T (2/2); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCA > NC_002947/4915612‑4915897 | 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 < 3:167332/133‑1 (MQ=255) ggAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTg > 4:167332/1‑133 (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 > 5:25612/1‑74 (MQ=38) cGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGCCatc < 6:25612/76‑3 (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) | GCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGG‑ATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCA > NC_002947/4915612‑4915897 |
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: ‑ |