| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,586,895 | 1 | . | G | 58.8% | 32.6 / 22.2 | 17 | coding (836/1590 nt) | hsdM | DNA methylase M |
| Reads supporting (aligned to +/- strand): ref base . (6/1); new base G (0/10); total (6/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.66e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCG > W3110S.gb/4586841‑4586942 | ataataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGt < 1:78663/64‑1 (MQ=255) taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGtt > 1:2720996/1‑64 (MQ=255) taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTg > 1:1055198/1‑65 (MQ=255) taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTg > 1:884591/1‑65 (MQ=255) taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTg > 1:510831/1‑65 (MQ=255) taataTGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTg > 1:1888321/1‑65 (MQ=255) taataTGCTGCATAAAGCACAACTGATTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTg > 1:1129319/1‑65 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1090887/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:933632/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:2382227/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:2010812/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1780590/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1740519/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1568743/65‑1 (MQ=255) cGGGTGAACAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1477819/65‑1 (MQ=255) cGGGTGAAAAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1642154/65‑1 (MQ=255) aaCAAAGGTGCGGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCg < 1:1796975/59‑1 (MQ=255) | ATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGC‑GGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCG > W3110S.gb/4586841‑4586942 |
| 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: ‑ |