| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,886,161 | 0 | T | G | 56.5% | 19.0 / 22.8 | 23 | intergenic (‑286/+73) | ygcB/cysH | conserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase |
| Reads supporting (aligned to +/- strand): ref base T (7/3); new base G (13/0); total (20/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.78e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATT > W3110S.gb/2886106‑2886200 | gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAgg < 1:4538605/60‑1 (MQ=255)gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAgg < 1:2523744/60‑1 (MQ=255) ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAAtgct > 1:98435/1‑62 (MQ=255) aCGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATgctgct > 1:172757/1‑62 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAAtgctgat > 1:1797400/2‑59 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:1648787/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:1062876/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:4590883/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:1103721/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:4149187/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:3156590/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:1133035/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:2195486/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:2109773/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:2048109/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:1924695/2‑61 (MQ=255) tgCCATTCTATACCTGTTAATTCTCCCAGCGACGTTCTCGTTTAAATGCCGGAATgctgct > 1:3966204/2‑61 (MQ=255) aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa > 1:3485424/1‑62 (MQ=255) aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa > 1:4529951/1‑62 (MQ=255) aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa > 1:725551/1‑62 (MQ=255) ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaa > 1:3419283/1‑46 (MQ=255) ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaa > 1:1549579/1‑46 (MQ=255) cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCAtt < 1:907924/60‑1 (MQ=255) | GCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATT > W3110S.gb/2886106‑2886200 |
| 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: ‑ |