Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 311,457 | A→G | 10.6% | I472I (ATT→ATC) | acrB ← | multidrug efflux system protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 311,457 | 0 | A | G | 10.6% | 230.8 / 17.3 | 85 | I472I (ATT→ATC) | acrB | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base A (33/43); new base G (3/6); total (36/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01 |
GTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA > minE/311389‑311525 | gtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa < 1:3944650/70‑1 (MQ=255) gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa > 1:2363516/1‑69 (MQ=255) gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa > 1:1001928/1‑69 (MQ=255) cacaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa < 1:2640482/66‑1 (MQ=255) gcgCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa > 1:181615/3‑71 (MQ=255) agCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc < 1:1056014/71‑1 (MQ=255) gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa > 1:342043/1‑65 (MQ=255) gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag < 1:115343/71‑1 (MQ=255) gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat < 1:3740230/69‑1 (MQ=255) aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag > 1:362022/1‑70 (MQ=255) aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag > 1:3684563/1‑70 (MQ=255) gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa < 1:2757028/49‑1 (MQ=255) tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATaga > 1:425825/1‑49 (MQ=255) tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc > 1:1124834/1‑69 (MQ=255) tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc > 1:3731864/1‑69 (MQ=255) tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa < 1:3324/71‑1 (MQ=255) cAATGCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga < 1:324753/50‑1 (MQ=255) cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga > 1:941713/1‑58 (MQ=255) aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAg > 1:1952394/4‑70 (MQ=9) aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGAGATAAATTGCCCCAg > 1:239808/4‑70 (MQ=9) aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag < 1:1619759/48‑1 (MQ=255) aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg < 1:3192083/70‑1 (MQ=255) aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACg < 1:2226893/55‑1 (MQ=255) tgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga < 1:3478613/46‑1 (MQ=19) cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg < 1:2575960/71‑1 (MQ=255) gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa > 1:398186/1‑66 (MQ=255) gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTa > 1:3591795/1‑69 (MQ=255) ccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag < 1:1045963/41‑1 (MQ=255) cagaacagaaAGCGCCATTGCCGAAACGATGGTGATAGAGAAc > 1:76246/1‑43 (MQ=16) caGTACCGACAGCGCCATTGCTGAAACAATGGTAATATAGAACTGACgat < 1:2672565/50‑1 (MQ=255) caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat < 1:2252326/50‑1 (MQ=255) caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga > 1:976084/1‑49 (MQ=255) aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc > 1:1847989/1‑46 (MQ=255) aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc > 1:3109767/1‑46 (MQ=255) gTACCGACAGCGCCATTGCTGAAACAATGGTAATag < 1:3117052/36‑1 (MQ=255) gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa < 1:798959/71‑1 (MQ=255) tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc < 1:742265/44‑1 (MQ=255) tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg < 1:1599318/62‑1 (MQ=255) tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg < 1:3384108/62‑1 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGCGAACt > 1:1596407/1‑39 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga > 1:3956652/1‑44 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCg > 1:2230114/1‑65 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa > 1:4080458/1‑68 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag < 1:2055908/71‑1 (MQ=255) ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag < 1:4008394/71‑1 (MQ=255) cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat < 1:3272867/44‑1 (MQ=255) aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag > 1:1475278/1‑71 (MQ=255) aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa > 1:942825/1‑70 (MQ=255) cgcCATTGCTGAAACAATTGTAATAGAGAACTGACGATAGATAGCACCa < 1:434122/49‑1 (MQ=255) gcCATTGCTGAAACAATGGTAATAGAGAACTGACgata < 1:2975359/38‑1 (MQ=255) gcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag < 1:3311858/63‑1 (MQ=12) ccATTGCTGAAACAATGGTAATAGAGAACTGACgata < 1:2301310/37‑1 (MQ=255) aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa < 1:217482/59‑1 (MQ=255) ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg < 1:2775942/48‑1 (MQ=255) ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg > 1:1240749/1‑70 (MQ=255) ttGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAAcc < 1:2349282/52‑1 (MQ=12) ttGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag < 1:2683586/59‑1 (MQ=9) tGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa < 1:4143381/57‑1 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATATCACCAg < 1:5109/43‑1 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg < 1:3609478/43‑1 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg < 1:3441354/43‑1 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa > 1:2858793/1‑48 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa > 1:3383940/1‑48 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc > 1:3593586/1‑71 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc > 1:3313949/1‑71 (MQ=255) gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc > 1:884535/1‑71 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc > 1:1052156/1‑70 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc > 1:1075745/1‑70 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg < 1:2634234/71‑1 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg < 1:2696760/71‑1 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg < 1:3604915/71‑1 (MQ=255) cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg < 1:445949/71‑1 (MQ=255) gAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaaga < 1:612709/55‑1 (MQ=12) aaaCAATGGTAATAGAGAACTGAGGATAGATAGCAc < 1:1026948/36‑1 (MQ=255) aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa < 1:3294276/70‑1 (MQ=255) aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa > 1:1214635/1‑70 (MQ=255) aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa > 1:2967883/1‑70 (MQ=255) aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt < 1:456851/71‑1 (MQ=255) aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt < 1:726494/71‑1 (MQ=255) aaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa < 1:3808073/49‑1 (MQ=255) aCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg < 1:928449/66‑1 (MQ=255) cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc < 1:603014/71‑1 (MQ=255) cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc < 1:2620980/71‑1 (MQ=255) tGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc < 1:1802132/71‑1 (MQ=255) gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc > 1:760279/1‑37 (MQ=255) gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc > 1:3733192/1‑37 (MQ=255) gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa < 1:73839/71‑1 (MQ=255) | GTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA > minE/311389‑311525 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |