| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 95,912 | T→C | 19.3% | G43G (GGT→GGC) | 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,912 | 0 | T | C | 19.3% | 47.1 / 5.4 | 36 | G43G (GGT→GGC) | mutT | nucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP |
| Reads supporting (aligned to +/- strand): ref base T (4/25); new base C (4/3); total (8/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.01e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01 | |||||||||||
CAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACC > minE/95854‑95966 | cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAGTTGAAATGGGCGAAACGCCgg < 1:2901013/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:1239362/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:3170618/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:3121630/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:1410600/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:641100/69‑1 (MQ=255)cagcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:749028/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:388092/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:2586467/68‑1 (MQ=255) agcagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCgg < 1:101280/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) gggTGAAACGCCTGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCgg < 1:3396749/51‑1 (MQ=255) | CAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACC > minE/95854‑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: ‑ |