Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,343 | 0 | C | A | 66.7% | ‑5.7 / 16.8 | 21 | intergenic (‑7/‑472) | mdtM/BW25113_RS25845 | multidrug efflux MFS transporter MdtM/Rpn family recombination‑promoting nuclease/putative transposase |
Reads supporting (aligned to +/- strand): ref base C (2/5); new base A (14/0); total (16/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACT > NZ_CP009273/4558271‑4558431 | gTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGgaga < 2:355456/90‑1 (MQ=255) acagtACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGgaga < 2:218863/86‑1 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACTCCTTTT‑CCGGAGAACCCCCACgg > 2:103645/1‑90 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAAATCCTTTT‑CCGGAGAACCCCAACgg > 2:599376/1‑90 (MQ=255) cAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTTT‑CCGGAGAACCCCCAAgggg > 2:77443/1‑90 (MQ=255) aaCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAACACCCTTT‑CCGGGGAACCCCCACGGGGa > 1:610150/1‑90 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGCGGGCTAAAAAAAAACTCCTTGT‑CCGGGGAACCCCCAAggggga > 2:148700/1‑88 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGGATAAAAAAAAAAACCTTTT‑TCGGGGAAACCCCACggggga > 1:318289/1‑88 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATAAAAAAAAACACCTTTTCCCGGGGAACCCCACCgctgt > 1:218863/1‑86 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAAACCCTTTT‑CCGGGGAACCCCCACggggga > 1:80602/1‑88 (MQ=255) cATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGGAAAAAAAAAAAAACCCTTT‑CCCGGGAAACCCCACGGGGaaaaag > 1:258668/1‑88 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGATAAAAAAAAAAACCCTGT‑CCGGGGAACCCCAACGGGGGAAAAGa > 1:527560/1‑90 (MQ=255) aCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATAc < 1:639158/90‑1 (MQ=255) aGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCg < 2:560873/90‑1 (MQ=255) aGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCAAAAAAAAAAAATCCTTTT‑CCGGAGAACCCCCACGGGGGAAAAGAAAATCAATAccc > 2:545031/1‑89 (MQ=255) aGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGAAAACAAAAAACTCCCTTT‑CCGGAGGAACCCAACGGGGGAAAAGAAAATCAAAACCg > 1:428153/1‑90 (MQ=255) gCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAACAACTCCTTTT‑CCGGGGAACCCCAACGGGGGAAAAGAAAATCAAAACCGt > 2:98257/1‑90 (MQ=255) gCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCCTGT‑CCGGAGAACCCCAACGGGGAAAAAGAAAATTAATAcccg > 2:197154/1‑88 (MQ=255) gcggcATGGCGGGTAAAAAAAAATGGGATAACAAAAAAACCCCTTT‑CCGGGGAACCCCAACGGGGGAAAAGAAAATAAATACCGTCggga > 1:67559/1‑88 (MQ=255) gcggcATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACCCCTTTT‑CCGGGGAACCCCCACGGGGAAAAAGAAATTCAAAACCGTCggga > 1:192022/1‑88 (MQ=255) cggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa < 2:428153/90‑1 (MQ=255) aacaaacaACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACg < 1:224713/90‑1 (MQ=255) caaacaaCTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag < 1:42116/90‑1 (MQ=255) caaacaaCTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag < 2:468787/90‑1 (MQ=255) aacaacTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc > 2:43451/1‑90 (MQ=255) acaacTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACt < 2:403524/90‑1 (MQ=255) acaacTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACt < 1:481338/90‑1 (MQ=255) | GTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACT > NZ_CP009273/4558271‑4558431 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |