Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,777	T→A	8.2%	T104T (ACA→ACT) ‡	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,777	0	T	A	8.2%	240.6 / 12.0	85	T104T (ACA→ACT) ‡	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base A (4/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313707‑313847                                                                       |                                                                       gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCATATCACCTTTCGCACTGTCGTAt                                                                        <  1:2998664/71‑1 (MQ=255) gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                                        >  1:3725478/1‑71 (MQ=255) gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                                        <  1:3816526/71‑1 (MQ=255) gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                                        <  1:761226/71‑1 (MQ=255)  aTTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATg                                                                       <  1:2943857/71‑1 (MQ=255)   ttCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                                      <  1:2425386/71‑1 (MQ=255)    tccccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                                      >  1:3718281/4‑70 (MQ=255)       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)             aTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                                      <  1:122118/61‑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)                        ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                                        >  1:2449599/1‑48 (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:1907769/66‑1 (MQ=255)                                         tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                     <  1:3297119/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:1322179/70‑1 (MQ=255)                                           tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                               <  1:3670336/70‑1 (MQ=255)                                           tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                               <  1:1184254/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:3671282/44‑1 (MQ=255)                                                         ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:2144366/44‑1 (MQ=255)                                                         ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:1586394/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)                                                                       |                                                                       GATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313707‑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: ‑