Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,464,518 | 0 | C | A | 32.4% | 49.2 / 10.5 | 35 | T193N (ACT→AAT) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base C (7/16); new base A (11/0); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAC > NZ_CP009273/4464430‑4464599 | aTGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt < 1:225887/90‑1 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAatat > 2:288063/1‑90 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGGTAAAAAatat > 1:297731/1‑90 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt < 2:173156/90‑1 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt < 1:254553/90‑1 (MQ=255) gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt < 1:276478/90‑1 (MQ=255) ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTgag > 1:447787/1‑88 (MQ=255) ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTgag < 2:447787/88‑1 (MQ=255) tCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAAAATTCTTtagggg > 2:530892/1‑88 (MQ=255) gAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATg > 1:87781/1‑90 (MQ=255) agagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTg < 1:530892/90‑1 (MQ=255) gagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGa < 1:240313/90‑1 (MQ=255) aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTTAGGAGGGATGCGGGGCTTGACCCGaaa > 2:460641/1‑90 (MQ=255) aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGaaa < 1:144179/90‑1 (MQ=255) ataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGgggg > 1:346247/1‑90 (MQ=255) gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGATGCGGGGCTTGACCCCAAAAGCGGGaa > 2:302162/1‑90 (MQ=255) caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGCCTTGACCCGAGAGGCGGGgaaag > 2:453502/1‑88 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAATATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGaaaaagc > 2:290784/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc > 1:585043/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGTGATGCGGGCCTTGGCCCGAAAGGCGGGGAAacgc > 2:22784/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAacgc > 2:509429/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGCAGTGACAAAAAGCATTTTTTAGGAGGGCGACAGGCCCTTCCCCGAAATGGGGGAAAacgc > 1:511222/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt < 1:51807/43‑1 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt > 2:51807/1‑43 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 1:76354/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 1:192220/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGGGGAGGGATGCGGGCCTTGACCCGAAAGGGGGGGAAgggc > 1:438003/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAAGCGGGGAAacgc > 2:212137/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGTGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAacgc > 2:375870/1‑90 (MQ=255) ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc < 2:76354/90‑1 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTTGGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAAAGCCCGcc > 2:424220/1‑90 (MQ=255) aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 1:22784/90‑1 (MQ=255) aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 2:87781/90‑1 (MQ=255) aaTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGAt < 1:290784/90‑1 (MQ=255) ggCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg < 2:212922/90‑1 (MQ=255) gggggAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACa < 1:218866/90‑1 (MQ=255) gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc < 2:192220/90‑1 (MQ=255) | ATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAC > NZ_CP009273/4464430‑4464599 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |