| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,289,607 | 1 | . | C | 73.3% | 76.6 / 25.2 | 30 | coding (319/2652 nt) | ppc | phosphoenolpyruvate carboxylase |
| Reads supporting (aligned to +/- strand): ref base . (0/8); new base C (22/0); total (22/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.71e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACAGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACT > minE/2289553‑2289638 | cacaGTATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:2617998/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCCGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1181455/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACt > 1:605335/1‑69 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:2947078/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:864460/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:725779/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:541332/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:479539/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:462797/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:391326/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:3255308/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:3248936/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1162776/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:264760/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:2578636/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:2058857/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1963189/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:185021/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1643678/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1632046/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCTGCGTAAACTg > 1:1174544/1‑70 (MQ=255)cacaGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCACCCCGGCGTAAACTg > 1:2950697/1‑70 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:2257642/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:3224587/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:3235770/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:446458/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:1619822/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:1499068/41‑1 (MQ=255) aTCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:1451492/41‑1 (MQ=255) tCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACt < 1:545759/40‑1 (MQ=255) | CACAGCATTTCGCCGAAAGGCGAAGCTGCCAGCAACCCGGAAGTGATCGCCCGCA‑CCCTGCGTAAACTGAAAAACCAGCCGGAACT > minE/2289553‑2289638 |
| 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: ‑ |