| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,327,153 | C→T | 12.1% | intergenic (‑143/+133) | wcaK ← / ← wzxC | predicted pyruvyl transferase/colanic acid exporter |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,327,153 | 0 | C | T | 12.1% | 90.3 / 3.7 | 33 | intergenic (‑143/+133) | wcaK/wzxC | predicted pyruvyl transferase/colanic acid exporter |
| Reads supporting (aligned to +/- strand): ref base C (15/14); new base T (2/2); total (17/16) | |||||||||||
| 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.75e-01 | |||||||||||
GCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAC > minE/1327084‑1327221 | gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa < 1:3342320/71‑1 (MQ=255)gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa > 1:2991447/1‑71 (MQ=255)gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa > 1:2346680/1‑71 (MQ=255) cTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCAt < 1:1629873/71‑1 (MQ=255) gACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATcc > 1:1040995/1‑71 (MQ=255) aCGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCg < 1:2459691/71‑1 (MQ=255) tCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTg > 1:3359663/1‑71 (MQ=255) gCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTg < 1:2173838/71‑1 (MQ=255) ccTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa < 1:2513445/71‑1 (MQ=255) tACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATg > 1:2605930/1‑71 (MQ=255) cAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCg < 1:258520/71‑1 (MQ=255) cAAGACACGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATc > 1:6785/1‑50 (MQ=255) aaGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGa > 1:2317013/1‑71 (MQ=255) aGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc > 1:768160/1‑61 (MQ=255) aGAACCGTAGGACGGATAAGGCGTTCACGCAGCATCCGGCATTCAGTGCCTGATGCGa > 1:2046003/1‑58 (MQ=255) ccccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTa > 1:3443401/3‑71 (MQ=11) aaCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTa > 1:2399944/1‑71 (MQ=255) cccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTAt < 1:581783/70‑1 (MQ=11) ccGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATc < 1:1761952/71‑1 (MQ=255) gTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa < 1:3559753/47‑1 (MQ=255) gTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAg < 1:64397/71‑1 (MQ=255) tAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAgg > 1:320066/1‑71 (MQ=255) caGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTc < 1:392608/62‑1 (MQ=9) caGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg > 1:1780504/2‑71 (MQ=14) ggACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc < 1:3437885/71‑1 (MQ=255) gACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc > 1:3628021/1‑70 (MQ=11) aCGGATAAGGCTTTCACGCCTTATCCGACAACAACTTCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa < 1:3546791/71‑1 (MQ=2) aCGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa < 1:809116/71‑1 (MQ=255) aCGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa < 1:1623898/71‑1 (MQ=255) aTAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAg > 1:3711482/1‑70 (MQ=255) aaGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTc > 1:574995/1‑71 (MQ=255) ggCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTccc < 1:1327311/71‑1 (MQ=255) gTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAg > 1:139580/1‑71 (MQ=255) cgcATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAc > 1:3232511/1‑71 (MQ=255) gcATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc < 1:911362/52‑1 (MQ=255) | GCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAC > minE/1327084‑1327221 |
| 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: ‑ |