| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 95,901 | A→G | 20.5% | I40V (ATT→GTT) | mutT → | nucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 95,901 | 0 | A | G | 20.5% | 44.2 / 5.6 | 34 | I40V (ATT→GTT) | mutT | nucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP |
| Reads supporting (aligned to +/- strand): ref base A (2/25); new base G (4/3); total (6/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.56e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
CTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACC > minE/95837‑95966 | cTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGa < 1:2994804/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCgg < 1:2901013/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:1239362/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:3170618/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:1410600/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:3121630/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:749028/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:641100/69‑1 (MQ=255) cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:700704/69‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCgg < 1:608040/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCgg < 1:2634986/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:101280/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:388092/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:2586467/68‑1 (MQ=255) tGCGCACATGTCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGc < 1:2489541/69‑1 (MQ=255) gcgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCg > 1:1565982/1‑69 (MQ=255) gcgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCg > 1:1382976/1‑69 (MQ=255) cgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGa < 1:2046269/52‑1 (MQ=255) cgcACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGa < 1:1302549/52‑1 (MQ=255) aTAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCGGAACAGGCGGTGGTGc > 1:696032/1‑63 (MQ=39) aTAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCGGAACAGGCGGTGGTGCGTGAAc > 1:666124/1‑69 (MQ=39) aTAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCGGAACAGGCGGTGGTGCGTGAAc > 1:632653/1‑69 (MQ=39) aTAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCGGAACAGGCGGTGGTGCGTGAAc > 1:3409005/1‑69 (MQ=39) ggcggTAAAATTGAAATGGGTGAAACGCCTGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCgg < 1:1607482/68‑1 (MQ=255) cggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGAt < 1:712490/69‑1 (MQ=255) cggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGAt < 1:172230/69‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:863378/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:1013270/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:2593476/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:2203269/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:2118647/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:1901469/45‑1 (MQ=255) gTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTg < 1:1430385/45‑1 (MQ=255) tAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTAc < 1:2247593/69‑1 (MQ=255) aaaaTTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTAcc > 1:3247248/1‑69 (MQ=255) aaaaTTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTAcc > 1:2596416/1‑69 (MQ=255) | CTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACC > minE/95837‑95966 |
| 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: ‑ |