Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 1,544,643 | 0 | T | G | 55.2% | 11.3 / 26.1 | 31 | Q412P (CAG→CCG) | pykA | Pyruvate kinase II |
Reads supporting (aligned to +/- strand): ref base T (2/11); new base G (9/7); total (13/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.24e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTG > KT2440/1544515‑1544765 | aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCTGCCGGCAGCGAGGCCGGGTCGAAGGCAATAGGGTACACGCCCCGGCACATCGAGGCGCGCGCCGGGGTgg > 6:253631/1‑135 (MQ=255) aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGCAGGCAATCGGGTACCCGCCACGGAACCTCGAGGCGCCCGCCCgggggg > 7:108083/1‑135 (MQ=255) ccTTGCTCCACCAGGCCAAGTTTGAGCAGTTCGTCGACCCCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAACGAAATCCGGTACACGCCCCGGAACATCCAGGCGCGCGCCAGGGGGgggccg > 8:234208/1‑133 (MQ=255) tCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTgcgg < 2:10203/134‑1 (MQ=255) aGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAgg < 3:273634/135‑1 (MQ=255) gCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGc > 7:228786/1‑135 (MQ=255) aGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGCTCGGTACa > 7:250072/1‑135 (MQ=255) gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGGGGCGCCAGGGCGAAGCTCGGTACAAGCgc > 8:20776/1‑134 (MQ=255) gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGTGCGGCGGCAGGGCGCAGATCGGTTCATGCGa > 3:103151/1‑135 (MQ=255) gtcgACCGCCGCCTGGCTCAACTTATCGGCCGGCCGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCCGAACATCGAGGCGCGCGCCGTGGTTGCGCGGGGGCGCGCCACGGCGAACATCGGGGCCTGCGa > 8:177619/1‑135 (MQ=255) gccTGGCTCACGTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGAACATGCGAACGCAGgcgc < 6:71294/135‑1 (MQ=255) ccTTTTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATGGGCTACACGCCACGGAACAGCGAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCGGGCGCGACATGATCa < 6:12079/135‑1 (MQ=255) gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCCAAGAACCGTTCATGGGAACGCAGGCGCGAAATTATACACGGGgt > 8:229642/1‑134 (MQ=255) gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGTGCGGCGCCACGGGGAAGAACCGTTCATGCGAACGCAGGGGCCGCAAGATAACCGGGgtg > 8:175922/1‑135 (MQ=255) aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCGGAACATCGAGGCGCGCGCCGGGGGGGGGCGGGGCGGCGCCCGGGCGCAGATCGGTAACTGCGCCCGCAGGCGCGACAAGATCAGCGGGGTGTAGCCg > 2:101974/1‑135 (MQ=255) ggCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTc < 2:273429/135‑1 (MQ=255) gCCGGGTCGTAGTCAATGGGGTAGAGGCGACGGAAGATCTAGGCCCGCGCGTGTGTGGCGCGGTGCGGGGACTGGGCGAAGATCGGTTCATGCGAACGCAGGCTCGACATGATCATCGGGGTGTAGCCGCTTTCg < 2:276904/135‑1 (MQ=255) ggcatgcGGGAACACGCCACGCACCAGCGCGGCGCGCGCCGGGGGGGGGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat < 1:253020/129‑1 (MQ=255) ggCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat < 7:26580/135‑1 (MQ=255) ggCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat < 8:122802/135‑1 (MQ=255) ggCAATCGGGAACACGCCAGGGAACTTCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGATGTTCGGTACATGCGAACGCAGGCGCGTCATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGTCgat < 7:177619/135‑1 (MQ=255) gCAAGCGGGAACCCGCCACGGAACATCGCGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgatg < 5:243052/135‑1 (MQ=255) gggTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGGGGGGCGCCAGGGCGAAGATCGGTTCATGCGAACCCAGGCGCGACATGATAAGCGGGGTGTAGCCGCGTTCGGGGACGGCGATGATCGcc > 2:47801/1‑135 (MQ=255) ggTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCt < 8:228786/135‑1 (MQ=255) gCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGcc > 7:194638/1‑135 (MQ=255) ccagcGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt < 1:116470/131‑1 (MQ=255) acctCGCGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt < 5:186683/132‑1 (MQ=255) aCATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt < 7:20776/135‑1 (MQ=255) tCGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt < 2:285880/135‑1 (MQ=255) gcggggcgcgcGGGGGGGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt < 2:297150/128‑1 (MQ=255) aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCTAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg < 8:44788/134‑1 (MQ=255) | AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTG > KT2440/1544515‑1544765 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |