Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,464,532 | 0 | T | G | 56.6% | 33.6 / 43.1 | 54 | *198G (TGA→GGA) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base T (6/17); new base G (30/0); total (37/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.12e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAA > NZ_CP009273/4464447‑4464617 | gAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGAAAAAAAATATTTTTTAGGGGGGATg > 1:453855/1‑90 (MQ=255) cGAGTGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCt < 1:297262/90‑1 (MQ=255) cGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCt < 1:399311/90‑1 (MQ=255) gCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGGGGGATGCGGGCCTTGACCCG‑a > 2:201003/1‑90 (MQ=255) gCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑a < 2:562032/90‑1 (MQ=255) aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTGACCCG‑aaa > 2:93153/1‑90 (MQ=255) acaacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTGAGGGGGGATGCGGGCCCTGACCCG‑AAAGGCGGGaaaa > 1:391572/1‑89 (MQ=255) caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTGGGGGGGATGCGGGGCTTGACCCG‑AAAGGCGGGgaaac > 1:376763/1‑88 (MQ=255) caacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGaaaac > 2:279850/1‑88 (MQ=255) acccctGTGCAACAGCTGCAATGGGGGGGGAAGCGAAAAAAACTATTCTTGAGGAGTGTTGCGGGCCTTGCCCCG‑AAGGGGGGGAAAgcg < 1:516428/86‑1 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAAGGGGGAAAgggc > 2:325219/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGGGGGATGCGGGCCTTGACCCG‑AAAGGCGGGGAAgcgc > 1:587378/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGGGAGATGCGGGCCTTGGCCCG‑CAAGGCGGGGAAgcgc > 1:369224/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTTACCCG‑AAAGGGGGGAAAgcgc > 2:583984/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTTACCCG‑AAAGGCGGGGAAacgc > 2:524413/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAacgc > 2:482385/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAgagc > 1:331996/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCCTGACCCG‑AAAGGCGGGgaaaagc > 1:97529/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGCCTTGACCCG‑AAAGGCGGGGAAgcgc > 1:92419/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTTTTTGGGGGGGATGGGGGCCTTGGCCCG‑AAAGGGGGGgaaaagc > 1:491865/1‑86 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGAAAAAAAATATTTTTGAGGAGGGATGGGGGCCTTGACCCGAAAAGGGGGGaaacgg > 2:863/1‑87 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGGTAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAgcgc > 2:126603/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTTACCCC‑AAAGGCGGGGAAacgc > 2:89099/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCC‑AAAGGGGGGAAAgcgc > 2:431255/1‑90 (MQ=255) acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTTAGGAGTGGTGTGGGGCTTGACCCG‑AAAGGCGGGgaaaagc > 1:435969/1‑86 (MQ=255) ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGcc < 2:470295/90‑1 (MQ=255) ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGcc < 2:319666/90‑1 (MQ=255) ctctGCGCAAAAGCTGAAATGGCGGGGGAAGCAAAAAAACCTATTCTTGGGGGGTGAGGGGGGCCTTGCCCCG‑AAAGGCGGGAACGCGcc < 1:561762/90‑1 (MQ=255) tctGTGCAACAGCTGCAATGGGGGGGGAAGCGAAAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCt < 1:186156/90‑1 (MQ=255) tctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGGCTTTACCCG‑AAAGGCGGGAAAGCGCCt > 2:209356/1‑90 (MQ=255) tctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTTTTTGGGGGGGATGCGGGGCTTTACCCG‑AAAGGGGGGGAAAGGCCt > 1:534580/1‑90 (MQ=255) ctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGAGGGATGCGGGGCCTGACCCG‑AAAGGGGGGAAAGCGcccg > 1:404200/1‑88 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTGGGAGGGATGGGGGCCTTGACCCG‑AAAGGCGGGAAAAcgcccgcc > 1:319666/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGAAGCGGGCCTTGACCCG‑AAAAGGGGGGAAACGccctcc > 1:603310/1‑86 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGGGGGGAAAGGCCCGcc > 2:286830/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAACGCCTGcc > 2:516428/1‑90 (MQ=255) gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGAGGGGGGATGGGGGGCCTTACCCG‑AAAGGCGGGGAAAcgcccgcc > 1:296919/1‑90 (MQ=255) aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATa < 1:524413/90‑1 (MQ=255) aaCAGCTGAAATGGCGGGGGAAGCGAAAAAAACTATTTTTGAGGAGTGATGCGGGTCTTGCCCCA‑AAAGGCGGGAAAGCGCCTGCCGATa < 1:4573/90‑1 (MQ=255) aCAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGGGGAGGGATGCGGGCCTTGACCCC‑AAAGGCGGGAAAGCGCCCGCCGGTaa > 2:584931/1‑90 (MQ=255) cAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGATGCGGGCCTTGACCCG‑AAAGGCGGGGAAACGGCCGCCGATAAg > 1:470295/1‑90 (MQ=255) cAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAg > 2:218715/1‑90 (MQ=255) aTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATa < 2:181818/90‑1 (MQ=255) tGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATaa < 2:603310/90‑1 (MQ=255) ggggggAGGGATAAAAACTATTTTTGAGGGGGGATGCGGGCCTTGACCCG‑AAAGGGGGGAAAGCGCCGGCCGGTAAGTTGTGGTAAAACa > 1:250852/1‑90 (MQ=255) gggggAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACa > 2:320210/1‑90 (MQ=255) gCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTc > 1:583405/1‑90 (MQ=255) gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc < 2:453855/90‑1 (MQ=255) ttCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCt > 2:208144/1‑90 (MQ=255) ttGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGc < 2:587378/90‑1 (MQ=255) tggggagTGATGCGGGCTTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCg < 1:583984/87‑1 (MQ=255) tGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCg > 2:637784/1‑90 (MQ=255) tGAGGAGTGATGCGGGCCTTGACCCA‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCg < 1:637784/90‑1 (MQ=255) ggagTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGaa < 2:514264/90‑1 (MQ=255) | GAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCG‑AAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAA > NZ_CP009273/4464447‑4464617 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |