Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,186,017 | 0 | G | T | 55.6% | 17.7 / 29.7 | 18 | T4K (ACA→AAA) | thiC | phosphomethylpyrimidine synthase ThiC |
Reads supporting (aligned to +/- strand): ref base G (1/7); new base T (10/0); total (11/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGA > NZ_CP009273/4185933‑4186106 | cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGGAGGTtttttgta > 2:663117/1‑88 (MQ=255) gTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGCAGGttttttttag > 1:673362/1‑86 (MQ=255) ggCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaa < 1:626283/90‑1 (MQ=255) cGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaaaa < 1:155574/90‑1 (MQ=255) cGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaaaa < 2:116990/90‑1 (MQ=255) gTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTGCAGACATAGCTCA‑TTCCaaaaaa > 1:105162/1‑89 (MQ=255) gTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTTGAGAAAAAGCTCATTTCCaaaaa > 1:656980/1‑90 (MQ=255) cTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTaa < 1:496955/90‑1 (MQ=255) gtgtCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCAGACATAACTCC‑TTTCAAAAAAATAAAGACgggg > 2:160560/1‑90 (MQ=255) gtCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCCGGTTTGTTGGAGAAATAGCTCA‑TTTCAAAAAGTTAAAGACgggggt > 1:440519/1‑88 (MQ=255) tCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTTTTGCAGACATAGCTCA‑TTCCCAAAAAATAAGGACGTGGGtt > 1:333810/1‑90 (MQ=255) cGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGGCGGGGGTTg > 1:119885/1‑90 (MQ=255) gggCGCGTTGTTCGCGGCGGGTCAGTTTTTTTTCAGACAAAGCTCA‑TTTCAAAAAATTAAGGGCGGGGGTTTTTAGACCACGGAGGaaaa > 2:73266/1‑88 (MQ=255) gggCGCGTTGTTCGCGGCGGGGCAGGTTTTTTGGGGACAAAGGTCC‑TTTCCAAAAGGTAAAGGCGGGGGTTGTTTGGCGAGGGGTGAAgt > 2:593102/1‑89 (MQ=255) ggCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCAGACATAGCTCA‑TTCCAAAAAAATAAGGACGGGGCTTTTTAGAAGACGGGTGAAgga > 2:127962/1‑88 (MQ=255) tCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATc < 2:656980/90‑1 (MQ=255) gggTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAg < 1:342850/90‑1 (MQ=255) gggcAGTTTTTTTTCAGACATAGCTCA‑TTCCAAAAAATTAAGGACGGGGGTTGTTAGACGACGGATGAAGCAAGAGAAGAACGCCGCAgg > 2:260313/4‑90 (MQ=255) gggcAGTTTTTTTGCAGACATAGCTCA‑TTTCAAAAAAATAAAGAGGGGGCTTTTTGGACGACGGGGGAAAGAAGAGAAGAACGCCGCAgg > 2:628439/4‑90 (MQ=255) cAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc < 1:274019/90‑1 (MQ=255) gttgCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcga > 2:756686/1‑90 (MQ=255) | CGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGA > NZ_CP009273/4185933‑4186106 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |