| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,535,985 | 1 | . | G | 65.0% | 43.4 / 21.3 | 20 | coding (248/846 nt) | nhoA | N‑hydroxyarylamine O‑acetyltransferase |
| Reads supporting (aligned to +/- strand): ref base . (5/2); new base G (13/0); total (18/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAG > W3110S.gb/1535936‑1536038 | ggTTACTGTTTTGAGCAGAATGGCGTGTTTTAGCGGGTGTTACGCGAGCT‑GGGGTTTAACg < 1:328185/61‑1 (MQ=255)ggTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACg < 1:3749751/61‑1 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGCGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:2308554/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:1246843/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:541718/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:4256210/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:4103522/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:4103253/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:391115/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:3421925/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:30896/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:2996662/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:279796/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:2167654/1‑62 (MQ=255) ttACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGGTTTAACGtt > 1:1547738/1‑62 (MQ=255) tACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCt < 1:365254/47‑1 (MQ=255) tACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCt < 1:3971728/47‑1 (MQ=255) tACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCt < 1:1295849/47‑1 (MQ=255) tACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGGGTTACGCGAGCt < 1:484794/47‑1 (MQ=255) tGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTa > 1:4193233/1‑61 (MQ=255) gcgAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAg > 1:4113455/1‑60 (MQ=255) gcgAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAg > 1:1221032/1‑60 (MQ=255) gcgAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAg > 1:689413/1‑60 (MQ=255) gcgAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAg > 1:935087/1‑60 (MQ=255) | GGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAG > W3110S.gb/1535936‑1536038 |
| 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: ‑ |