Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,546,603 | 0 | C | T | 26.7% | 43.8 / 17.2 | 30 | L15L (CTC→CTT) | gntT | high‑affinity gluconate transporter |
Reads supporting (aligned to +/- strand): ref base C (15/7); new base T (3/5); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACGTTACCCATAACAAATGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGC > NC_000913/3546523‑3546682 | gACGTTACCCATAACAAATGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTc > 2:681679/1‑81 (MQ=255) gACGTTACCCATAACAAATGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTc > 2:681667/1‑81 (MQ=255) aCCCATAACAAATGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCtgat > 2:681666/1‑80 (MQ=255) aaTGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTc > 1:681677/1‑80 (MQ=255) aTGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCa > 2:681686/1‑80 (MQ=255) tGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCaa > 2:681662/1‑80 (MQ=255) gCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAAc < 1:681628/81‑1 (MQ=255) ggTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCt > 2:681688/1‑81 (MQ=255) gTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCtt < 2:681629/81‑1 (MQ=255) tAAATCATGCCATTAGTCATTGTTGCCATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTc < 2:681630/81‑1 (MQ=255) aaaTCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCa > 1:681673/1‑81 (MQ=255) aaTCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCTTCAt < 1:681631/81‑1 (MQ=255) cATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGc > 2:681685/1‑81 (MQ=255) aTGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGct > 2:681676/1‑81 (MQ=255) ccATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCtctc > 1:681675/1‑81 (MQ=255) gTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCt < 2:681633/80‑1 (MQ=255) gTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTc < 1:681632/81‑1 (MQ=255) tgttgCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGgc > 1:681670/1‑81 (MQ=255) ttgCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGgcgc > 1:681671/1‑81 (MQ=255) cTATCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCtt < 2:681634/80‑1 (MQ=255) tCGGTGTAATCTTGTTGTTGCTTCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCt < 2:681636/80‑1 (MQ=255) tCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCt < 2:681635/80‑1 (MQ=255) tCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTg > 1:681682/1‑81 (MQ=255) ggTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGTCGCTTGCTGtt < 1:681637/81‑1 (MQ=255) gtgtAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTg > 1:681711/1‑81 (MQ=255) gtgtAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTg > 2:681691/1‑81 (MQ=255) cTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAAt > 2:681687/1‑81 (MQ=255) gttgttgCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAAt < 2:681638/78‑1 (MQ=255) gttgCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAg > 1:681674/1‑80 (MQ=255) gCTTCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAAt > 2:681680/1‑81 (MQ=255) cTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAAt < 1:681639/80‑1 (MQ=255) tCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGc > 1:681681/1‑81 (MQ=255) | GACGTTACCCATAACAAATGAAAGGCCAGGTAAATCATGCCATTAGTCATTGTTGCTATCGGTGTAATCTTGTTGTTGCTCCTGATGATCCGCTTCAAAATGAACGGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGC > NC_000913/3546523‑3546682 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |