Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,292,475 | 0 | T | G | 55.0% | 15.1 / 17.7 | 20 | T141P (ACC→CCC) | mdtO | multidrug efflux transporter permease subunit MdtO |
Reads supporting (aligned to +/- strand): ref base T (3/6); new base G (11/0); total (14/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.32e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292391‑4292552 | gCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTGc > 2:309005/1‑90 (MQ=255) cAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGccc < 2:53468/90‑1 (MQ=255) aaGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTGCCCg > 2:110616/1‑90 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTGGCCGGAAAAAACGAc > 2:440528/1‑90 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGGTTTCCCGTAAAAAGCGAc > 1:186241/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTTCCCCTAAATAGCGACaa > 1:350551/1‑90 (MQ=255) aacaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGGTGGCCGTAAAAAACGACAAAGGc > 2:445372/1‑90 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc < 2:21022/90‑1 (MQ=255) aaGCGCACGACCGCTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc < 2:314020/90‑1 (MQ=255) gcgcACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCaa < 1:44496/90‑1 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTTCCCCTAAATAACGACAAAGGCGGCGGCGAAAAAAAACAGCCCCCAGCg > 1:418142/1‑90 (MQ=255) cGACCACTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTGCCCGGTAAAAAAGACAATGGGGACGGCGAAAAAAAACCGCCCCCAggggt > 1:335053/1‑86 (MQ=255) aCTTCCGGATAGTCGAGCATGGC‑GGGGGAGGGTTTCCCGTAAAAAACGACAATGGCGGCGGCGAAAAAAAACAGCCCCAAGCGATggggg > 1:428412/1‑88 (MQ=255) gTCGAGCATGGC‑GGGGGGGGGTTTCCCGGAAAAAACGACAAAGGGGGCCGCGGAAAAAACCAACCCCCAACGAAGGGGGCGGCTCaaaaa > 2:413503/1‑90 (MQ=255) gTCGAGCATGGC‑GGGGGAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGGAAAAAAACAGCCCCCAGCGATGGGGGGGCATCaaaaa > 2:242458/1‑90 (MQ=255) gTCGAGCATGGC‑GGGGGAGGGTTTCCCGGAAAAAACGACAAAGGGGGCGGCGGAAAAAAACAGCCCCCAACGATGGGGGGGCCTCaaaaa > 2:398410/1‑90 (MQ=255) gTCGAGCATGGC‑GGGGGAGGGTTGCCCGTAAAAAGCGACAAAGGCGGCGGCGAAAAAAACCAGCCCCCAACGATGGGGGGGCCTCaaaaa > 1:255727/1‑90 (MQ=255) cGAGCATGGC‑GGGGGAGGTGTTCCCCTAAAAAGCGACAAAGGCGGTGGCGGAAAAGAACCGCCCCCAGCGGTGGGGGCGCCTGaaaaaaa > 1:448786/1‑88 (MQ=255) gCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGc < 2:74649/90‑1 (MQ=255) tGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 1:419151/90‑1 (MQ=255) | GCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGC‑GGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292391‑4292552 |
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: ‑ |