Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 2,865,292 | G→A | 19.6% | Q267* (CAG→TAG) | rpoS ← | RNA polymerase, sigma S (sigma 38) factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,865,292 | 0 | G | A | 19.6% | 123.6 / 23.7 | 56 | Q267* (CAG→TAG) | rpoS | RNA polymerase, sigma S (sigma 38) factor |
Reads supporting (aligned to +/- strand): ref base G (26/19); new base A (4/7); total (30/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.13e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01 |
GCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC > W3110S.gb/2865229‑2865362 | gcgcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg > 1:4375915/1‑70 (MQ=255) gcTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg < 1:1527983/71‑1 (MQ=255) gATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAacac < 1:2524628/71‑1 (MQ=255) atatTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg < 1:3113269/71‑1 (MQ=255) aGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCAc < 1:1767045/70‑1 (MQ=255) gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg > 1:1083917/1‑71 (MQ=255) cTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa < 1:377133/71‑1 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGAAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg > 1:1777521/1‑71 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt > 1:2924433/1‑46 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg > 1:2217424/1‑71 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg > 1:304162/1‑70 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg > 1:1786782/1‑70 (MQ=255) cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa > 1:3918867/1‑61 (MQ=255) cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa > 1:1539130/1‑61 (MQ=255) cGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa > 1:3596099/1‑61 (MQ=255) tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa < 1:450880/71‑1 (MQ=255) tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGcc > 1:2666497/1‑70 (MQ=255) tttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAAAACGTTCACGGGTGAGGcc > 1:392943/1‑70 (MQ=255) gCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa > 1:1043569/1‑69 (MQ=255) cAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt > 1:114975/1‑71 (MQ=255) cAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt < 1:4043481/71‑1 (MQ=255) tgATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg < 1:1878271/56‑1 (MQ=255) gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc < 1:2141719/53‑1 (MQ=255) gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt < 1:1097314/57‑1 (MQ=255) gATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg < 1:234287/58‑1 (MQ=255) tttCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg > 1:595327/1‑52 (MQ=255) tCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa > 1:4237505/1‑48 (MQ=255) tCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc > 1:3873890/1‑71 (MQ=255) cgcgCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc > 1:3452134/1‑70 (MQ=255) cgcAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc < 1:2378390/71‑1 (MQ=255) gcAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:3875656/71‑1 (MQ=255) gcAAACGGCGCAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:3383390/71‑1 (MQ=255) aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAacac > 1:4363206/1‑36 (MQ=255) aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt < 1:55978/57‑1 (MQ=255) aaCGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCt > 1:3369568/1‑71 (MQ=255) aaCGGCGCAGGCCTTCAACCTGAATCTG‑‑GAACACGTTCACGGGTGAGGcc < 1:3198942/50‑1 (MQ=255) ggCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc < 1:4275449/64‑1 (MQ=255) ggCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc < 1:3960717/71‑1 (MQ=255) gcgcAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCa < 1:136746/71‑1 (MQ=255) gcAGGCCTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAgt < 1:3269529/71‑1 (MQ=255) gCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa > 1:1674942/1‑52 (MQ=255) ccTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc > 1:1395355/1‑71 (MQ=255) cTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc > 1:1320624/1‑71 (MQ=255) cTTCAACCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc > 1:231901/1‑62 (MQ=255) ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:3599599/47‑1 (MQ=255) ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:3131810/47‑1 (MQ=255) ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc < 1:4005885/61‑1 (MQ=255) tCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc > 1:3927816/1‑59 (MQ=255) tCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGtt < 1:1227117/66‑1 (MQ=255) cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa > 1:546218/1‑44 (MQ=255) cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt > 1:2002739/1‑71 (MQ=255) cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt > 1:1618751/1‑71 (MQ=255) cAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCt > 1:4534652/1‑71 (MQ=255) aaCCTAAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCtt < 1:1753398/71‑1 (MQ=255) aaCCTAAATCTGGCGAACACCTTCACGGGTGAGGCCAATTTCa > 1:532393/1‑43 (MQ=38) cTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAgt < 1:1535188/58‑1 (MQ=255) gAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCgatc > 1:2922923/1‑68 (MQ=255) gAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCgatc > 1:882879/1‑68 (MQ=255) | GCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAC > W3110S.gb/2865229‑2865362 |
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: ‑ |