Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,732,535 | 0 | A | C | 55.6% | 19.7 / 22.4 | 19 | Y560S (TAC→TCC) | malS | alpha‑amylase |
Reads supporting (aligned to +/- strand): ref base A (5/3); new base C (0/10); total (5/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.54e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGC > NZ_CP009273/3732446‑3732611 | ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTa > 1:230178/1‑90 (MQ=255) aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCtcgt < 2:357293/90‑1 (MQ=255) gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCAt < 2:220507/90‑1 (MQ=255) gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCAt < 1:101551/90‑1 (MQ=255) cTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGGGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTCCCTCTCGTCGCATGAtcc < 1:251755/90‑3 (MQ=255) tGGCGCAGATGGATACGCCCTGGCAGAAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTCCCTCTCGTCGCATGATACCCGcc < 1:200709/90‑1 (MQ=255) tGGCGCAGATGGATACGACCGGGCAGCAAATGGCGGAGAATTTGCAGGGTTTCACCGTGTTGAGCTCCCTCTCGTCGCATGATACCCGcc < 1:388532/90‑1 (MQ=255) gATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCg > 1:413897/1‑90 (MQ=255) aTACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAg > 1:84857/1‑90 (MQ=255) cGACCTGGCAGCAAATGGCGGAAAAATTGCAGGTTTTCAACGTGTTGAGCTCCCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgggg < 1:41257/90‑1 (MQ=255) aCCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGc > 1:442138/1‑90 (MQ=255) aCCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGc > 1:379085/1‑90 (MQ=255) tGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGAc > 2:440445/1‑90 (MQ=255) gcagaaAATGGCGGAAAAATTGCGGGTTTTCAACGTGTTGAGCTCCCTCTCGTCGCATGATCCCCGCCTGTTCCGTGAAGGGGGCGACaa < 2:379085/90‑1 (MQ=255) caacaACTGGGGGAGAATTTCCAGTGTTTCAACGTTTTAAGCTGCCCCTCGTCGAATGTAACCCGCCTGTTCCGTGAAGGGGGCGACaaa < 1:265069/87‑1 (MQ=255) aTGGCGGACAAATTCCGGGGTTTAAAGGTGTGGGCCTCCCTCCCCTCGCTTGATCCCCGCCTGTTCCGTGAAGGGGGCGACAAAgcagca < 1:299246/90‑1 (MQ=255) cGGAAAAATTGCAGGGTTTCAACGTGTTGAGCTCCCTCTCCTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGt < 1:422146/90‑1 (MQ=255) cAGGTTTTCAACGTGTTGAGCTCCCTCTCGTCGCATGATCCCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTa < 2:430839/90‑1 (MQ=255) cAGGTTTTAAACGTGTTGAGCTCCCTCTCGTCGCATGATCCCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTa < 2:84857/90‑1 (MQ=255) aaCGTGTGGGGCTCCCCCTCGTCGCATGATCCCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGgc < 2:414698/90‑1 (MQ=255) | GGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGC > NZ_CP009273/3732446‑3732611 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |