Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,464,532 | 0 | T | G | 53.1% | 25.4 / 29.8 | 33 | *198G (TGA→GGA) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base T (5/10); new base G (17/0); total (23/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.65e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATT > NZ_CP009273/4464450‑4464619 | ggCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCgg < 1:170420/90‑1 (MQ=255) cGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCt < 1:288797/90‑1 (MQ=255) gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGGCCTGAccc > 2:52503/1‑90 (MQ=255) ggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGAGGAGTGATGCGGGCCTTGACCCg > 1:360437/1‑90 (MQ=255) ctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGGGGGATGCGGGGCTTGACCCGAAAGGc > 2:150857/1‑90 (MQ=255) gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTTTTTTGGGGGGATGCGGGGCTTGACCCGAAAGGCgggga > 2:443958/1‑88 (MQ=255) acaacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAg > 1:4839/1‑90 (MQ=255) acaacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGgaaa > 2:61547/1‑89 (MQ=255) caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTTACCCCAAAGGCGGGaaaac > 1:389917/1‑88 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTCTTGAGGAGAGAGGCGGGCCTTGACCCGAAAGGGGGGGAAgggc > 1:374653/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTTACCCGAAAGGCGGGgaaaccc > 1:359921/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc > 2:57242/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGGGGGGAAacgc > 2:211598/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGTGGGGGGATGCGGGGCCTGGCCCGAAAGGCGGGGAAacgc > 2:321150/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGGGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 2:91135/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc > 2:204567/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt > 1:5290/1‑56 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt < 2:5290/56‑1 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc > 1:207531/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGAGGCGGGCCTTGAGCCGAAAGGCGGGAAAgcgc > 1:35060/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGACATAAAAAACATTCTTGAGGAAGGGTGTGGGTCTTAAGCCGAAAGGGGGGaaaatgg > 2:319197/1‑86 (MQ=255) tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAACGccc > 1:332166/1‑89 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGGGGGGAAGcgcccgcc > 1:448341/1‑90 (MQ=255) aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 1:7880/90‑1 (MQ=255) aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 1:241533/90‑1 (MQ=255) aCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAACGCCCGCCGATaa > 2:197709/1‑90 (MQ=255) tGCAATGGCGGGGGAAGGGATAAAAAATATTCTTGAGGGGGGATGCGGGGCTTGACCCGAAAGGGGGGGAAGCGCCTGCCGGTAAGTtgg > 2:328999/1‑89 (MQ=255) tGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa > 1:180142/1‑84 (MQ=255) tGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa < 2:180142/84‑1 (MQ=255) aagggCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGGGGGCCTTGCCCCGAAAGGGGGGAAAGCGCCTGCCGATAAGTTGTGAt < 2:339725/87‑1 (MQ=255) gCGATAAAAACTATTCTTGAGGAGTGAT‑‑GGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc < 1:491506/90‑1 (MQ=255) gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc < 1:81119/90‑1 (MQ=255) ggagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGaa < 1:149651/90‑1 (MQ=255) gagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAAt < 1:283538/90‑1 (MQ=255) agTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAAtt > 1:336457/1‑90 (MQ=255) | GGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATT > NZ_CP009273/4464450‑4464619 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |