| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,246,521 | C→T | 33.8% | G78G (GGC→GGT) | thiC → | thiamin (pyrimidine moiety) biosynthesis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,246,521 | 0 | C | T | 33.8% | 45.5 / 28.8 | 59 | G78G (GGC→GGT) | thiC | thiamin (pyrimidine moiety) biosynthesis protein |
| Reads supporting (aligned to +/- strand): ref base C (33/6); new base T (8/12); total (41/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.05e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTA > minE/2246456‑2246575 | ggCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCg < 1:556623/69‑1 (MQ=255) aGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGc > 1:356174/1‑60 (MQ=255) aGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGc > 1:3039100/1‑60 (MQ=255) aGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGt > 1:447449/1‑69 (MQ=255) aaGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCTGTCTACGACACCTCCGGCCCGTATGGTGATc < 1:1463555/69‑1 (MQ=255) aaGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc < 1:1506843/69‑1 (MQ=255) aGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc < 1:1668164/68‑1 (MQ=255) aGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc < 1:3427070/68‑1 (MQ=255) aGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc < 1:2594984/68‑1 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCg > 1:3356496/1‑56 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:3096188/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:1307383/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:2999168/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:3387573/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:2887821/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:2737740/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:1283786/1‑66 (MQ=255) aaCAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATc > 1:3311238/1‑66 (MQ=255) gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg > 1:2248704/1‑68 (MQ=255) gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg > 1:2683070/1‑68 (MQ=255) gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg > 1:3337803/1‑68 (MQ=255) gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg > 1:1569090/1‑68 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:207332/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:3140486/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:1291207/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:1079652/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:1330964/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:3398720/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:171432/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:2627979/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:1928998/69‑1 (MQ=255) aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:2141788/69‑1 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:1771395/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:2671465/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:1773763/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:1211756/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:3079204/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:1535707/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:232537/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa > 1:1510760/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:1177248/68‑1 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa < 1:534856/68‑1 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCa > 1:2710793/1‑57 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa > 1:837904/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa > 1:298021/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa > 1:2810329/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa > 1:280636/1‑69 (MQ=255) cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa > 1:1939589/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:858596/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:2567655/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:2448481/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:1841543/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:1391099/1‑69 (MQ=255) gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg > 1:1319477/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:3099152/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:3349629/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:3097816/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:2982035/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:2690262/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:262835/1‑69 (MQ=255) tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa > 1:1633548/1‑69 (MQ=255) | GGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTA > minE/2246456‑2246575 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |