Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 262,073 | 0 | A | G | 52.9% | 14.5 / 22.4 | 17 | G64G (GGT→GGC) | BW25113_RS25425 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (1/7); new base G (9/0); total (10/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGA > NZ_CP009273/261992‑262155 | cTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAAGGGGGGCCGGATaa > 1:249099/1‑90 (MQ=25) tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCa < 1:306549/90‑1 (MQ=255) aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAGGGAAAACCCCCAACGTGGGGCCGGGTAAGGAGTCAAcc > 2:48879/1‑90 (MQ=25) tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGCGTAAAAACCGCAAAGGGGGGGCGGGTAAATAGTCAACCcgtggc > 2:255166/1‑88 (MQ=16) agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTGCAAAGGGGGGCCGGGTAAAAAGTAAACCTGTGAcc > 2:521435/1‑90 (MQ=18) aaCTTCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGGCCAttttt > 2:169009/1‑90 (MQ=255) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAAGGGGGGCCGGGTAAAAAGTCAACCCGTGACCcttctt > 2:141432/1‑90 (MQ=255) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTGCAAAGGGGGGCCGGGGAAGAAGTCAAACCGTGACCAttctt > 1:237028/1‑90 (MQ=255) tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGc < 1:102866/90‑1 (MQ=255) gCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCt < 2:249099/90‑1 (MQ=255) cGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAAGGGGGGCCGGAGAAGGAGTGAACCCGGGGCCCTTCTTTTCTGGcatcat > 1:20565/1‑90 (MQ=255) cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggt < 1:48879/90‑1 (MQ=255) aTCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGGTAAAAAGTCAAACTGTGACCATTCTTTGCTGACATCATCCtgggg > 2:567430/1‑88 (MQ=16) gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATc < 2:41025/90‑1 (MQ=255) aaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCgg < 1:449046/90‑1 (MQ=255) gCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgggg > 1:323178/1‑90 (MQ=255) aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa < 1:567430/90‑1 (MQ=255) | CTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGA > NZ_CP009273/261992‑262155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |