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.3% | 24.3 / 34.1 | 30 | G67G (GGT→GGG) | hemW | radical SAM family heme chaperone HemW |
Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (16/0); total (20/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTA > NZ_CP009273/3090745‑3090914 | tCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAg > 2:992556/1‑90 (MQ=255) tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCCAGCCCGCt > 2:108667/1‑90 (MQ=255) tGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGATGTAAAGACAATCTTTATTGGCGGGGGGACGCCCAGGCCGCttt > 2:301216/1‑90 (MQ=255) aaCGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTcccgc > 2:8426/1‑90 (MQ=255) aTCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCGAGCCCGCTTTtcggcccgg > 2:325527/1‑90 (MQ=255) cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGAGCCTGCCTTCCGGCCCcgcg > 2:133311/1‑90 (MQ=255) cAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCGCGCCCGCTTTTCGGCCCCGCGCTgga > 1:492090/1‑88 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGGGCCCGCTTTCCGGCCCGGCGGAGCaaa > 1:1465284/1‑90 (MQ=255) aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGTAAGCCGAGCCTGCTTTCCCGCCCCGCGGTGCaaa > 2:91972/1‑90 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:546588/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 2:1552119/90‑1 (MQ=255) cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc < 1:95840/90‑1 (MQ=255) gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCCGCCCCGCGATGGAAACgctg > 2:445387/1‑90 (MQ=255) gCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCCGCCCGGCGGTGCAAACGctgctg > 2:1130536/1‑90 (MQ=255) aCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGACGCCGGGCCCGCTTTCCCGGCCCGCGGTGCAAACGCTGCTggggg > 1:482531/1‑87 (MQ=255) gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTTCTTTTCCGCCCGGCGGTTCAAAAACTGCTGGGCGGCGGGCGTgggc > 2:1392542/1‑90 (MQ=255) aaGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGtt < 2:492090/90‑1 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCGGCCCGGCGGTGCAAAAGCGGCTGGGCGGCGGGCGGGCGCGtttt > 1:538543/1‑89 (MQ=255) gTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCCGCTTTTCGGCCCCGCGATGCAAAAGCTGCTGGGCGGCGGGCGGGCGCGGTTg > 1:1292223/1‑90 (MQ=255) aaTCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGc > 2:1395551/1‑90 (MQ=255) aaTCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCCGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGGTTGCCGCTggg > 2:606093/1‑89 (MQ=255) aTCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGGGCGCGGTTGCCGCTGGCa > 2:812603/1‑90 (MQ=255) tCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAg < 1:992556/90‑1 (MQ=255) tCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCGGCCCCGCGGTGCAAACGCTGCTGGGCGGCGGGCGGGCGCGGTTGCCGCCGGcgg > 1:286130/1‑88 (MQ=255) tATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGa < 1:1347948/90‑1 (MQ=255) tATTGGCGGTGGGACGCCGAGCCCGCTTTCCCGCCCGGCGGGGCAAAAGCTGCTGGGCGGCGGGGGTGCGCGTTTGCCGCTGGGCGCggg > 1:68675/1‑89 (MQ=255) ttGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATg < 2:538543/90‑1 (MQ=255) tGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGc < 1:785141/90‑1 (MQ=255) gCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAg < 1:301216/90‑1 (MQ=255) cGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGa < 1:721880/90‑1 (MQ=255) ggtggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGaa < 1:75831/90‑1 (MQ=255) gtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTa < 1:1392542/90‑1 (MQ=255) | TCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTA > NZ_CP009273/3090745‑3090914 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |