Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,795 A→G 6.8% P98P (CCT→CCC)  acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,7950AG 6.8% 253.5 / 8.4 88P98P (CCT→CCC) acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base G (3/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.42e-01

GATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCA  >  minE/313726‑313859
                                                                     |                                                                
gATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                                                 >  1:309744/1‑71 (MQ=255)
 atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                <  1:2961427/71‑1 (MQ=255)
   atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                               >  1:4093586/1‑70 (MQ=255)
   atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                              >  1:1425047/1‑71 (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)
           cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                      <  1:2514709/71‑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:1907769/66‑1 (MQ=255)
                     tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                 <  1:3297119/66‑1 (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:1184254/70‑1 (MQ=255)
                       tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                           <  1:3670336/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:2144366/44‑1 (MQ=255)
                                     ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                                       <  1:1586394/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:2159283/36‑1 (MQ=255)
                                         cACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                           <  1:2013728/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:1259639/71‑1 (MQ=255)
                                                   tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa              <  1:2436995/71‑1 (MQ=255)
                                                     tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                <  1:434208/67‑1 (MQ=255)
                                                      cGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                             >  1:4047242/1‑53 (MQ=255)
                                                       gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGc                          >  1:1688524/1‑55 (MQ=255)
                                                        ccTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGa                                >  1:2331952/1‑48 (MQ=255)
                                                          tGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                           <  1:2989797/35‑1 (MQ=255)
                                                           gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc      >  1:86362/1‑71 (MQ=255)
                                                            aTAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa              <  1:2817295/62‑1 (MQ=255)
                                                            aTAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTa         >  1:2662533/1‑67 (MQ=255)
                                                             tAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCt     >  1:3936640/1‑70 (MQ=255)
                                                               ggTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTc   <  1:225172/70‑1 (MQ=255)
                                                               ggTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCa  <  1:3618414/71‑1 (MQ=255)
                                                                   gCAGGATCAATCTGATAGAGAGAGACACCTTCTTCGAt                               <  1:1585523/38‑1 (MQ=255)
                                                                   gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGAt                               <  1:2488304/38‑1 (MQ=255)
                                                                     |                                                                
GATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCA  >  minE/313726‑313859

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: ‑