Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 4,807,656 | 0 | G | C | 27.1% | 99.6 / 10.2 | 49 | S66T (AGC→ACC) | mcbC | Microcin B17‑processing protein McbC |
Reads supporting (aligned to +/- strand): ref base G (27/8); new base C (0/13); total (27/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCTTTCCATGCCAAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACT > KT2440/4807522‑4807789 | ggCTTTCCATGCCAAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAg < 8:58429/135‑1 (MQ=255) cTTTCCATGCCAAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCGCCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGc > 1:102789/1‑134 (MQ=255) atttCCATGCCAAAGTCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGc < 4:36477/133‑1 (MQ=255) aTGCCAAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAAc > 1:126009/1‑134 (MQ=255) aaaGGAAACTATACGAACCACAAACCAACTCAACACTTGAGCGCCCTCCACAAAATCCCACCAAACCACCTCTAAAAACTTACCGGAAACGAACAAAAACTCAATACA‑ATTGTAGCCCCAACAAAAAACaaacacc < 8:281615/136‑5 (MQ=255) agtcAAATATACTATAAACAAAGCAACTCAACACTTGAGCGCCCTGCACAAAATACCACCAAAAGACCTCTAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCACAAGAAAAAACAAAGAAcc < 1:210629/132‑1 (MQ=255) ataCGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCtt > 5:49845/1‑134 (MQ=255) tACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACctc > 2:197259/1‑135 (MQ=255) acaAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTATgg > 5:191240/1‑135 (MQ=255) aaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGa > 1:185251/1‑135 (MQ=255) aaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGa > 3:290241/1‑135 (MQ=255) aaCTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGccc > 2:18580/1‑135 (MQ=255) tCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCacac > 7:59159/1‑136 (MQ=255) accacTTCCGCCCCCTCCACAAAAAACAACAAAAACCCCGCCAAAAACTTCCCGGAAACCAACAAAAACTCAAAACC‑TTTGTACCCCCAAGCAAAAACAAAGACCCGACTCTTACACCTCTGGCGAGCCCACAcc < 3:76957/133‑1 (MQ=255) aaCACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACAcc > 2:203033/1‑135 (MQ=255) acacTTCAGCGCCCTCCACAAAATCCCACCAAAAGACCTCGAAAAACTTCCCGCAAACCAACTAAAACTCAAAACA‑ATTGAACCCCCAACCAAAAACAAACAACCGACTCTTACACCTCTGGCGAGCCCACACCCg < 7:23623/136‑1 (MQ=255) cTTCAGCGCCCTGCAACAAAAACCACAAAACCACCTCAAAAAACTTCCCGGAACCCAAATAAAACTCAAAACC‑ATTGAAGCCCCAACAAAAAACAACCAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCt < 3:52945/135‑1 (MQ=255) ttGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCCCACCCGCTa > 1:213487/1‑135 (MQ=255) ttGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTa > 1:39334/1‑135 (MQ=255) cgcCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGCCTCTTACACATCTGGCGAGCCCCCACCCGCTACgaa > 6:40039/1‑134 (MQ=255) tGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATa > 1:288068/1‑135 (MQ=255) aaacaaccaAAACACCTCCAAAAACTTACCGCAAACCAACAAAAACACAAAACC‑ATTCTACCCCCAACCAAAAACAAACACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGc < 1:197259/130‑1 (MQ=255) accaccAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGa > 2:117643/1‑135 (MQ=255) caAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGAt > 4:139811/1‑135 (MQ=255) gACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGaa > 5:299458/1‑135 (MQ=255) gACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGaa > 5:159747/1‑135 (MQ=255) cGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACg > 6:56432/1‑135 (MQ=255) aaaaCTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGccc > 5:146369/1‑135 (MQ=255) aaaaTTACCGTAACCGAACAAAACCTCAAAACA‑ATTTTACCCCCAACCAAAAAAAAACACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct < 5:78490/135‑1 (MQ=255) aaaCTTCCCGGACACGAACTAAAACTCAAAACC‑ATTGTAGCCCCAAACAAAAACAAACACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct < 2:185251/135‑1 (MQ=255) aaaCTTCCCGGAACCAAACTAACACTAAATACC‑AATTACCCCCCAACCAAAAACACCCACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct < 2:101810/135‑1 (MQ=255) tacccGCAAACCAACAAAAACTCAAAACC‑ATTGTACCCCCAAGCAAAAACAAACAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTc < 3:197330/132‑1 (MQ=255) acccGCAAACTAACTAAAAC‑CAAAACCAATTGAACCCCCAACCAAAAACAACCCCCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCa < 2:288068/133‑1 (MQ=255) gcaaaCTAACTAAAACTCTAAACC‑TTTTTAGCACCAGCAAAAAACAATCAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg < 8:125952/134‑1 (MQ=255) gAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGAAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg < 4:290241/135‑1 (MQ=255) gAAACCAACTAAAACTCAATACC‑ATTGAAGCCCCAACCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg < 6:49845/135‑1 (MQ=255) gAAACCAACAAAAACACAATAAC‑ATTGTACCCCCAACCAAAAACAAATAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg < 6:146369/135‑1 (MQ=255) aaaaCCAACAAAAAATAAAACCA‑ATTGTACCCCCAACAAAACAAACACACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg < 2:126009/134‑1 (MQ=255) ccaaCTAAAACACTATACA‑ATTGACCCCCCCCGAAAAAAAAATCACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGatat < 3:254875/133‑1 (MQ=255) tAAAACTCAATACC‑ATTGTACCCCCAAGAAAAAACAAACAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCa < 8:23232/136‑1 (MQ=255) aaaaCTCAATACC‑TTTGTAGCCCCACCCAAAACCAAACACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCa < 2:53160/135‑1 (MQ=255) tCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCAt < 6:299458/135‑1 (MQ=255) cac‑ATTGTACCCCCCACAAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCAttgctt < 3:76107/133‑1 (MQ=255) c‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCAttgcttgc > 4:78170/1‑135 (MQ=255) tAGCCCCAAGCAAAAACAAAGCACCGACTCTTACACCTCTGGCGAGCCCCCACCCGCTACGAAGAAATAGTTCATGAGAACGATTTGAAGAACGCCCTCTTCACTTTGATATAGTCCAAGCCTTGCTTGCACCgat > 8:1812/1‑134 (MQ=255) tAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGct > 8:252586/1‑136 (MQ=255) tAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGc > 3:304806/1‑135 (MQ=255) ccAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTaa > 7:74561/1‑136 (MQ=255) aGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAAc > 1:319019/1‑134 (MQ=255) aGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACt > 7:165155/1‑135 (MQ=255) | GGCTTTCCATGCCAAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACC‑ATTGTAGCCCCAAGCAAAAACAAAGAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACT > KT2440/4807522‑4807789 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |