Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 911,817 | 0 | G | A | 56.2% | 0.3 / 16.8 | 16 | A252A (GCC→GCT) | hcr | HCP oxidoreductase, NADH‑dependent |
Reads supporting (aligned to +/- strand): ref base G (2/5); new base A (9/0); total (11/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.81e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAGACCGCTGGTC > W3110S.gb/911762‑911868 | cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:1646864/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:1670135/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:2050231/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:2482493/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:2487909/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:3716829/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:4354524/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:447518/1‑62 (MQ=255) cGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGAGCGTaa > 1:4700562/1‑62 (MQ=255) tCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGttt > 1:3355419/1‑62 (MQ=255) gcgcCTCCAGTAGCGGGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGt < 1:674047/56‑1 (MQ=255) gTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAg < 1:2364926/61‑1 (MQ=255) gTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAg < 1:4191379/61‑1 (MQ=255) gTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAg < 1:505386/61‑1 (MQ=255) tAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAg < 1:950690/60‑1 (MQ=255) tgCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAGACCGCTGGTc > 1:1029456/1‑62 (MQ=255) | CGACAACCGGAACGTTATTGCTTTCCAGCGCCTCCAGTAGCGTGGTGCCAACCGGGGCGTAAAATTCTCGTGCCGGTTGCAGTTTGGTGAATTTCAGACCGCTGGTC > W3110S.gb/911762‑911868 |
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: ‑ |