| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,336,585 | 0 | C | A | 20.0% | 101.6 / 13.6 | 40 | C24* (TGC→TGA) | fdoI | formate dehydrogenase‑O, cytochrome b556 subunit |
| Reads supporting (aligned to +/- strand): ref base C (13/19); new base A (8/0); total (21/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.63e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACGTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAACTGGTTG > minE/2336522‑2336642 | acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:3188966/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:289681/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:593750/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:707112/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:3030222/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:748238/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:1673093/1‑67 (MQ=255)acgTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCAGATTcgg > 1:334314/1‑67 (MQ=255) tGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTgg > 1:1924171/1‑71 (MQ=255) acacCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTggcg < 1:1425443/71‑1 (MQ=255) cacCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTggcgg < 1:159521/71‑1 (MQ=255) cacCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTggcgg < 1:349313/71‑1 (MQ=255) cacCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTggcgg < 1:2474017/71‑1 (MQ=255) ccATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGgcggcg < 1:25682/71‑1 (MQ=255) cATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGcggcgg > 1:1646130/1‑71 (MQ=255) aTCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGt < 1:2417919/71‑1 (MQ=255) aTCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGt < 1:54712/71‑1 (MQ=255) tCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGt > 1:1328055/1‑70 (MQ=255) cgcTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCg < 1:2667207/70‑1 (MQ=255) tACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggct > 1:338923/1‑71 (MQ=255) tACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggct > 1:2422430/1‑71 (MQ=255) tACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggct > 1:1511882/1‑71 (MQ=255) tACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggc < 1:764203/70‑1 (MQ=255) tACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggc < 1:11247/70‑1 (MQ=255) acacCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCtt > 1:359078/1‑46 (MQ=255) acacCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTc > 1:2081543/1‑47 (MQ=255) acacCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTc > 1:964030/1‑47 (MQ=255) cgcgCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCgggctgggct > 1:1334331/1‑71 (MQ=255) cgcCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCttt < 1:3096142/71‑1 (MQ=255) cgcCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCttt < 1:1005519/71‑1 (MQ=255) gAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGtt < 1:2527888/71‑1 (MQ=255) cGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGtt < 1:1743054/68‑1 (MQ=255) tCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGc < 1:985769/47‑1 (MQ=255) ccACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTc > 1:169574/1‑70 (MQ=255) cTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAACt < 1:2679626/71‑1 (MQ=255) tGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTc < 1:1258660/66‑1 (MQ=255) aTCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAACTGGtt < 1:1911230/71‑1 (MQ=255) tCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAACTGGTTg < 1:1076659/71‑1 (MQ=255) cTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCtt > 1:2981605/1‑36 (MQ=255) cTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAAc > 1:2085329/1‑55 (MQ=255) | ACGTGACACCATCGTGCGCTACACCGCGCCGGAACGTATCAACCACTGGATCACCGCCTTCTGCTTCATCCTGGCGGCGGTGAGCGGGCTGGGCTTTTTGTTCCCGTCCTTCAACTGGTTG > minE/2336522‑2336642 |
| 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: ‑ |