| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,292,475 | 0 | T | G | 53.6% | 22.2 / 27.2 | 28 | T141P (ACC→CCC) | mdtO | multidrug efflux transporter permease subunit MdtO |
| Reads supporting (aligned to +/- strand): ref base T (4/9); new base G (15/0); total (19/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.04e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292399‑4292552 | ggATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATa < 1:333473/90‑1 (MQ=255) gATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAg < 2:539342/90‑1 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGGGGGGTTGCCGTAAAAATCGACaa > 1:96872/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAAAAGCGACaa > 2:67387/1‑90 (MQ=255) gATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGc < 1:229000/90‑1 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAAAGGCGACGGCGaaaaa > 2:103209/1‑90 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAGCGGCAATGGGGACGGCGgaaaaaa > 2:229000/1‑88 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAAAGACAAAGGGGGCGGGGaaaaaaa > 1:128400/1‑88 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTTAAAAAGGGCAATGGCGGCGGCGgaaaaaa > 2:457499/1‑88 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGACAAAGGGGACGGCGaaaaaaa > 1:159180/1‑88 (MQ=255) caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTTCCCCTAAAAAGCGACAAAGGCGGCGGCGAAAAAGaa > 2:365944/1‑89 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc < 2:435618/90‑1 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAAcc > 2:26425/1‑90 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAGCGACAAAGGGGACGGCGGAAAAAACCa > 2:255982/1‑90 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGACAAAGGCGGCGGGGAAAAAGACCa > 1:506284/1‑90 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCa > 2:41039/1‑90 (MQ=255) gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAAAACGACAAAGGCGGCGGGGAAAAAAACCa > 2:156344/1‑90 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGGGGGTTGCCCCTAAAAAACGACAAAGGGGGCGGGGGAAAAAAACAGCCCCCAGCg > 1:248005/1‑90 (MQ=255) gcccACTCCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATg < 2:273312/88‑1 (MQ=255) tCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgcgc < 1:451211/90‑1 (MQ=255) ggATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATc < 1:28142/90‑1 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAACGGCAAAGGGGGCGGGGGAAAAGAACAGCCCCCAACGATGGGGGGGCATCaaaaa > 1:3395/1‑90 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTTCCCCTAAATAGCGACAATGGGGGCGGCGAAAAAAACCAGCCCCCAGCGATGGGGGGGCCTCaaaaa > 2:166270/1‑90 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAGCGACAAAGGGGGCGGCGAAAAAAACCAGCCCCCAACGAAGGGGGCGCATCaaaaa > 1:105398/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATCaaaaa > 2:102912/1‑90 (MQ=255) gAGCATGGCGGGGAAGGGTTGCCCGTAAAAAACGACAAAGGCGGCGGCGGAAAAAACCAACCCCCGGCGAAGGGGGCGCATCaaaaaaaa > 1:117615/1‑87 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 1:301792/90‑1 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 2:96872/90‑1 (MQ=255) | GGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC > NZ_CP009273/4292399‑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: ‑ |