Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,090,826 | 0 | T | G | 53.8% | 24.8 / 29.1 | 27 | G67G (GGT→GGG) | hemW | radical SAM family heme chaperone HemW |
Reads supporting (aligned to +/- strand): ref base T (4/8); new base G (14/0); total (19/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.17e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTAC > NZ_CP009273/3090742‑3090915 | tGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGcc < 2:918423/90‑1 (MQ=255) tCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCCAg > 2:371744/1‑90 (MQ=255) gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGCGCCt > 2:367776/1‑90 (MQ=255) cATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTg < 2:315316/90‑1 (MQ=255) aTCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCCG‑CTTTtcggccccg > 2:135943/1‑88 (MQ=255) tCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTccggcccggc < 1:201765/90‑1 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGAAGCCGAGCCCGCCTTTccggcccgcg > 2:289164/1‑88 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCaaa > 1:31913/1‑90 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAAc < 2:1286965/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAAc < 2:1030819/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAAc < 2:1339249/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGCCTG‑CTTTTCGGCCCGGGGATGCAAAc > 1:24172/1‑90 (MQ=255) aTGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACgc < 1:574033/90‑1 (MQ=255) aTGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACgc > 1:320900/1‑90 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGGGCCTG‑CTTTCCGGCCCCGCGGTGCAAACgctg > 2:77078/1‑90 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCTG‑GTTTCCCGCCCGGCGATGCAAAccccg > 1:1103554/1‑86 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGAc < 1:381225/90‑1 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCGGCCCG‑CTTTCCCGGCCGGCGCTGCCAACACTGCCgggc > 2:944531/1‑88 (MQ=255) aCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGCCTG‑CTTGCCGGCCCGGCGGTGCAAACCCCGCTGGACg > 1:1306408/1‑90 (MQ=255) cAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGTACGCCGAGCCTG‑CTTTCCGGCCCCGCGGGGCAAACGCTGCTGGGCGgcgtg > 1:585940/1‑90 (MQ=255) cAGGGCCGCGAAGTAAAGACAATCTTTATTGGCggtggt < 2:1284539/39‑1 (MQ=38) cAGGGCCGCGAAGTAAAGACAATCTTTATTGGCggtggt > 1:1284539/1‑39 (MQ=38) ccGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGGCCG‑CTTTTCCGCCCGGCGGTGCAAACGCTGCTGGGCGGCgggcgggc > 2:179686/1‑90 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTG‑CTTTCCCGCCCGGCGATGCAAACGCTGCTGGACGGCGGGCGTgcgc > 1:987538/1‑90 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGGGGAACGCCGAGACTG‑CTTTCCGGCGCCGCCGGGGAAAAACTGCTGCACGGCGTGGGTgggc > 2:1057793/1‑90 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGGCGGCGTGCGTGCGCGTTTg > 1:200910/1‑90 (MQ=255) aaaGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGcc > 1:37030/1‑90 (MQ=255) tGGCGGTGGGACGCCGAGCCTG‑CTTTTCGGCCCCGCGGTGCAAACGCTGCTGGGCGGCGGGCGGGCGCGTTTGGCGCTGGCAGCgggggc > 2:510607/1‑86 (MQ=255) cGGTGGTACGCCGAGCCCG‑CTTTCCCGCCCGGCGATGCAAACCCTGCTTGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGa > 1:709898/1‑90 (MQ=255) gtggtACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGaaa < 1:1261044/90‑1 (MQ=255) ggtACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:350176/90‑1 (MQ=255) tACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTAc < 2:1346658/90‑1 (MQ=255) | TGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTG‑CTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTAC > NZ_CP009273/3090742‑3090915 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |