Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,339 | 0 | C | A | 61.9% | ‑3.3 / 16.6 | 21 | intergenic (‑3/‑476) | mdtM/BW25113_RS25845 | multidrug efflux MFS transporter MdtM/Rpn family recombination‑promoting nuclease/putative transposase |
Reads supporting (aligned to +/- strand): ref base C (3/5); new base A (13/0); total (16/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.75e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAG > NZ_CP009273/4558251‑4558428 | gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaaca < 1:269560/90‑1 (MQ=255) 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) | GTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAG > NZ_CP009273/4558251‑4558428 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |