Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,768	A→G	7.1%	S107S (AGT→AGC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,768	0	A	G	7.1%	239.6 / 8.6	84	S107S (AGT→AGC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (32/46);  new base G (3/3);  total (35/49)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTC  >  minE/313698‑313833                                                                       |                                                                  tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:3409927/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:3108919/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:3844733/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:3306858/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:3180895/1‑71 (MQ=255)  cTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   <  1:2717168/70‑1 (MQ=255)      tAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACt                                                                 <  1:3814121/68‑1 (MQ=255)       aaCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTTGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                             >  1:1387003/1‑71 (MQ=255)        aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGt                                                            <  1:2587173/71‑1 (MQ=255)        aCGATTCAACGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                             <  1:2306680/70‑1 (MQ=255)          gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCATATCACCTTTCGCACTGTCGTAt                                                          <  1:2998664/71‑1 (MQ=255)          gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                          <  1:761226/71‑1 (MQ=255)          gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                          <  1:3816526/71‑1 (MQ=255)          gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                          >  1:3725478/1‑71 (MQ=255)           aTTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTg                                                                <  1:2562443/64‑1 (MQ=255)           aTTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATg                                                         <  1:2943857/71‑1 (MQ=255)            ttCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                        <  1:2425386/71‑1 (MQ=255)             tccccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                        >  1:3718281/4‑70 (MQ=255)               aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                           >  1:1301805/1‑65 (MQ=255)                ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                   >  1:2321227/1‑56 (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)                               tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                             <  1:1465081/47‑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)                                   ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                           >  1:2902012/1‑45 (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)                                          gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                             <  1:507384/36‑1 (MQ=255)                                          gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGcc                                                    <  1:2279407/45‑1 (MQ=255)                                           ccTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                             <  1:3489359/35‑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: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:2229540/69‑1 (MQ=255)                                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                  <  1:316824/69‑1 (MQ=255)                                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                 <  1:1322179/70‑1 (MQ=255)                                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                 <  1:1184254/70‑1 (MQ=255)                                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                 <  1:3670336/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:3671282/44‑1 (MQ=255)                                                                  ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                             <  1:2144366/44‑1 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc  >  1:1552036/1‑71 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc  >  1:1162506/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)                                                                       |                                                                  TCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTC  >  minE/313698‑313833

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