Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,321 | 0 | T | G | 53.1% | 38.0 / 46.3 | 32 | T6P (ACC→CCC) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base T (1/14); new base G (17/0); total (18/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.18e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC > NZ_CP009273/4558235‑4558409 | aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaa > 2:491284/1‑90 (MQ=255) ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 1:31202/1‑90 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 2:68365/1‑89 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaag > 2:181944/1‑88 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaaa > 2:439630/1‑88 (MQ=255) aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGt < 2:380316/90‑1 (MQ=255) aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCAt < 1:660422/90‑1 (MQ=255) aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCAt < 2:416056/90‑1 (MQ=255) ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa < 2:280715/90‑1 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATaac > 2:419368/1‑90 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaaaaaa > 2:631253/1‑78 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAc < 1:631253/79‑1 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaaca > 2:42202/1‑90 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGGATaaaa > 1:362726/1‑88 (MQ=255) cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAACGGCATAaaaaaaaa > 1:145425/1‑90 (MQ=255) gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAacaaaaaaac > 2:811161/1‑88 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGGATAAAAAAAAACACCCTTTCCCGGGaaac > 1:393912/1‑88 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGGATAAAAAAAAACCCCCTGTCCGGGGAACCCCAACgg > 1:289237/1‑90 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAACTCCTTGTCCGGAGAACCCCCACcgggga > 2:741849/1‑88 (MQ=255) cATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACCCCCTTTCCGGGGAACCCCCACGGGGgaaaag > 2:481415/1‑88 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:440807/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:456408/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:73488/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 1:256904/90‑1 (MQ=255) aaaaaCAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTTTCCGGGGAACCCCCACGGGGGAAAAGAAAATCa > 1:330595/1‑90 (MQ=255) aaaCAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAAt > 2:272891/1‑90 (MQ=255) gCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAAATCAATACCGt > 2:590963/1‑90 (MQ=255) cggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa < 2:362726/90‑1 (MQ=255) ggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAAt < 2:420474/90‑1 (MQ=255) ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGa < 1:644292/90‑1 (MQ=255) cGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa < 2:804131/90‑1 (MQ=255) gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc < 2:32307/90‑1 (MQ=255) | AGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC > NZ_CP009273/4558235‑4558409 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |