Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,292,475 | 0 | T | G | 56.5% | 26.1 / 28.0 | 23 | T141P (ACC→CCC) | mdtO | multidrug efflux transporter permease subunit MdtO |
Reads supporting (aligned to +/- strand): ref base T (2/8); new base G (13/0); total (15/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.18e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC > NZ_CP009273/4292402‑4292563 | tAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCg < 1:462780/90‑1 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGAc > 2:14338/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGGCCGTAAATAACGACaa > 2:462780/1‑90 (MQ=255) caacaaCGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCGAAAAAACGACAAAgg > 1:158286/1‑90 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGCAAAAAACGGCAAAGGCGGCGGCGgaaaa > 1:342960/1‑90 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGGAAAAAAGGACAATGGCGACGGCGaaaaaaa > 2:387687/1‑88 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCTAAAAAACGACAAAGGCGGCGGCGGAAAAGACCa > 1:282895/1‑90 (MQ=255) ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAAAAACGACAAAGGCGGCCGCGGAAAAAAACAGcc > 1:268423/1‑90 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc < 1:209905/90‑1 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAACGACAAAGGCGGCCGCGGAAAAAACCAGcccc > 2:44865/1‑90 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAATAACGGCAAAGGCGGCGGGGGAAAAAAACCGCCCCAAGCg > 1:429515/1‑90 (MQ=255) cGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGAt < 1:288585/90‑1 (MQ=255) ggATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATc > 2:245360/1‑90 (MQ=255) tAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCaaa < 2:420106/90‑1 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAAGGACAAAGGCGACGGGGAAAAAAACCAGCCCCCAACGGAGGGGGGGCGTCaaaaa > 1:250223/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGTTTGCCCCGAAAAAGCGGCAATTGCGACGGCGAAAAAAACCCGCCCCAAGGGATGGGGGCGCCTCaaaaa > 2:291419/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATCaaaaa > 2:226885/1‑90 (MQ=255) tGGCGGGGGAGGGTTGCCCGGAAAAAGCGACAAAGGCGACGGCGAAAAAAAACAGCCCCCAGCGATGGGGGCGCATCAAAAAAAagcagc > 1:25069/1‑90 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 1:14338/90‑1 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 2:214262/90‑1 (MQ=255) ggAAGGGTTGCCCGTAAAAAGCGACAAAGGCGACGGCGAAAAAGACCAGCCCCCAGCGATGGGGGGGCATCAAAAAAAAGCAGCCCATCa > 1:48947/1‑90 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 1:492843/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 1:291419/90‑1 (MQ=255) | TAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC > NZ_CP009273/4292402‑4292563 |
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: ‑ |