Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,774	G→A	7.8%	Y105Y (TAC→TAT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,774	0	G	A	7.8%	281.8 / 10.2	89	Y105Y (TAC→TAT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (38/44);  new base A (4/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTT  >  minE/313704‑313844                                                                       |                                                                       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:3725478/1‑71 (MQ=255)    gATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                                     <  1:3816526/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)         aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                                      >  1:1301805/1‑65 (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:2953146/70‑1 (MQ=255)                             ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                             <  1:3812138/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: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: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: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)                                                                       |                                                                       AACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTT  >  minE/313704‑313844

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