Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,464,506 | 0 | A | G | 53.7% | 32.4 / 48.5 | 41 | E189G (GAA→GGA) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (6/13); new base G (22/0); total (28/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACA > NZ_CP009273/4464422‑4464590 | gTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCgg > 1:518175/1‑89 (MQ=255) gTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGa > 1:695144/1‑90 (MQ=255) gTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGa > 1:369595/1‑90 (MQ=255) aCCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATa < 2:83560/90‑1 (MQ=255) gATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAc > 1:405325/1‑90 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAACAAt > 2:309801/1‑90 (MQ=255) cTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATAttt > 2:616870/1‑89 (MQ=255) tttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTgagg < 1:513312/90‑1 (MQ=255) ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGAAAAAAAAAATTCTTgaggg > 2:48893/1‑89 (MQ=255) ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGATGCAATGGCGGGGGGAGGGGTAAAAAATATTTTTtaggg > 2:144243/1‑89 (MQ=255) tGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTg < 2:369595/90‑1 (MQ=255) cGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCt < 1:280680/90‑1 (MQ=255) agagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTg < 1:256724/90‑1 (MQ=255) gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGGCTTGAccc > 2:280680/1‑90 (MQ=255) gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTTTGGAGGGATGCGGGGCTTTAccc > 2:716151/1‑90 (MQ=255) aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGaaa > 1:426679/1‑90 (MQ=255) ctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGCCTTTACCCCAAAAGc > 1:683127/1‑90 (MQ=255) ctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGGGGAGGGATGCGGGCCGTGACCCGAAAGGc > 2:529677/1‑90 (MQ=255) acaacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTTAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAg > 1:328913/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAacgc > 2:295317/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc > 1:324324/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGATGCGGGCCTTTACCCGAAAAGGGGGGAAacgc > 2:261153/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGgaaaggc > 2:215025/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 2:716108/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGAAAgcgc > 1:529756/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAACGGTTAAAACTATTTTTTTGGAGAGATGCGGGGCGTTACCCCAAAGGAGGGaaaagag > 2:306640/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc > 1:600147/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 1:134023/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 1:33199/1‑90 (MQ=255) tctGGGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCt < 2:683127/90‑1 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc > 2:200670/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc > 1:275030/1‑90 (MQ=255) cAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGAt < 1:113365/90‑1 (MQ=255) aaTGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAGCGCCCGCCGATAAATtgtggt > 2:507015/1‑88 (MQ=255) aaTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGAt < 1:361126/90‑1 (MQ=255) aTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGGGGGCCTTGACCCGAAAGGCGGGAACGCGCCTGCCGATAAGTTGTGATa < 2:529107/90‑1 (MQ=255) aTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATa < 1:144243/90‑1 (MQ=255) aTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCAAAAGGGGGGAAAGCGCCTGCCGATAAGTTGTGATa < 2:405325/90‑1 (MQ=255) ggCGGGGGAAGCGATAAAACCTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg < 1:254304/90‑1 (MQ=255) ggCGGGGGAAGCGATAAAAACTATTCTTGAGGATTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg < 1:200613/90‑1 (MQ=255) gggggAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACa > 1:234514/1‑90 (MQ=255) aaGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa > 2:682010/1‑70 (MQ=255) aaGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 1:682010/70‑1 (MQ=255) | GTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACA > NZ_CP009273/4464422‑4464590 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |