Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,090,826 | 0 | T | G | 55.8% | 47.2 / 52.7 | 43 | G67G (GGT→GGG) | hemW | radical SAM family heme chaperone HemW |
Reads supporting (aligned to +/- strand): ref base T (7/12); new base G (24/0); total (31/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.29e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT > NZ_CP009273/3090738‑3090913 | aTTTTGTTCAGCATCTGCTTAACGATGTGGACAACGATGTGGCTTACGCTCTGGGCCGTGAAGTAAAGACAATCTTTATTCGCGGGGGTa > 1:763425/1‑90 (MQ=255) gTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCg > 2:159105/1‑90 (MQ=255) tCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAg > 2:1424193/1‑90 (MQ=255) cAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGc > 2:1596246/1‑90 (MQ=255) gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCt > 2:1634119/1‑90 (MQ=255) gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGccg > 1:1297468/1‑89 (MQ=255) aTCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGc < 1:1534676/90‑1 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGCCTGcc > 1:1417040/1‑89 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCCGCt > 2:401645/1‑90 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCTGCt > 2:1382032/1‑90 (MQ=255) tGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCAGCCCGCttt > 2:1492709/1‑90 (MQ=255) aTCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTccggcccgg < 1:47581/90‑1 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCCGCTTTCCCGCCCGGCg > 2:120835/1‑90 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCCCGCCTGCTTTTCCGCCCcgcg > 2:1028011/1‑90 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGAACCCGCGCCCGGCTTccgggcccgtg > 1:1527093/1‑86 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCCGCTTTTCCGCCCcgcg > 1:58278/1‑90 (MQ=255) gACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATg < 2:482728/90‑1 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCCAGCCTGCTTTTCGGCCCGGCGGTGCaaa > 2:352667/1‑90 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGGCGCCGGGCCCGCCTTTCCGCCCGGCGGTGGaaa > 1:1310856/1‑90 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCGGCCTGCTTTTCGGCCCCGCGATGCaaa > 2:1472989/1‑90 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:1063466/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:657199/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:1028011/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:1424193/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 2:516183/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:1634119/90‑1 (MQ=255) aTGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACgc < 2:279969/90‑1 (MQ=255) aTGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGCCCGCTTTTCCGCCCCGCGGTGCAAAAgc > 2:1184544/1‑90 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCCGCGATGCAAACgctg > 2:342400/1‑90 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGTACGCCGAGCCCGCTTTTCGGCCCCGCGATGCAAACgccg > 1:203250/1‑88 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACgctg > 1:648011/1‑90 (MQ=255) ttACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAAGCTGCTggg > 2:225109/1‑89 (MQ=255) ttACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCCGCGGTGCAAAAGCTGCTggg > 1:554284/1‑89 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGAc < 1:1184544/90‑1 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGAc < 1:509761/90‑1 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCGGCCCCGCGGTGCAAAAGCTGCTgggc > 2:1374624/1‑88 (MQ=255) tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACCCCGGGGCTGCTTTCCGGCCCGGCGGTGCAAACGCTGCTgggg > 2:1353103/1‑88 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCGGCCCGGCGGTGCAAACACCGCTGGGCGGGGGGGGTgcgc > 1:1340492/1‑90 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCCGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTgcgc > 2:609969/1‑90 (MQ=255) aTCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCCGCCCCGCGATGCAAACGCTGCTGGGCGGGGGGGGGGCGCGGTTGGCGCCGGcg > 1:1077140/1‑89 (MQ=255) ggCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGGGGGCGTTTGCCGCTGGCAGCGGGTGCa > 2:607152/1‑90 (MQ=255) cGGTGGTACGCCGCGACTGCTTTTCGGCCCCGCGGTGCAAAAGCGGCTGGGCGGGGGGGCCGGGCGCTTTTCGCCGGCGGCGGATGCAGa > 1:40230/1‑90 (MQ=255) cGGTGGGACGCCGAGCCTGCTTTTCGGCCCGGCGGTGGAAAAGCTGGTGGGGGGCGGGGGGGGGGGGTTGCCGCTGGGCGCGGGTGCAGa > 1:783511/1‑90 (MQ=255) ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt < 1:493076/90‑1 (MQ=255) | ATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT > NZ_CP009273/3090738‑3090913 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |