| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 262,051 | 0 | T | G | 52.6% | 17.3 / 22.7 | 19 | T72P (ACC→CCC) | BW25113_RS25425 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (10/0); total (11/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATC > NZ_CP009273/261962‑262135 | cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGt < 1:187925/90‑1 (MQ=255) ttCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGt > 2:211560/1‑90 (MQ=255) tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTa < 1:211560/90‑1 (MQ=255) tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTaaaa < 1:78515/90‑1 (MQ=255) caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTAGGTaaaa > 1:41025/1‑90 (MQ=25) aGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAAcc < 2:23187/90‑1 (MQ=255) cTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGTAAAAACTGCAAAGGGGGGCCGGATaa > 1:249099/1‑90 (MQ=25) tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCa < 1:306549/90‑1 (MQ=255) aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAGGGAAAACCCCCAACGTGGGGCCGGGTAAGGAGTCAAcc > 2:48879/1‑90 (MQ=25) tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGGCGTAAAAACCGCAAAGGGGGGGCGGGTAAATAGTCAACCcgtggc > 2:255166/1‑88 (MQ=16) agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGGAAAAACTGCAAAGGGGGGCCGGGTAAAAAGTAAACCTGTGAcc > 2:521435/1‑90 (MQ=18) aaCTTCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGTAAAAACCGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGGCCAttttt > 2:169009/1‑90 (MQ=255) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGTAAAACCTGCAAAGGGGGGCCGGGTAAAAAGTCAACCCGTGACCcttctt > 2:141432/1‑90 (MQ=255) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGACGGAAAAACTGCAAAGGGGGGCCGGGGAAGAAGTCAAACCGTGACCAttctt > 1:237028/1‑90 (MQ=255) tGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGc < 1:102866/90‑1 (MQ=255) gCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCt < 2:249099/90‑1 (MQ=255) cGCAGCCATCGGGCGACTGGC‑GGGGACGTAAAAACTGCAAAGGGGGGCCGGAGAAGGAGTGAACCCGGGGCCCTTCTTTTCTGGcatcat > 1:20565/1‑90 (MQ=255) cATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggt < 1:48879/90‑1 (MQ=255) aTCGGGCGACTGGC‑GGGGACGTAAAAACCGCAAAGGGGGGCCGGGTAAAAAGTCAAACTGTGACCATTCTTTGCTGACATCATCCtgggg > 2:567430/1‑88 (MQ=16) gC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATc < 2:41025/90‑1 (MQ=255) | CGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATC > NZ_CP009273/261962‑262135 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |