| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 470,861 | 0 | C | A | 68.4% | 10.8 / 17.1 | 19 | P76Q (CCG→CAG) | sucD | succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (0/6); new base A (13/0); total (13/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.69e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATCGACGCAGGCATCAAACTG > minE/470822‑470910 | ccACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATc < 1:3138251/71‑1 (MQ=255)ccACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATc < 1:2788908/71‑1 (MQ=255)ccACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATc < 1:2738867/71‑1 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCAt > 1:1965354/1‑69 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCAt > 1:2834338/1‑69 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:1364762/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:426037/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:3066660/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:2845498/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:2405620/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:2386816/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:2090082/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:1613991/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:16071/1‑70 (MQ=255) cACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:1340053/1‑70 (MQ=255) cACTGGAGCTACCGCTTCTGTTATCTACGTACCAGCACAGTTCTGCAAAGACTCCATTCTGGAAGCCATc > 1:419749/1‑70 (MQ=255) cgcTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATc < 1:3199751/64‑1 (MQ=255) ttATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATCGACGCAGGCATCAAACTg < 1:2144992/68‑1 (MQ=255) ttATCTACGTACCAGCACCGTTCTGCAAAGACTACATTCTGGAAGCCATCGACGCAGGCATCAAACTg < 1:643168/68‑1 (MQ=255) | CCACTGGCGCTACCGCTTCTGTTATCTACGTACCAGCACCGTTCTGCAAAGACTCCATTCTGGAAGCCATCGACGCAGGCATCAAACTG > minE/470822‑470910 |
| 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: ‑ |