Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,464,532 | 0 | T | G | 61.1% | 9.4 / 12.8 | 18 | *198G (TGA→GGA) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (11/0); total (13/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.45e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACG > NZ_CP009273/4464456‑4464600 | agagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTg < 1:145098/90‑1 (MQ=255) gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGGGGGATGCGGGGCTTGAccc > 2:145098/1‑90 (MQ=255) ctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGc > 1:295139/1‑90 (MQ=255) tataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCg < 1:302199/90‑1 (MQ=255) caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGAGGGATGCGGGCCTTGACCCGGAAGGCGGGgaaac > 2:302199/1‑88 (MQ=255) aacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGGGGGCCCTGACCCGAAAGGCGGGGAAacg > 1:361027/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAacgc > 2:129428/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc > 2:33240/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTTGGAGGGATGCGGGGCTTGACCCGAAAGGCGGGAAAgcgc > 2:403787/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGAGGGGTGCGGGCCCTGACCCGAAAAGCGGGgaaaagc > 2:449546/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc > 2:36330/1‑90 (MQ=255) ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc < 1:403787/90‑1 (MQ=255) tctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGG‑CTTGACCCGAAAGGCGGGGAAGCGCCTg > 2:41968/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc > 1:54784/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc > 1:58885/1‑90 (MQ=255) gCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGATTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGa < 2:295139/90‑1 (MQ=255) aCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAACGCCCGCCGATaa > 2:150395/1‑90 (MQ=255) aTAAAACCTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAcg < 1:385231/90‑1 (MQ=255) | AGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACG > NZ_CP009273/4464456‑4464600 |
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: ‑ |