| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,951,645 | T→C | 16.7% | Q456Q (CAA→CAG) | ppk ← | polyphosphate kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,951,645 | 0 | T | C | 16.7% | 49.6 / 4.4 | 24 | Q456Q (CAA→CAG) | ppk | polyphosphate kinase |
| Reads supporting (aligned to +/- strand): ref base T (12/8); new base C (2/2); total (14/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01 | |||||||||||
TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCC > NC_002947/5951497‑5951771 | tgtACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGGCTTGAAGCCAACTACACCGTAGATCACTACCGCACCAGCCGCtt < 1:114606/150‑1 (MQ=255) gCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGc > 1:392090/1‑150 (MQ=255) cGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCa < 2:377316/150‑1 (MQ=255) ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaa > 1:494319/1‑150 (MQ=255) gTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtg < 2:505045/115‑3 (MQ=255) gTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtg > 1:505045/1‑113 (MQ=255) cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTg < 2:409079/121‑1 (MQ=255) cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTg > 1:409079/1‑121 (MQ=255) aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCg > 2:248476/1‑150 (MQ=255) ctcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAc < 2:98190/150‑1 (MQ=255) gcggcgTGGAATCATCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGATTCTTCGTCAAAACGCTCACGCAACTCGATCACCGCAGTGCcctc < 1:26292/150‑1 (MQ=255) ggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctcct > 1:209834/1‑150 (MQ=255) gAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGt > 1:273500/1‑150 (MQ=255) aGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCCCACGCAACTCGATCACCGACGTGACCTCCTTGCCGTTAcgc > 2:59144/1‑150 (MQ=255) gCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGt > 1:43263/1‑87 (MQ=255) gCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGt < 2:43263/87‑1 (MQ=255) tcatcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTg < 2:364098/79‑1 (MQ=255) tcatcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTg > 1:364098/1‑79 (MQ=255) tGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAgg > 2:51049/1‑150 (MQ=255) ggTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGc < 2:494319/150‑1 (MQ=255) ttGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcg > 2:101085/1‑150 (MQ=255) cacaCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa < 2:392090/150‑1 (MQ=255) tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGcc > 1:244336/1‑150 (MQ=255) cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGAtt < 1:78125/122‑1 (MQ=255) cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGAtt > 2:78125/1‑122 (MQ=255) | TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCC > NC_002947/5951497‑5951771 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |