Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 312,225 | C→A | 31.1% | A216A (GCG→GCT) | acrB ← | multidrug efflux system protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 312,225 | 0 | C | A | 31.1% | 90.7 / 44.8 | 77 | A216A (GCG→GCT) | acrB | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base C (45/8); new base A (11/13); total (56/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.17e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-03 |
TGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGT > minE/312163‑312292 | tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCa > 1:1139815/1‑69 (MQ=255) tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCa > 1:2364843/1‑69 (MQ=255) gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc > 1:372151/1‑69 (MQ=255) gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc > 1:2975600/1‑69 (MQ=255) gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc > 1:725117/1‑69 (MQ=255) gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc > 1:2028324/1‑69 (MQ=255) ataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCaa > 1:2658836/1‑61 (MQ=255) gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc > 1:430561/1‑69 (MQ=255) gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc > 1:3264492/1‑69 (MQ=255) gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc > 1:2786224/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCt > 1:2250333/1‑68 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:2259199/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:2164412/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:1431834/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:3295229/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:2181951/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:835399/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:2696083/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:484426/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:1320212/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:2878777/1‑69 (MQ=255) ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg > 1:1499852/1‑69 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTTACCCGCCGCAACCTGGGCGTTCTGCGCttt < 1:809402/69‑1 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCACCCTGGGCGTTCTGCGCttt > 1:1705023/1‑69 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt < 1:2277039/69‑1 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt > 1:1302374/1‑69 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt < 1:2255684/69‑1 (MQ=255) aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt > 1:2232589/1‑69 (MQ=255) tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg > 1:3118144/1‑69 (MQ=255) tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg > 1:1897063/1‑69 (MQ=255) tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg > 1:600790/1‑69 (MQ=255) tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg > 1:3252444/1‑69 (MQ=255) gCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt > 1:3324381/1‑69 (MQ=255) gCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt > 1:2686746/1‑69 (MQ=255) tttCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAAt > 1:2842782/1‑69 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTc > 1:2855617/1‑45 (MQ=38) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:2105950/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:889310/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:676711/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:1694619/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:696110/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:2452437/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:1816651/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:1844992/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:1248518/1‑68 (MQ=255) aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc > 1:3344012/1‑68 (MQ=255) ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa > 1:1422177/1‑69 (MQ=255) ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa > 1:2179174/1‑69 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:3185091/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:1478262/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:603127/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:101998/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:2941744/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:1410470/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:2375443/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:2178183/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:2112361/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:145919/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:1957006/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:1604366/53‑1 (MQ=255) tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa < 1:863649/53‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt < 1:3017226/50‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt < 1:2763983/50‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt < 1:1736255/50‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt < 1:3425509/50‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt < 1:1647967/50‑1 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTaa > 1:3361060/1‑52 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:89342/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:366514/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:3324596/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:3138870/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:2531617/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:216625/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:1992864/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt > 1:1114657/1‑69 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGt > 1:1848539/1‑68 (MQ=255) ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGt > 1:1507303/1‑68 (MQ=255) ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt < 1:455265/69‑1 (MQ=255) ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt < 1:3401244/69‑1 (MQ=255) ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt < 1:1071213/69‑1 (MQ=255) | TGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGT > minE/312163‑312292 |
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: ‑ |