Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,090,826 | 0 | T | G | 68.6% | 4.1 / 18.8 | 35 | G67G (GGT→GGG) | hemW | radical SAM family heme chaperone HemW |
Reads supporting (aligned to +/- strand): ref base T (2/9); new base G (24/0); total (26/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.79e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-07 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTA > NZ_CP009273/3090737‑3090914 | gaTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCggtggt < 1:1006221/90‑1 (MQ=255) aTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCgggggtg > 1:240075/1‑89 (MQ=255) aGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGcc < 2:1306945/90‑1 (MQ=255) cATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTg < 1:1755532/90‑1 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGAAGCCGGGCCTGCt > 2:1773419/1‑90 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACCCCGAGCCTGCt > 2:1491204/1‑90 (MQ=255) ttAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCttttcc > 2:629582/1‑87 (MQ=255) ttAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGAACCCGAGCCCGCTTTccg > 2:763603/1‑90 (MQ=255) tAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGGGCCCGCGTTtcgg > 1:441291/1‑90 (MQ=255) tCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCGGCCTGCCTTccggcccggc > 2:761966/1‑90 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCCAGCCCGCGTTCCGGGCCcgcg > 1:1573101/1‑90 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGAGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCaaa > 2:503527/1‑90 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:761966/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTACGGCCCGGCGATGCAAAc < 2:992288/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACCCCCAGCCCTCCTTCCTGGCCGGGGAGGGAAAc > 2:63112/1‑90 (MQ=255) aTGTGGCTTCCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACgc < 2:522915/90‑1 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGAACCCCAGCCCGCCTTTCGGCCCGGGGATGGAAACgctc > 1:856096/1‑89 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTCCTGGAc < 2:1115473/90‑1 (MQ=255) caggCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACTCTGCTGGAc < 2:626222/87‑1 (MQ=255) ccGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCCGCCCGGCGGTGCAAAAGCTGCTGGACGGCgggcgggc > 2:1713962/1‑90 (MQ=255) ccGTGAAGTAAAGACAATCTTTATTGGCGGGGGGCGCCCCGGCCCGCTTTTCCGCCCGGCGCCGCCAACGCCGTTGGTGGgggtgcgggt > 2:1166453/1‑87 (MQ=255) ccGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGGCCGCTTTTCCGCCCCGCGATGCAAAAACTGCTGGGCggcgggcgggg > 1:1706022/1‑87 (MQ=255) cGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGGGCCCGCTTTCCGGCCCCGCGATGGAAAAGCTGCTGGGCggcgggcgggcg > 1:799044/1‑90 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCGGGCCTGCTTTTCGGCCCGGCGGAGCAAAAACTGCTGGACGGCGggcgtgcgg > 1:33509/1‑89 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGCGCCCGCTTTTCCGCCCGGGGGGGCAAAACCTGCTGGGAGGCGGGGGTgggc > 1:1786436/1‑90 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCGGCGGTGCAAACGCTGCTGGGCGGGGGGCGGGGGCGgttt > 2:1357697/1‑89 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGGACGCCCAGCCTGCTTTCCGGCCCGGCGGTGCAAACGCTGCTGGACGGGGTGCGTGCGCGGTTg > 2:245126/1‑90 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGGACGCCCAGCCCGCTTTTCCGCCCCGCGATGCAAACGCTGCTGGACGGCGTGCGGGGGCGtttt > 1:1540348/1‑89 (MQ=255) gTAAAGACAATCTTTATTGGCGGGGGGACGCCCAGCCCGCTTTCCGGCCCGGCGGTGCAAACGCTGCTGGGCGGGGGGCGGGCGCGGTTg > 2:1297654/1‑90 (MQ=255) tAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGATGCTGGACGGCGTGCGTGCGCGTTTGc < 1:1660029/90‑1 (MQ=255) aaTCTTTATTGGCGGTGGGGCGCCGAGCCTTCTTTTCGGCCCCGCGGTGCAAAAGCTGCTGGGCGGGGGGGGGGCGCGGTTGGCGCCggg > 1:1276974/1‑89 (MQ=255) aTCTTTATTGGCGGTGGGACGCCGGGCCTGCTTTCCGGCCCGGCGATGCAAAAACCGCTGGGCGGGGGGCGTGGGGGGTTTGCGCTGGcc > 1:507020/1‑89 (MQ=255) tCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCCGCGATGCAAACACTGCTTTACGGCGGGCGGGCGCGGTTTCCGCTGGGAg > 2:872976/1‑90 (MQ=255) tCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCCGCGATGCAAAAGCTGCTGGGCGGCGGGCGGGCGCGGTTTCCCCTGGCAg > 2:1213319/1‑90 (MQ=255) gCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAg < 1:1627590/90‑1 (MQ=255) cGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGGTGCAGa > 1:1309144/1‑90 (MQ=255) cGGTGGGACGCCGAGCCCGCTTTCCGGCCCGGCGATGCAAAACCTGCTGGGCGGGGGGGGGGCGCGGTTGCCGCCGGGCGCGGGTGCCGa > 1:1330265/1‑90 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:33509/90‑1 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:456108/90‑1 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:1540349/90‑1 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:886442/90‑1 (MQ=255) gtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTa < 2:934711/90‑1 (MQ=255) | GATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTA > NZ_CP009273/3090737‑3090914 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |