Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,445 T→C 8.7% S476S (TCA→TCG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4450TC 8.7% 258.0 / 13.9 92S476S (TCA→TCGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (37/47);  new base C (3/5);  total (40/52)
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.82e-01

GGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACG  >  minE/311377‑311514
                                                                    |                                                                     
ggTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTg                                                                      >  1:876342/1‑70 (MQ=255)
ggTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTg                                                                      >  1:1291885/1‑70 (MQ=255)
    tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  >  1:3133498/1‑70 (MQ=255)
    tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  <  1:4109940/70‑1 (MQ=255)
    tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  >  1:2111685/1‑70 (MQ=255)
    tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  <  1:1858629/70‑1 (MQ=255)
    tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  <  1:1169174/70‑1 (MQ=255)
      aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                                               >  1:268024/1‑71 (MQ=255)
      aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                                               <  1:1121023/71‑1 (MQ=255)
      aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                                               >  1:1570549/1‑71 (MQ=255)
         aTGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATgg                                                             >  1:792499/1‑70 (MQ=255)
         aTGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                                            <  1:3270665/71‑1 (MQ=255)
            gtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                                          <  1:3944650/70‑1 (MQ=255)
               gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                                               >  1:1696424/1‑62 (MQ=255)
               gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        >  1:1001928/1‑69 (MQ=255)
               gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        >  1:2363516/1‑69 (MQ=255)
                cacaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                                          <  1:2640482/66‑1 (MQ=255)
                 acaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                                     <  1:1246194/54‑1 (MQ=255)
                 acaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                                     <  1:1627369/54‑1 (MQ=255)
                   aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                                     <  1:3044367/52‑1 (MQ=255)
                   aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                                     <  1:3397926/52‑1 (MQ=255)
                   aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  <  1:958860/55‑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)
                             aGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                                               <  1:10949/48‑1 (MQ=255)
                                 aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                     >  1:362022/1‑70 (MQ=255)
                                 aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                     >  1:3684563/1‑70 (MQ=255)
                                   gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        <  1:2757028/49‑1 (MQ=255)
                                     tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                                            <  1:543723/43‑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:3383940/1‑48 (MQ=255)
                                                                  gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                          >  1:2858793/1‑48 (MQ=255)
                                                                  gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:884535/1‑71 (MQ=255)
                                                                  gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:3593586/1‑71 (MQ=255)
                                                                  gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:3313949/1‑71 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:1075745/1‑70 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:1052156/1‑70 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg  <  1:3604915/71‑1 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg  <  1:445949/71‑1 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg  <  1:2696760/71‑1 (MQ=255)
                                                                   cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg  <  1:2634234/71‑1 (MQ=255)
                                                                    |                                                                     
GGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACG  >  minE/311377‑311514

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: