Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,213,568 | C→T | 33.3% | intergenic (+24/+55) | rpoD → / ← ygjF | RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,213,568 | 0 | C | T | 33.3% | 28.1 / 21.7 | 27 | intergenic (+24/+55) | rpoD/ygjF | RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base C (6/12); new base T (3/6); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01 |
ACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAC > W3110S.gb/3213501‑3213637 | aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc > 1:46357/1‑70 (MQ=255) cTGAAGTGCTGCGTAGCTTCCTTGACGATTAATCGGTAGGCCGGATCAGGCGTTAtgccg > 1:130452/1‑60 (MQ=38) cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGTATCAGGCGTTAcgccg > 1:704431/1‑60 (MQ=255) cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCAc < 1:631274/70‑1 (MQ=255) aaGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAgg < 1:3888897/71‑1 (MQ=255) ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc < 1:208632/71‑1 (MQ=255) ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc < 1:2496611/71‑1 (MQ=255) ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc < 1:3096667/71‑1 (MQ=255) cGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTg > 1:2254215/1‑70 (MQ=255) tAGCTTCCTGGACGATTACTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCa > 1:4034317/1‑56 (MQ=255) gCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCTTTACGCCGCACCCGGCACTAGGCCCTCTGcaca < 1:480106/70‑1 (MQ=255) cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCAc < 1:1996497/54‑1 (MQ=255) tCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCg > 1:1667732/1‑48 (MQ=255) cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGc > 1:4117606/1‑70 (MQ=255) aCGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCt < 1:516411/71‑1 (MQ=255) aCGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCt < 1:3492871/71‑1 (MQ=255) aTTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt < 1:3686332/71‑1 (MQ=255) ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc > 1:4252673/1‑51 (MQ=255) aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg > 1:4489749/1‑71 (MQ=255) aTCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg < 1:2238486/70‑1 (MQ=255) aTCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg < 1:350027/70‑1 (MQ=255) aTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg < 1:4204658/70‑1 (MQ=255) aTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCgg < 1:2692453/70‑1 (MQ=255) ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGttttt > 1:272940/1‑71 (MQ=255) cGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttttt < 1:3723022/71‑1 (MQ=255) aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGc < 1:3540886/70‑1 (MQ=255) ttATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGc < 1:4146577/71‑1 (MQ=255) ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGc < 1:2916615/71‑1 (MQ=255) aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc > 1:296779/1‑71 (MQ=255) aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc > 1:356418/1‑71 (MQ=255) cgccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc < 1:683913/70‑1 (MQ=255) | ACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAC > W3110S.gb/3213501‑3213637 |
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: ‑ |