Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,292,475 | 0 | T | G | 51.7% | 35.1 / 38.2 | 29 | T141P (ACC→CCC) | mdtO | multidrug efflux transporter permease subunit MdtO |
Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (15/0); total (18/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292404‑4292552 | gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGGTGCCCCGAAATAACGAc > 2:249374/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAGCGACaa > 1:394037/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACaa > 1:124241/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCGAAATAGCGACaa > 1:103735/1‑90 (MQ=255) aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGAc < 1:309871/90‑1 (MQ=255) aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGAc < 1:384215/90‑1 (MQ=255) gATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAACGACAAAGGGggcggc > 2:55805/1‑90 (MQ=255) gATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAACGACAAAGGGGGCggg > 1:59195/1‑89 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACAAAGGGGGCGGC‑Ggaaaa > 2:108974/1‑90 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCCAAAAAACGACAAAGGGGGCGGC‑Ggaaaa > 2:304933/1‑90 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGcc < 2:29750/57‑1 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGcc > 1:29750/1‑57 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAATAGCGACAATGGCGACGGC‑GGAAAAGa > 1:58332/1‑90 (MQ=255) caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCCAAAAAACGACAAAGGCGGCGGG‑Ggaaaaaaa > 2:296328/1‑87 (MQ=255) caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGAc < 1:174556/90‑1 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGAcc < 1:317442/90‑1 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGACAAAGGCGGCGGC‑GGAAAAAACCa > 2:174556/1‑90 (MQ=255) ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGcc < 1:351356/90‑1 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGcccc < 1:296328/90‑1 (MQ=255) gcgcACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAAAAGCGACAAAGGCGGCGGC‑GGAAAAAACCAAccccca > 2:128760/1‑88 (MQ=255) gcgcACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACAAAGGCGGCGGC‑GGAAAAAAACAGccccca > 2:317442/1‑88 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAACGCCAAAGGCGGCGGCGGAAAAAATCCACCCCCaggg > 2:20829/1‑87 (MQ=255) aCGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAAAAACGACAAAGGCGGCGGG‑GGAAAAAACCAACCCCCAGCGa > 2:28551/1‑90 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTgc < 1:28551/90‑1 (MQ=255) gTCGAGCATGGCGGGGGAGGTTTTCCCGTAAATAGCGACAAAGGCGACGGC‑GAAAAAAACCAGCCCCAAACGATGGGGGGGCATCaaaaa > 1:220323/1‑90 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAACGACAATGGCGGCGGC‑GGAAAAAACCCGCCCCCAACGGTGGGGGCGCCTCaaaaa > 1:309591/1‑90 (MQ=255) gAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACAt < 2:309591/90‑1 (MQ=255) gCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGc < 1:332578/90‑1 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 2:67898/90‑1 (MQ=255) | GAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGC‑GAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292404‑4292552 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |