Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,090,826 | 0 | T | G | 57.1% | 8.3 / 10.7 | 14 | G67G (GGT→GGG) | hemW | radical SAM family heme chaperone HemW |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (8/0); total (10/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.50e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT > NZ_CP009273/3090746‑3090913 | cAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCAGc > 2:440752/1‑90 (MQ=255) gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCt > 2:154753/1‑90 (MQ=255) cATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTg < 2:219660/90‑1 (MQ=255) tGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTTGTACGCCGAGCCCGCctt > 1:376548/1‑90 (MQ=255) tGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCTGCttt > 1:86073/1‑90 (MQ=255) ttAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTccc > 2:451162/1‑89 (MQ=255) tCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTccggcccggc < 2:265258/90‑1 (MQ=255) gATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCGGCCTGCTTTTCGGCCCCGGGGTGCaaaag > 1:372719/1‑88 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCCGCTTTTCGGCCCGGCGATGCAAACgctg > 1:237619/1‑90 (MQ=255) cTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGgcg < 2:35158/90‑1 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:110690/67‑1 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc > 2:110690/1‑67 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCCGCTTTCCGGCCCCGCGATGCAAAAGCTGCTGGGCGGCGGGGGTgcgc > 1:100409/1‑90 (MQ=255) gTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCCGCTTTTCCGCCCGGCGGTGCAAACGCCGCTGGGCGGCGGGCGGGGGCGGTTg > 1:393665/1‑90 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 1:154753/90‑1 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 2:284208/90‑1 (MQ=255) | CAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT > NZ_CP009273/3090746‑3090913 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |