Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 405,632 | A→G | 19.4% | T2A (ACA→GCA) | aroL → | shikimate kinase II |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 405,632 | 0 | A | G | 19.4% | 74.4 / 13.0 | 36 | T2A (ACA→GCA) | aroL | shikimate kinase II |
Reads supporting (aligned to +/- strand): ref base A (16/13); new base G (5/2); total (21/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATT > W3110S.gb/405577‑405701 | ttCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGgcac < 1:2524007/59‑1 (MQ=255) ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACaa > 1:221117/1‑59 (MQ=255) ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCttttt > 1:3487738/1‑68 (MQ=255) ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCt < 1:2788388/70‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTg < 1:2485056/68‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTg < 1:1511783/68‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg < 1:3391263/68‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg < 1:3169897/68‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg < 1:1664992/68‑1 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg < 1:821638/68‑1 (MQ=255) cTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGgca > 1:3118074/1‑43 (MQ=255) tCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGc < 1:1633181/70‑1 (MQ=255) gcATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg < 1:1000915/50‑1 (MQ=255) aTTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:1987805/71‑1 (MQ=255) aTTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:513987/71‑1 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCttttt > 1:3087748/1‑42 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGa > 1:310985/1‑46 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt > 1:1490487/1‑47 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt > 1:2372418/1‑47 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt > 1:855786/1‑47 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt > 1:2536503/1‑47 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa > 1:1090635/1‑70 (MQ=255) gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa > 1:244907/1‑70 (MQ=255) aTTGGGGAAAACCCACGATGGCACAACCTCttttt > 1:2393117/1‑35 (MQ=255) aTTGGGGAAAACCCACGATGGCACAACCTCTTTTTCTGAt > 1:2867570/1‑40 (MQ=255) ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg < 1:2811872/70‑1 (MQ=255) ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg < 1:2425266/70‑1 (MQ=255) ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg > 1:1298047/1‑70 (MQ=255) gAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATg > 1:3629263/1‑71 (MQ=255) aaaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCtg > 1:2548552/1‑48 (MQ=255) aaaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:711388/54‑1 (MQ=255) aaCCCACGATGGCACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGcc > 1:939197/1‑71 (MQ=255) aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa > 1:1994743/1‑54 (MQ=255) aCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGccc < 1:974577/71‑1 (MQ=255) cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATgg > 1:2832447/1‑67 (MQ=255) cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGc > 1:1913676/1‑68 (MQ=255) ccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCtt > 1:3068813/1‑71 (MQ=255) gACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGAtt > 1:415551/1‑71 (MQ=255) acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:1482170/41‑1 (MQ=255) acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:632895/41‑1 (MQ=255) acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:404688/41‑1 (MQ=255) acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa < 1:1224598/41‑1 (MQ=255) acacAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGaa > 1:2927815/1‑55 (MQ=255) | TTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATT > W3110S.gb/405577‑405701 |
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: ‑ |