Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,778	G→T	7.7%	T104K (ACA→AAA) ‡	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,778	0	G	T	7.7%	302.6 / 6.6	90	T104K (ACA→AAA) ‡	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base T (4/3);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

ATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313708‑313847                                                                       |                                                                      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)                       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:430368/60‑1 (MQ=255)                                   tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                <  1:999900/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:3178778/1‑60 (MQ=255)                                     ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                              >  1:3545509/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:2229540/69‑1 (MQ=255)                                          tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                <  1:316824/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:3671282/44‑1 (MQ=255)                                                        ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:1586394/44‑1 (MQ=255)                                                        ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                           <  1:2144366/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)                                                                       |                                                                      ATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGA  >  minE/313708‑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: ‑