Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,321 | 0 | T | G | 57.1% | 22.9 / 33.0 | 28 | T6P (ACC→CCC) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base T (3/9); new base G (16/0); total (19/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGA > NZ_CP009273/4558232‑4558404 | cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGt < 2:464409/90‑1 (MQ=255) cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaa > 1:212127/1‑90 (MQ=255) cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaa > 1:374483/1‑90 (MQ=255) ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 2:604385/1‑90 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 2:448255/1‑89 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 2:290637/1‑89 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaaaaaag > 2:555108/1‑88 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaaa > 1:574686/1‑88 (MQ=255) aTCAGATCCGTCGACAGATACTCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGt < 1:495879/90‑1 (MQ=255) cAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg < 2:614128/90‑1 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaaa > 2:634681/1‑89 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGGATaaaa > 2:409481/1‑88 (MQ=255) aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca < 1:409481/90‑1 (MQ=255) cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa < 1:333508/90‑1 (MQ=255) cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAaaaaaaaa > 1:704758/1‑90 (MQ=255) aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCCAAAAAAAAAACTCCt > 2:604499/1‑90 (MQ=255) aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACCGGGCATAAAAAACAACACCCTTTcccg > 2:521722/1‑88 (MQ=255) aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCAAAAAAAAAAAACCCTTTTcccg > 1:687496/1‑88 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACTCCTTTTCCCGGGAACCCCAACgg > 2:334998/1‑90 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCCTGTCCGGAGAACCCCCACggggga > 2:367192/1‑88 (MQ=255) cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACCCCTTTTCCGGGGAACCCCCACGGGGgaa > 2:481353/1‑90 (MQ=255) caggggAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:485255/87‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 1:363031/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:105725/90‑1 (MQ=255) cGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGaaa < 1:555108/90‑1 (MQ=255) cAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACTCCTTTTCCGGGGAAACCCCACGGGGGAAAAGAAAATCaaaaac > 2:194270/1‑86 (MQ=255) cGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACTCCTTGTCCGGGGAACCCCCAACGGGAGAACGAAAGTCAATAACgac > 1:727140/1‑88 (MQ=255) cGTGGCGGCATGGCGGGGAAAAAAAAGGGGAGGAACAAACAAAATCTTGTACCGGGAAACCCCACGGGGGAAAAGAAAATTAATACCgcg > 2:582204/1‑88 (MQ=255) ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGa < 1:367115/90‑1 (MQ=255) | CCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGA > NZ_CP009273/4558232‑4558404 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |