Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,780 C→T 7.9% A103A (GCG→GCA)  acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,7800CT 7.9% 276.7 / 10.5 88A103A (GCG→GCA) acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base T (4/3);  total (44/44)
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

CCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313713‑313847
                                                                   |                                                                   
ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGcc                                                                  >  1:93367/1‑71 (MQ=255)
 cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                                                 <  1:3959604/71‑1 (MQ=255)
   tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                               <  1:2368864/71‑1 (MQ=255)
     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                               >  1:908203/1‑69 (MQ=255)
     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATa                                                             <  1:1722934/71‑1 (MQ=255)
       ttGCGCGATATTGGCTGCAGCCTGGGCCTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAg                                                            >  1:3445291/1‑70 (MQ=255)
        tGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGt                                                          >  1:3636597/1‑71 (MQ=255)
             gATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                                     >  1:309744/1‑71 (MQ=255)
              atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                    <  1:2961427/71‑1 (MQ=255)
                atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                                                 <  1:1913551/56‑1 (MQ=255)
                atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                   >  1:4093586/1‑70 (MQ=255)
                atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                  >  1:1425047/1‑71 (MQ=255)
                 ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                               <  1:2362564/57‑1 (MQ=255)
                 ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                  >  1:1820507/1‑70 (MQ=255)
                  tGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                 >  1:1241011/1‑70 (MQ=255)
                   ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                               >  1:3447952/1‑71 (MQ=255)
                   ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                <  1:3812138/70‑1 (MQ=255)
                   ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                <  1:2953146/70‑1 (MQ=255)
                      tGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                   >  1:2681059/1‑64 (MQ=255)
                        cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGc                                                                   >  1:869541/1‑46 (MQ=255)
                        cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                          <  1:2514709/71‑1 (MQ=255)
                          gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGcc                                                                  <  1:2279407/45‑1 (MQ=255)
                           ccTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                          <  1:3646586/68‑1 (MQ=255)
                            cTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                        >  1:2684716/1‑69 (MQ=255)
                             tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                <  1:999900/60‑1 (MQ=255)
                             tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                <  1:430368/60‑1 (MQ=255)
                             tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCATTCTGATa                                        >  1:1592507/1‑68 (MQ=255)
                             tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGAt                                         >  1:575804/1‑67 (MQ=255)
                              gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                  >  1:3755217/1‑57 (MQ=255)
                               ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                              >  1:3545509/1‑60 (MQ=255)
                               ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                              >  1:3178778/1‑60 (MQ=255)
                               ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                      >  1:1097008/1‑68 (MQ=255)
                                gCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                  <  1:2717014/71‑1 (MQ=255)
                                  tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                    <  1:2007923/51‑1 (MQ=255)
                                  tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                     <  1:3297119/66‑1 (MQ=255)
                                  tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                     <  1:1907769/66‑1 (MQ=255)
                                   ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGcgg                                                     >  1:2355189/1‑47 (MQ=255)
                                   ttCGCCAGATCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                               >  1:1811369/1‑68 (MQ=255)
                                   ttCGCCAGATCACCTTTCGCACTGTCGTATGTCCCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                >  1:2040572/1‑70 (MQ=255)
                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                <  1:316824/69‑1 (MQ=255)
                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                <  1:2229540/69‑1 (MQ=255)
                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                               <  1:3670336/70‑1 (MQ=255)
                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                               <  1:1184254/70‑1 (MQ=255)
                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                               <  1:1322179/70‑1 (MQ=255)
                                      gCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                       <  1:4120862/60‑1 (MQ=255)
                                        cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                          <  1:1232969/55‑1 (MQ=255)
                                        cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                >  1:3914454/1‑65 (MQ=255)
                                        cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                           >  1:514151/1‑70 (MQ=255)
                                          gATCACCTTTCGCACTGTCGTATGTCGCCTTATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCt                        <  1:441301/71‑1 (MQ=255)
                                          gATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCt                        <  1:2203957/71‑1 (MQ=255)
                                           aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                            <  1:3090864/50‑1 (MQ=255)
                                            tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                      <  1:536999/71‑1 (MQ=255)
                                            tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                      <  1:2338972/71‑1 (MQ=255)
                                             cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                       >  1:775562/1‑53 (MQ=255)
                                             cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                      >  1:713602/1‑70 (MQ=255)
                                               ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                     <  1:2225833/69‑1 (MQ=255)
                                                cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATg                  <  1:2848394/71‑1 (MQ=255)
                                                  ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:1586394/44‑1 (MQ=255)
                                                  ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:2144366/44‑1 (MQ=255)
                                                  ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:3671282/44‑1 (MQ=255)
                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                >  1:1162506/1‑71 (MQ=255)
                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                >  1:1552036/1‑71 (MQ=255)
                                                   tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTa                                        >  1:1278226/1‑46 (MQ=255)
                                                    cgcgCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTggt                                         >  1:1986979/1‑42 (MQ=255)
                                                      cACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                               <  1:2013728/36‑1 (MQ=255)
                                                      cACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                               <  1:2159283/36‑1 (MQ=255)
                                                       aCTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                           >  1:1316272/1‑55 (MQ=255)
                                                       aCTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                           >  1:3968059/1‑55 (MQ=255)
                                                        cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                             >  1:4027660/1‑36 (MQ=255)
                                                        cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAcc          >  1:2066732/1‑71 (MQ=255)
                                                        cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAc           >  1:3722110/1‑70 (MQ=255)
                                                         tGTCGTATGTCGCCTGATAGGTCGCAGGATCTATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa            <  1:428906/68‑1 (MQ=255)
                                                         tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                     <  1:2301160/59‑1 (MQ=255)
                                                         tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa            <  1:1771386/68‑1 (MQ=255)
                                                           tCGTATGTCGCCTGATAGGTCTCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                      <  1:1791863/56‑1 (MQ=255)
                                                           tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa            >  1:3876393/1‑66 (MQ=255)
                                                           tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttc       >  1:1580177/1‑71 (MQ=255)
                                                             gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                        >  1:2030076/1‑36 (MQ=255)
                                                             gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttctt     >  1:1450486/1‑71 (MQ=255)
                                                              tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTTGCTACCTTCttt    >  1:1451079/1‑71 (MQ=255)
                                                              tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt    <  1:2620215/71‑1 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATagagagag                                 >  1:1421552/1‑40 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                    >  1:858950/1‑37 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGc              >  1:2204003/1‑59 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg   >  1:1325165/1‑70 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa  <  1:2436995/71‑1 (MQ=255)
                                                                tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa  <  1:1259639/71‑1 (MQ=255)
                                                                  tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt    <  1:434208/67‑1 (MQ=255)
                                                                   cGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                 >  1:4047242/1‑53 (MQ=255)
                                                                   |                                                                   
CCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313713‑313847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑