| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 618,746 | 0 | C | . | 57.1% | 74.4 / 119.0 | 49 | coding (53/1674 nt) | rpsA | 30S ribosomal subunit protein S1 |
| Reads supporting (aligned to +/- strand): ref base C (19/2); new base . (0/28); total (19/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.11e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCT > minE/618709‑618783 | gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCTGGTTCTATCGTTCGTg < 1:3117533/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:185695/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2316232/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2274152/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2373049/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2211665/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2196107/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2523075/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2839541/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2077395/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:301642/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2066809/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2064953/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:2057865/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:19339/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1828792/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:310715/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1663246/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1656958/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:3125137/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1463701/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1461883/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:129963/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1294823/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:545841/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:642274/60‑1 (MQ=255)gCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1003936/60‑1 (MQ=255) cAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA‑CCGCCCGGGTTCTATCGTTCGTg < 1:1036076/57‑1 (MQ=255) tGAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCgttgttgt < 1:2514789/60‑1 (MQ=255) tGAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCgttgttgt < 1:2888665/60‑1 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTgg > 1:912818/1‑40 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2071421/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:861796/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:751339/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:704882/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:1074790/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:1075303/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:366712/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:1513557/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:1706787/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:1750293/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2228596/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2906254/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2775234/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:275406/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2632773/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2161195/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2386601/1‑53 (MQ=255) tAAAAGAAATCGAAACCCGACCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCt > 1:2148922/1‑53 (MQ=255) | GCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCT > minE/618709‑618783 |
| 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: ‑ |