| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,639,432 | C→T | 33.9% | T119T (ACC→ACT) | ugpQ → | glycerophosphodiester phosphodiesterase, cytosolic |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,639,432 | 0 | C | T | 33.9% | 36.4 / 24.3 | 50 | T119T (ACC→ACT) | ugpQ | glycerophosphodiester phosphodiesterase, cytosolic |
| Reads supporting (aligned to +/- strand): ref base C (18/15); new base T (2/15); total (20/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.36e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01 | |||||||||||
GCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGG > minE/2639369‑2639491 | gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:876559/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:2753103/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:3232474/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:1526657/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:3281723/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:1604085/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:2104880/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:2029074/67‑1 (MQ=255)gCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACcaccggc < 1:1949893/67‑1 (MQ=255) tGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGc < 1:937724/69‑1 (MQ=255) tGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGc < 1:3190390/69‑1 (MQ=255) tGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGc < 1:2808542/69‑1 (MQ=255) tGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGc < 1:457214/69‑1 (MQ=255) tGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGc < 1:1858142/69‑1 (MQ=255) ggCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCgggg < 1:3043485/68‑2 (MQ=255) ggCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGcc > 1:600992/1‑69 (MQ=255) ggCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGcc > 1:3095273/1‑69 (MQ=255) ggCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGcc > 1:1141418/1‑69 (MQ=255) ggCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGcc > 1:3245103/1‑69 (MQ=255) gcgAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCaaa < 1:2910759/68‑1 (MQ=255) gcgAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCaaa < 1:2907558/68‑1 (MQ=255) gcgAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCaaa < 1:1644597/68‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCtaaa < 1:2451053/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:2353229/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:3097285/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1273041/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1011869/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1397476/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1876873/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:791859/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:383045/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:2213155/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:738804/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:175656/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:2011092/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1886570/66‑1 (MQ=255) gcaaaCACGGGATGATGGCGAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCaaaa < 1:1883534/66‑1 (MQ=255) gAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCaa > 1:267839/1‑67 (MQ=255) gAATATCGAAATCAAACCCACCACTGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCAc > 1:1972051/1‑68 (MQ=255) gAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAAt > 1:1124014/1‑51 (MQ=255) gAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCAc > 1:2032019/1‑68 (MQ=255) gAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCAc > 1:241725/1‑68 (MQ=255) gAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACAGGCAAAATGGTGGCGCTGGCGGCAc > 1:3183335/1‑68 (MQ=255) atatCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCAcgc > 1:1653892/1‑68 (MQ=255) atatCGAAATCAAACCCACCACCGGCAACGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCAcgcg > 1:285397/1‑69 (MQ=255) tCGAAATCAAAGCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:1136157/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGgcg > 1:2883640/1‑53 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:2633181/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:2279618/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:908866/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:3017442/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:2894076/1‑69 (MQ=255) tCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGc > 1:678729/1‑69 (MQ=255) aCCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:1540961/69‑1 (MQ=255) aCCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:2269945/69‑1 (MQ=255) aCCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:58772/69‑1 (MQ=255) aCCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:3361266/69‑1 (MQ=255) aCCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:2900645/69‑1 (MQ=255) cccACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCgg < 1:728819/68‑1 (MQ=255) | GCAGGTGGCGGAACGTTGCCGCGAACACGGGATGATGGCGAATATCGAAATCAAACCCACCACCGGCACCGGGCCATTAACGGGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGG > minE/2639369‑2639491 |
| 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: ‑ |