Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,089,567 | 0 | T | G | 56.2% | 9.7 / 17.3 | 16 | G81G (GGT→GGG) | hisF | imidazole glycerol phosphate synthase subunit HisF |
Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (9/0); total (10/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCA > NZ_CP009273/2089480‑2089654 | cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTgg > 2:458673/1‑90 (MQ=255) cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTa < 2:156739/90‑1 (MQ=255) tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg < 1:363662/90‑1 (MQ=255) tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg < 2:556363/90‑1 (MQ=255) tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTggg > 1:225515/1‑89 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGT‑TCTGGGAGATGc > 1:462335/1‑90 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGGAGATg > 1:839/1‑90 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAATCTCTCGaaaaag > 1:328689/1‑86 (MQ=255) gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGcc < 2:186184/90‑1 (MQ=255) gATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGATGCCGCGaaa > 2:53060/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTCTCTGGGAGAAGCCGCGAAAATTTTTTTCt > 2:130161/1‑90 (MQ=255) ctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcg < 1:516511/90‑1 (MQ=255) ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCTTGGAGAAGCCCCGAAAAATTTTTTCTTTTGGGCGGGTaaaaa > 2:516511/1‑89 (MQ=255) ttGTGTGGCGGGGGGGGTTAAATCTCTGGAAGATGCCCCGAAAATTTTTTTCTTTTGGGCGGGTAAAAATTTCATCAAATCTCCTgcgcc > 2:424201/1‑89 (MQ=255) ttGTGTGGCGGGGGGGATTAAATCTCTGGAAGATGCCCCGAAAAATTTTTCCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTGCGCt > 1:506451/1‑90 (MQ=255) ggTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCa < 1:53060/90‑1 (MQ=255) | CACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCA > NZ_CP009273/2089480‑2089654 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |