| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 226,884 | T→A | 8.1% | M327K (ATG→AAG) | frsA → | hydrolase, binds to enzyme IIA(Glc) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 226,884 | 0 | T | A | 8.1% | 170.6 / 4.9 | 60 | M327K (ATG→AAG) | frsA | hydrolase, binds to enzyme IIA(Glc) |
| Reads supporting (aligned to +/- strand): ref base T (33/22); new base A (2/3); total (35/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTG > minE/226819‑226951 | tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatg < 1:513269/71‑1 (MQ=255)tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCat > 1:124862/1‑70 (MQ=255)tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCat > 1:3640510/1‑70 (MQ=255) aTTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat > 1:3055332/1‑71 (MQ=255) ttttAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat > 1:2625742/1‑70 (MQ=255) tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa > 1:3584917/1‑71 (MQ=255) tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa < 1:2590749/71‑1 (MQ=255) gCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGAt > 1:995053/1‑71 (MQ=255) aGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGaa < 1:1438712/71‑1 (MQ=255) gCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAg > 1:2315251/1‑71 (MQ=255) gCAACACGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAg > 1:3582979/1‑71 (MQ=255) aaCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg > 1:2155474/1‑71 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATg > 1:3080399/1‑65 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGa > 1:3386003/1‑66 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGa > 1:3642173/1‑66 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt > 1:1931874/1‑71 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt > 1:3575954/1‑71 (MQ=255) aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt > 1:2315821/1‑71 (MQ=255) aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc < 1:1588636/59‑1 (MQ=38) aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc < 1:566931/59‑1 (MQ=38) aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc < 1:439593/59‑1 (MQ=38) tGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTgcgc < 1:1669603/71‑1 (MQ=255) ccGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt < 1:2830726/71‑1 (MQ=255) cGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTg < 1:26798/71‑1 (MQ=255) ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg < 1:3401596/61‑1 (MQ=255) ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg < 1:1379236/61‑1 (MQ=255) ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg < 1:2792925/61‑1 (MQ=255) tGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTg < 1:3258260/71‑1 (MQ=255) tATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt < 1:805079/62‑1 (MQ=255) tATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGaa < 1:2627397/71‑1 (MQ=255) aTCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt < 1:1546280/61‑1 (MQ=255) tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg > 1:1829862/1‑62 (MQ=255) tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg > 1:1456278/1‑62 (MQ=255) cGTTCTGGCCAGTCGTTTGGGGATGCATGATGATTCCGATGAAGCGt < 1:2527099/47‑1 (MQ=255) gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc < 1:1050072/71‑1 (MQ=255) gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc < 1:2502745/71‑1 (MQ=255) ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg < 1:2814034/44‑1 (MQ=255) ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg < 1:3092831/44‑1 (MQ=255) cTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGc > 1:911384/1‑55 (MQ=255) tccAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa < 1:2131884/68‑1 (MQ=255) gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa < 1:2772505/70‑1 (MQ=255) gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa < 1:2335458/70‑1 (MQ=255) gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa > 1:349469/1‑70 (MQ=255) gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa < 1:827774/70‑1 (MQ=255) aaGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATAGCTATTCATTAAAAg > 1:2134162/2‑71 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg > 1:1743050/1‑69 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:2212472/1‑71 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:880445/1‑71 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:2182875/1‑71 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:427046/1‑71 (MQ=255) gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCCTGGAg > 1:269663/1‑47 (MQ=255) tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:794893/1‑70 (MQ=255) tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:1706128/1‑70 (MQ=255) tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg > 1:513288/1‑70 (MQ=255) ttGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAg > 1:597724/1‑70 (MQ=255) tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg > 1:840908/1‑63 (MQ=255) tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg > 1:1657831/1‑63 (MQ=255) tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGa > 1:1645794/1‑71 (MQ=255) ggAAGCATCATGCATCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTg > 1:2191541/1‑71 (MQ=255) ggAAGCATCATGCATCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTg > 1:973165/1‑71 (MQ=255) | TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTG > minE/226819‑226951 |
| 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: ‑ |