Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,321 | 0 | T | G | 59.1% | 18.0 / 22.4 | 22 | T6P (ACC→CCC) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (13/0); total (14/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC > NZ_CP009273/4558236‑4558409 | ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 1:135672/1‑90 (MQ=255) ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 2:107131/1‑90 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 2:410588/1‑89 (MQ=255) gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTggga > 2:285966/1‑88 (MQ=255) ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa < 2:462016/90‑1 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATaaa > 2:456062/1‑89 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATaaaa > 2:235053/1‑88 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaaaa > 2:87312/1‑88 (MQ=255) gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaaaa > 2:250051/1‑88 (MQ=255) caggATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca < 1:238197/87‑1 (MQ=255) aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCt < 2:142207/90‑1 (MQ=255) aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTcctg > 2:238197/1‑88 (MQ=255) aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAACTCCTTGTCCGga > 2:261380/1‑90 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAAAAAAAAACACCCTGTCCGGGGAAcc > 2:638320/1‑90 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAAAACCTTTTTCGGGGaagc > 1:617238/1‑88 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAAAACCCTTTCCCGGGAAcc > 1:459466/1‑90 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATAAAAAAAAAAACCCTTTCCGGAGAAACCCCACggggga > 1:166290/1‑88 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 2:120392/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa < 1:61120/90‑1 (MQ=255) aCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATAc < 2:158322/90‑1 (MQ=255) cggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa < 2:723472/90‑1 (MQ=255) gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc < 1:456062/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 < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |