| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,160,604 | 0 | A | C | 28.0% | 47.8 / 10.5 | 25 | Y740S (TAC→TCC) | mdtC | multidrug efflux pump RND permease subunit MdtC |
| Reads supporting (aligned to +/- strand): ref base A (18/0); new base C (0/7); total (18/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.08e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCG > NC_000913/2160537‑2160671 | aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:1004571/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:240539/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:215845/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:2727718/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:3019697/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:3338143/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:394133/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:726537/1‑68 (MQ=255)aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa > 1:1274779/1‑68 (MQ=255) cGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGCCAGCGGCAAATCTCGACCATTTACCAGc > 1:400681/1‑68 (MQ=255) aaGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACAAGCCGAta > 1:2332350/1‑67 (MQ=255) aGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGa > 1:1608923/1‑68 (MQ=255) ccgccAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATCAAc > 1:3302159/1‑68 (MQ=255) ccAACAGTCTTTTAAATACCCCCTTCGGTCAGCGGCAAATCTCGACCTTTTCCCAGCCGATGAACCAg < 1:2368572/68‑1 (MQ=255) cAACGGTCTGTTAAATACCGCCGGCGGTCAGGGGAAATCCTCGGCCATTTCCCAGCCGATGAACCAGt < 1:1090801/68‑1 (MQ=255) cAACATTCTGTTAAATCCCCCCTGCGGCCAGCGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGt < 1:1478985/68‑1 (MQ=255) cAACATTCTGTTAAATAACGCTTGCGGGCGGGGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGt < 1:2537333/68‑1 (MQ=255) aCAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGcat > 1:2667945/1‑68 (MQ=255) aCAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATATCGACCATTTACCAGCCGATGAACCAGtat > 1:2765554/1‑68 (MQ=255) cAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGtata > 1:1923195/1‑68 (MQ=255) aGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATaa > 1:278907/1‑68 (MQ=255) cTGTTAAATACCGCCTTCGGTCAGCGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGTATAAAgt < 1:1823714/68‑1 (MQ=255) cTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgt < 1:745605/68‑1 (MQ=255) aTAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGCTg > 1:3154704/1‑68 (MQ=255) aTAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATg > 1:255044/1‑68 (MQ=255) aTAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATg > 1:2336502/1‑68 (MQ=255) aTAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATg > 1:3169522/1‑68 (MQ=255) aaCGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATgct > 1:1020865/1‑66 (MQ=255) aaCGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGa > 1:2453022/1‑68 (MQ=255) aaCGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGa > 1:1881213/1‑68 (MQ=255) aCGCCTTCGGTCGGCGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGaa < 1:464476/68‑1 (MQ=255) gCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGt > 1:3540697/1‑68 (MQ=255) gCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGt > 1:1608124/1‑68 (MQ=255) ggTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATcc > 1:3164952/1‑68 (MQ=255) tCAGCGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgc < 1:770326/68‑1 (MQ=255) tCAGCGGCAAATCTCGACCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgc < 1:1791580/68‑1 (MQ=255) tCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgc > 1:2608075/1‑68 (MQ=255) tCACCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgc < 1:1359588/68‑1 (MQ=255) gcAGCGGGAACTCTGGGCCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgc < 1:3292835/67‑1 (MQ=255) cAGCGCCAACTCGCGTCCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcgcg < 1:148744/68‑1 (MQ=255) cgCGGCAACTCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTNATGGAAATGGATCCgcgc < 1:2705294/67‑1 (MQ=255) ggCAAATCTCGCCCTTTCCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCtat < 1:1606212/68‑1 (MQ=255) ggCAAATCTCGCCCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCtat < 1:3007691/68‑1 (MQ=255) ggCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCtat < 1:328643/68‑1 (MQ=255) ggCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGCGGATCCGCGCtat > 1:2377393/1‑68 (MQ=255) ggCAAATCGCGCCCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCtat < 1:224530/68‑1 (MQ=255) cAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATAc < 1:1891404/68‑1 (MQ=255) aaaaTCTCGACCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATAc < 1:1148198/67‑1 (MQ=255) aTCGCGACAATTCCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCa < 1:1072558/68‑1 (MQ=255) ctcGCCGATTAACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAgg < 1:1532753/68‑1 (MQ=255) ctcGACCATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAgg < 1:1689078/68‑1 (MQ=255) ctcGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAgg > 1:61830/1‑68 (MQ=255) gcccATTTCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACa < 1:1560788/66‑1 (MQ=255) gACCATTCCCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACa < 1:1147997/68‑1 (MQ=255) cATTTACCAGCCGATGAACCAGTATAACGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCa > 1:3014947/1‑68 (MQ=255) aCCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTgcg > 1:1893175/1‑68 (MQ=255) | ATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCG > NC_000913/2160537‑2160671 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |