| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,363,111 | G→T | 11.7% | L235I (CTC→ATC) | glmM ← | phosphoglucosamine mutase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,363,111 | 0 | G | T | 11.7% | 83.5 / 3.1 | 35 | L235I (CTC→ATC) | glmM | phosphoglucosamine mutase |
| Reads supporting (aligned to +/- strand): ref base G (12/18); new base T (2/2); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01 | |||||||||||
CATCAGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGC > NC_002947/5362962‑5363247 | catcaGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCgag < 2:833981/150‑1 (MQ=255) atcaGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCgaga > 1:147964/1‑150 (MQ=255) atcaGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCgaga > 2:377056/1‑150 (MQ=255) cGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTc > 1:739784/1‑150 (MQ=255) ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAttgcctt < 1:711990/150‑7 (MQ=255) ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGa < 1:234562/150‑1 (MQ=255) ttGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGc > 1:242911/1‑150 (MQ=255) cTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTg < 1:377056/150‑1 (MQ=255) gacATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAg < 1:585313/148‑1 (MQ=255) tGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGAt < 1:529963/150‑1 (MQ=255) gCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt < 2:517027/150‑1 (MQ=255) tCGCGAGCAATGACGAACAGCACCTCGTCCCCGTCGATGATGGCGCCCGTGTGGTCGGACATCAGGCCGCGATCGCCGTCACCGTCAAAGGCAATGCCGAAATCGGCATGTCCGACCAGTACCGCGGCCTGCAGCGATTCGATGTGTGTc < 2:815088/150‑1 (MQ=255) gcgAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGTGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGa < 2:37331/150‑1 (MQ=255) gcgAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGa < 2:637848/150‑1 (MQ=255) tgatgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAAc < 2:739784/150‑1 (MQ=255) gaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTc > 2:575758/1‑150 (MQ=255) aaCAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCTACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCg > 2:694389/1‑150 (MQ=255) aCAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGt < 1:575758/150‑1 (MQ=255) gcagcTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGa < 2:242911/150‑1 (MQ=255) cTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGtt > 2:146238/1‑150 (MQ=255) cGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGGTCGGCATGTCGGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAATCTTCGTTGATGTTCa < 2:147964/150‑1 (MQ=255) aCCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAcc > 2:150317/1‑88 (MQ=255) aCCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAcc < 1:150317/88‑1 (MQ=255) gTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTc > 2:41308/1‑150 (MQ=255) ggCGCCCGCGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTcc < 1:221483/72‑1 (MQ=255) ggCGCCCGCGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTcc > 2:221483/1‑72 (MQ=255) cacggTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTgcgc > 2:13804/4‑142 (MQ=255) cacggTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTgcgc < 1:13804/139‑1 (MQ=255) tgGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACg < 1:41308/150‑1 (MQ=255) gCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCAc > 2:23426/1‑150 (MQ=255) gcgATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCAcaccg > 1:358007/1‑146 (MQ=255) gcgATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCaa > 2:348152/1‑149 (MQ=255) cGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCAc < 2:174270/150‑1 (MQ=255) gCCGTCACCGTCAAAGGCAATGCCGAGATCGGCACGTCCGATCTACCGGCCGTCTTCCAGCGATTCGATGTGTGTCGGGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACg < 1:405697/150‑1 (MQ=255) cgtcaccgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGa < 1:23426/150‑1 (MQ=255) tcggacaTGCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGgc > 2:682960/7‑150 (MQ=255) caTGCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAg < 2:151855/41‑1 (MQ=255) caTGCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAg > 1:151855/2‑42 (MQ=255) | CATCAGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGC > NC_002947/5362962‑5363247 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |