Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,592,341 | T→G | 10.3% | L185R (CTG→CGG) ‡ | mccA → | methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,592,341 | 0 | T | G | 10.3% | 88.1 / 5.3 | 39 | L185R (CTG→CGG) ‡ | mccA | methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha |
Reads supporting (aligned to +/- strand): ref base T (17/18); new base G (2/2); total (19/20) | |||||||||||
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.71e-01 |
TACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGT > NC_002947/4592199‑4592486 | tACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGTGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCg > 1:589454/1‑150 (MQ=255) gAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGt < 1:277225/150‑1 (MQ=255) ggACCGGGACACCGTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGcaca < 2:229313/150‑1 (MQ=255) cTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAgtcg > 1:314315/1‑150 (MQ=255) acCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTg < 2:8246/117‑1 (MQ=255) acCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTg > 1:8246/1‑117 (MQ=255) gcGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGc < 2:213403/149‑1 (MQ=255) gAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtg > 2:92697/1‑150 (MQ=255) gCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggtgg < 2:589454/150‑1 (MQ=255) aTCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGaa > 1:443199/1‑150 (MQ=255) tACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATAAAAGTGGTCGAGGAACAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTAc < 2:531680/150‑1 (MQ=255) aCCCGGTGCTGCTCGAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACg > 1:228677/1‑150 (MQ=255) cccGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGt > 1:572896/1‑150 (MQ=255) ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGc < 1:564551/111‑1 (MQ=255) ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGc > 2:564551/1‑111 (MQ=255) ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg > 2:340745/1‑123 (MQ=255) ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg < 1:340745/123‑1 (MQ=255) tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt < 1:556920/87‑5 (MQ=37) tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt > 2:556920/1‑83 (MQ=37) tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAgccac < 1:92697/150‑1 (MQ=255) cAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGCGGGCTAGAGCCAAATGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGt > 2:21578/1‑150 (MQ=255) aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGcca > 2:190391/1‑75 (MQ=37) aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGcca < 1:190391/76‑2 (MQ=37) aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTg < 2:228677/150‑1 (MQ=255) ggCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg < 1:20032/150‑1 (MQ=255) gTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGgtgt > 2:75045/1‑150 (MQ=255) gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggt > 2:272500/1‑97 (MQ=255) gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggt < 1:272500/97‑1 (MQ=255) gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg > 1:591102/1‑137 (MQ=255) gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg < 2:591102/137‑1 (MQ=255) gAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGc > 2:555287/1‑150 (MQ=255) gAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGc > 1:129626/1‑150 (MQ=255) ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg < 1:516953/121‑1 (MQ=255) ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg > 2:516953/1‑121 (MQ=255) gaaAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTAc < 1:555287/150‑1 (MQ=255) tGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCaacgaac < 2:491639/150‑1 (MQ=255) gCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGt < 2:572896/150‑1 (MQ=255) cccTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGgaga > 1:331670/1‑93 (MQ=255) cccTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGgaga < 2:331670/93‑1 (MQ=255) | TACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGT > NC_002947/4592199‑4592486 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |