Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 312,513 | T→C | 7.5% | Q120Q (CAA→CAG) | acrB ← | multidrug efflux system protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 312,513 | 0 | T | C | 7.5% | 191.8 / 6.4 | 67 | Q120Q (CAA→CAG) | acrB | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base T (32/30); new base C (3/2); total (35/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01 |
GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTA > minE/312444‑312582 | gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt > 1:3118476/1‑71 (MQ=255) aaCGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCg > 1:3789940/1‑71 (MQ=255) cGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg > 1:631382/1‑70 (MQ=255) cGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgc < 1:3168788/71‑1 (MQ=255) gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg > 1:1275268/1‑69 (MQ=255) cGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc < 1:2261467/71‑1 (MQ=255) aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa > 1:2210745/1‑71 (MQ=255) aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa > 1:3850119/1‑71 (MQ=255) aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa > 1:2944373/1‑71 (MQ=255) aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa < 1:3189533/71‑1 (MQ=255) aCCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACgg > 1:2683621/1‑71 (MQ=255) aCCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACgg > 1:3066913/1‑71 (MQ=255) ccATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGc < 1:1405669/71‑1 (MQ=255) cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa < 1:129090/71‑1 (MQ=255) cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa > 1:3013521/1‑71 (MQ=255) cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa < 1:3217671/71‑1 (MQ=255) tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc < 1:4006811/71‑1 (MQ=255) tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc < 1:339736/71‑1 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:3920429/1‑69 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:2771831/1‑69 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:2661605/1‑69 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:2656583/1‑69 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:3842665/1‑69 (MQ=255) gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc > 1:1051993/1‑69 (MQ=255) aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGccct > 1:2937306/1‑68 (MQ=255) aGCTGCAGGATGATTTCTCAACGATCACCCCTTGCTGCTGAActtctt > 1:1505627/1‑48 (MQ=255) ctgctgTATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg < 1:2824893/71‑1 (MQ=255) ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg < 1:1574427/71‑1 (MQ=255) tgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc < 1:1488798/53‑1 (MQ=255) tgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc < 1:3339193/53‑1 (MQ=255) tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg < 1:2469988/70‑1 (MQ=255) tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg < 1:46484/70‑1 (MQ=255) tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt < 1:2443771/71‑1 (MQ=255) tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt < 1:263968/71‑1 (MQ=255) atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt < 1:3458726/70‑1 (MQ=255) tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt > 1:1570324/1‑49 (MQ=255) tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg < 1:3633194/70‑1 (MQ=255) ctcAACGCTCACCCCTTGCTGCTGGACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa > 1:3768052/1‑71 (MQ=255) ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc < 1:1285477/70‑1 (MQ=255) ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc < 1:2309462/70‑1 (MQ=255) ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa > 1:3244722/1‑71 (MQ=255) aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt > 1:737035/1‑71 (MQ=255) aCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTg < 1:2484024/71‑1 (MQ=255) cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg > 1:1001038/1‑52 (MQ=255) cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa < 1:3679459/71‑1 (MQ=255) gatccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCAcc > 1:2370368/4‑64 (MQ=255) cACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt > 1:2061855/1‑57 (MQ=255) ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc < 1:3909276/71‑1 (MQ=255) ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc > 1:2161056/1‑71 (MQ=255) ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc > 1:3584377/1‑71 (MQ=255) cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc > 1:144054/1‑70 (MQ=255) ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt < 1:2385447/48‑1 (MQ=255) tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg < 1:983296/48‑1 (MQ=255) tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTg < 1:2691196/48‑1 (MQ=19) gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt > 1:2949376/1‑51 (MQ=255) gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat < 1:1439362/70‑1 (MQ=255) gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata > 1:4121072/1‑71 (MQ=255) gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt > 1:4055741/1‑44 (MQ=16) gctgctGAACATCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt > 1:1610271/1‑44 (MQ=16) ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGat < 1:3652408/69‑1 (MQ=255) gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc > 1:3016428/1‑71 (MQ=255) ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc > 1:3629191/1‑57 (MQ=255) ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat < 1:882196/66‑1 (MQ=255) ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg > 1:2171676/1‑71 (MQ=255) gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg < 1:510482/58‑1 (MQ=255) aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc < 1:1552607/69‑1 (MQ=255) aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt < 1:3180559/71‑1 (MQ=255) aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAc > 1:3997405/1‑47 (MQ=255) aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt > 1:4026884/1‑70 (MQ=255) cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat > 1:1272280/1‑61 (MQ=255) cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat > 1:1338591/1‑61 (MQ=255) cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat > 1:2182012/1‑61 (MQ=255) ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGtt < 1:3660103/35‑1 (MQ=255) tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg > 1:35390/1‑70 (MQ=255) cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc > 1:1235457/1‑49 (MQ=255) cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc > 1:1240221/1‑49 (MQ=255) cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa > 1:92574/1‑71 (MQ=255) | GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTA > minE/312444‑312582 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |