Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,321 | 0 | T | G | 56.4% | 35.7 / 46.5 | 39 | T6P (ACC→CCC) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base T (2/15); new base G (22/0); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.41e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC > NZ_CP009273/4558236‑4558409 | ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 2:121686/1‑90 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa > 2:401911/1‑89 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa > 1:621581/1‑89 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 2:279793/1‑89 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaag > 1:201348/1‑88 (MQ=255) gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaa > 2:629820/1‑79 (MQ=255) gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaa < 1:629820/79‑1 (MQ=255) gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTggga > 1:33427/1‑88 (MQ=255) gATCCGTCCACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg < 2:202596/88‑1 (MQ=255) gATCCGTCCACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTgg > 1:202596/1‑88 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaac > 1:162030/1‑90 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATaaa > 1:343603/1‑89 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATaaa > 2:501735/1‑89 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATaaaa > 1:56462/1‑88 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATaaaa > 2:340593/1‑88 (MQ=255) cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa < 1:340593/90‑1 (MQ=255) aGATACGCAGCAAAGTCATACAAGATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAaacaac < 2:276730/90‑1 (MQ=255) gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCGTAaaaaaaaaaa > 1:493799/1‑88 (MQ=255) gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCAaaaaaaaaaaaa > 2:371844/1‑88 (MQ=255) gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACACCTTTTcc > 2:270565/1‑90 (MQ=255) gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGGATAAAAAAAAACCCCCTTTcc > 2:83735/1‑90 (MQ=255) aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAAACCCTTTCCgg > 2:286082/1‑90 (MQ=255) aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGATAAAAAAAAAAACCCTTTCCGga > 2:7569/1‑90 (MQ=255) cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAc < 1:121686/90‑1 (MQ=255) aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc < 2:615985/90‑1 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg < 2:61055/53‑1 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTgg > 1:61055/1‑53 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGGATAAAAAAAAACTCCTTTTCCGGAGAACCCCAACgg > 1:132377/1‑90 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAAACCCTTTCCCGGGGACCCCCACgg > 2:436670/1‑90 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACggggga > 1:172424/1‑88 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAAAAACTCCTTTTCCGGGGGACCCCCACggggga > 1:506062/1‑88 (MQ=255) aTCGGGAAAAACAGCGTGGCTGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 1:691453/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:518179/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:544421/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 1:286082/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:33427/90‑1 (MQ=255) aaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATa < 2:599788/90‑1 (MQ=255) gggTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg < 1:114454/90‑1 (MQ=255) gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc < 1:220411/90‑1 (MQ=255) | GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC > NZ_CP009273/4558236‑4558409 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |