Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,292,475 | 0 | T | G | 52.0% | 64.6 / 70.6 | 50 | T141P (ACC→CCC) | mdtO | multidrug efflux transporter permease subunit MdtO |
Reads supporting (aligned to +/- strand): ref base T (4/20); new base G (26/0); total (30/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC > NZ_CP009273/4292388‑4292563 | tcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCTAGCATGGCGGGGAAGGttt < 1:866607/90‑1 (MQ=255) ggATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGTCCCTaaaaa > 1:238186/1‑88 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGCAAAAAGCGAc > 1:528138/1‑90 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCCAAAAAGCGAc > 1:473789/1‑90 (MQ=255) gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGAc > 2:383501/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGACaa > 1:878890/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACaa > 2:287707/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACaa > 2:351622/1‑90 (MQ=255) ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACaa > 2:156117/1‑90 (MQ=255) acaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAGCGACAAAgggg > 2:143536/1‑88 (MQ=255) aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGAc < 2:138088/90‑1 (MQ=255) gATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTGCCCGTAAATAACGACAATggcggcggc > 2:678198/1‑90 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACAAAGGCGGCGGCGgaaaa > 1:315189/1‑90 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTGCCGTAAAAAAGGGCAATGGCGGCGGGGgaaaa > 2:520507/1‑90 (MQ=255) ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAAAGGCGGCGGCGgaaaa > 2:468490/1‑90 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAATAGCGACAAAGGCGGCGGCGgaaaaaa > 2:407498/1‑88 (MQ=255) acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCTAAAAAGCGACAAAGGCGACGGCGaaaaaaa > 2:809553/1‑88 (MQ=255) caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAACGACAAAGGCGACGGGGGAAAAAAc > 1:210327/1‑90 (MQ=255) caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCCAAAAAACGACAAAAGGGGCGGCGgaaaaaaa > 1:902913/1‑87 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc < 1:678198/90‑1 (MQ=255) aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc < 2:156721/90‑1 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc < 2:939718/90‑1 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCCGTAATAGGAACAAATGGGAGGGGGGAAAAGACCAAcccc > 1:751938/1‑90 (MQ=255) aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGcccc > 2:103933/1‑90 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAGCGACAAAGGGGGCGGCGGAAAAAACCAGCCCCCAGCg > 1:67853/1‑90 (MQ=255) cACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACAATGGCGACGGCGGAAAAGACCAGCCCCAAGCg > 2:125912/1‑90 (MQ=255) cGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGAt < 1:577394/90‑1 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc < 2:7231/90‑1 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc < 2:504349/90‑1 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc < 2:205022/90‑1 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc < 1:383501/90‑1 (MQ=255) cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc < 1:905992/90‑1 (MQ=255) ttCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgcg < 2:645352/90‑1 (MQ=255) aGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCaaaa < 2:473789/90‑1 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAACGACAAAGGCGGCGGGGGAAAAAACCAGCCCCCAGCGATGGGGGGGCCTCaaaaa > 1:205022/1‑90 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGGGACGGCGAAAAAAACCAGCCCCCAACGATGGGGGCGCATCaaaaa > 2:864646/1‑90 (MQ=255) gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAAAGGGGACGGCGAAAAAAACCAGCCCCCAGCGGAGGGGGGGCATCaaaaa > 1:346248/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAAACAGCCCCAAACGATGGGGGGGCATCaaaaa > 2:94972/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGATGGGGGCGCATCaaaaa > 2:682631/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGAGGGGGGCGCATCaaaaa > 2:531961/1‑90 (MQ=255) gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAACGATGGGGGGGGATCaaaaa > 2:829893/1‑90 (MQ=255) gAGCATGGCGGGGAAGGGTTGCCCGTAAATAACGACAATGGGGACGGCGAAAAAAAACAGCCCCCAACCATGGGGGCGGAGCaaaaaaaa > 2:540617/1‑87 (MQ=255) gCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGc > 2:395282/1‑90 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 1:540617/90‑1 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 2:528138/90‑1 (MQ=255) tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc < 1:459600/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 2:232180/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 2:552483/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 1:463070/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 1:287707/90‑1 (MQ=255) gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc < 1:130167/90‑1 (MQ=255) | TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC > NZ_CP009273/4292388‑4292563 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |