| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 186,630 | C→T | 41.7% | H894Y (CAT→TAT) | dnaE → | DNA polymerase III alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 186,630 | 0 | C | T | 41.7% | 16.0 / 40.1 | 36 | H894Y (CAT→TAT) | dnaE | DNA polymerase III alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (10/11); new base T (5/10); total (15/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.01e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01 | |||||||||||
CAAAAAGTTGAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGCG > minE/186563‑186695 | cAAAAAGTTGAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACAt > 1:1806775/1‑70 (MQ=255) aaaaaGTTGAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCATATc > 1:623178/1‑70 (MQ=255) aGTTGAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCa < 1:3267578/71‑1 (MQ=255) gTTGAACCGTCGCGTGCTGGAAAAACTTATCATGTCCGGGGCGTTTGACCGTCTTGGGCCATATCGCGCAg < 1:3798075/71‑1 (MQ=255) gAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCATATCGCGCAgcg > 1:2133385/1‑70 (MQ=255) aaCCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCt < 1:308347/71‑1 (MQ=255) aCCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCt < 1:2621784/70‑1 (MQ=255) tGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAAc < 1:3774144/68‑1 (MQ=255) gAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCcac < 1:976899/43‑1 (MQ=255) aaaaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCATATCGCGCAGCGCtgatg < 1:3572544/59‑1 (MQ=255) aaaaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGc > 1:4046490/1‑71 (MQ=255) aaaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGc < 1:3161049/70‑1 (MQ=255) aaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGa > 1:3440485/1‑71 (MQ=255) aaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGa > 1:2484426/1‑71 (MQ=255) aaaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGCCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGa > 1:1440241/1‑71 (MQ=255) aaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTggg > 1:2313908/1‑67 (MQ=255) aaCTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGAt < 1:747570/71‑1 (MQ=255) aCTGATCATGTACGGGGCGTTTGACCGTCTTGGGCCACATCGCGCa > 1:505444/1‑46 (MQ=255) gATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCgcgc < 1:3563431/42‑1 (MQ=255) aTGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTaaa < 1:3850894/71‑1 (MQ=255) gCGTTTGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGa < 1:2763669/71‑1 (MQ=255) gTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGAt > 1:1701065/1‑70 (MQ=255) tttGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGc < 1:1802806/48‑1 (MQ=255) tttGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGc < 1:3817019/48‑1 (MQ=255) tGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGc < 1:259768/46‑1 (MQ=255) tGACCGTCTTGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCa > 1:3914109/1‑64 (MQ=255) tGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCaa > 1:1495324/1‑70 (MQ=255) ccGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACAcg < 1:2767069/71‑1 (MQ=255) ttGGGCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATc < 1:852955/60‑1 (MQ=255) gggCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAg < 1:1855615/70‑1 (MQ=255) gggCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGc > 1:3314798/1‑71 (MQ=255) gggCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAg < 1:1145984/70‑1 (MQ=255) gggCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGc > 1:2137022/1‑71 (MQ=255) ggcccACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAg < 1:2438148/67‑1 (MQ=255) ggCCATATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGaaa < 1:2963667/68‑1 (MQ=255) ggCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGCg > 1:618808/1‑71 (MQ=255) ggCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGCg < 1:3148624/71‑1 (MQ=255) acaTCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGaa < 1:250650/63‑1 (MQ=255) | CAAAAAGTTGAACCGTCGCGTGCTGGAAAAACTGATCATGTCCGGGGCGTTTGACCGTCTTGGGCCACATCGCGCAGCGCTGATGAACTCGCTGGGCGATGCGTTAAAAGCGGCAGATCAACACGCGAAAGCG > minE/186563‑186695 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |