| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,193,337 | C→A | 44.4% | intergenic (+156/+639) | yqiK → / ← rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,193,337 | 0 | C | A | 44.4% | 23.3 / 24.0 | 27 | intergenic (+156/+639) | yqiK/rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (10/5); new base A (6/6); total (16/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.52e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01 | |||||||||||
TTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTC > W3110S.gb/3193276‑3193388 | ttttGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCt < 1:1636041/65‑1 (MQ=255) tGCGTGACGGCTCGGGGAATGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATcaacaa > 1:2658352/1‑71 (MQ=255) tGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCgg > 1:2446020/1‑71 (MQ=255) aCGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCggg > 1:459500/1‑70 (MQ=255) aCGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCggg > 1:1909818/1‑70 (MQ=255) cGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGa < 1:2358775/71‑1 (MQ=255) ggggAAGGCAAAAGTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagc < 1:2677162/71‑1 (MQ=255) ggggAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCGGGGAGTcagc < 1:566859/71‑1 (MQ=255) ggCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAggggg > 1:2667692/1‑71 (MQ=255) ggCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAggggg > 1:2554332/1‑71 (MQ=255) gCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGa < 1:2361946/58‑1 (MQ=255) cAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCt < 1:1302067/41‑1 (MQ=255) aaaaCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagcag > 1:706396/1‑64 (MQ=255) aaaCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagcag < 1:2202101/63‑1 (MQ=255) aaaCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagcag < 1:2042922/63‑1 (MQ=255) aaaCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCGGGGAGTcagc < 1:1145810/61‑1 (MQ=255) aCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAAt > 1:2026107/1‑47 (MQ=255) aCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagcag > 1:386120/1‑61 (MQ=255) aCTGCATCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATc > 1:999169/1‑48 (MQ=255) cTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAggggg > 1:2111405/1‑64 (MQ=255) gCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAAt > 1:336183/1‑44 (MQ=255) tCAGAGAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGaa < 1:2511175/71‑1 (MQ=255) agagAAAGCATTGCGTGGGGTCTGGAGCTCATcaacaa > 1:2358373/1‑38 (MQ=255) agagAAAGCATTGCGTGGGGTCTGGAGCTCATcaacaa > 1:2153737/1‑38 (MQ=255) agagAAAGCATTGCGTGGGGTCTGGAGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAg < 1:2789621/71‑1 (MQ=255) gaAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTcagca > 1:199351/1‑49 (MQ=255) aGCATTGCGTGGGGTCCGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTc > 1:1513874/1‑71 (MQ=255) | TTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTC > W3110S.gb/3193276‑3193388 |
| 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: ‑ |