Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,801	A→G	8.6%	I96I (ATT→ATC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,801	0	A	G	8.6%	195.4 / 7.6	71	I96I (ATT→ATC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base G (3/3);  total (36/35)
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.97e-01

GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  minE/313732‑313865                                                                      |                                                                 ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                  <  1:1700155/70‑1 (MQ=255)   cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                <  1:2781175/70‑1 (MQ=255)   cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                <  1:3621205/70‑1 (MQ=255)   cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                                               >  1:3065537/1‑71 (MQ=255)      cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                                               <  1:3575031/68‑1 (MQ=255)      cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                              <  1:3543853/69‑1 (MQ=255)      cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            <  1:685667/71‑1 (MQ=255)       aGCCTGGGCTTTCGACAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            >  1:3426661/1‑70 (MQ=255)        gCCTGGGTTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                          <  1:657635/71‑1 (MQ=255)        gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            >  1:2569564/1‑69 (MQ=255)        gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            >  1:620324/1‑69 (MQ=255)           tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                          >  1:3057524/1‑68 (MQ=255)            gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                  <  1:458518/59‑1 (MQ=255)               cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:1896816/65‑1 (MQ=255)                 ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggac                                                   <  1:1106380/69‑2 (MQ=255)                 ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagagag                                                   >  1:2727577/1‑69 (MQ=255)                  tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        >  1:1044064/1‑63 (MQ=255)                  tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                  <  1:2579346/69‑1 (MQ=255)                  tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                 <  1:1487958/70‑1 (MQ=255)                  tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                 <  1:3551790/70‑1 (MQ=255)                   cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGaca                                                <  1:3874872/70‑1 (MQ=255)                    gCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                                               >  1:1605930/1‑54 (MQ=255)                     ccAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGAt                                                           <  1:198998/57‑1 (MQ=255)                      cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTg                                            <  1:3451331/71‑1 (MQ=255)                        gATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTg                                                             >  1:1723289/1‑52 (MQ=255)                          tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                       >  1:2331124/1‑56 (MQ=255)                          tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:1729312/71‑1 (MQ=255)                           cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:241718/1‑70 (MQ=255)                            aCCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                                       <  1:848699/70‑1 (MQ=255)                            aCCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGa                                      <  1:3254427/71‑1 (MQ=255)                             ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCtt                                         >  1:992920/1‑67 (MQ=255)                              cTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                                             >  1:911433/1‑46 (MQ=255)                               tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                                   >  1:3272140/1‑71 (MQ=255)                               tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                                   >  1:3695322/1‑71 (MQ=255)                                ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:920251/1‑41 (MQ=255)                                ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:2717618/1‑41 (MQ=255)                                ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:3783257/1‑41 (MQ=255)                                ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:2754584/1‑41 (MQ=255)                                ttCGCAATGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:3799638/1‑65 (MQ=255)                                 tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                                >  1:2542452/1‑40 (MQ=255)                                 tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                                >  1:3475039/1‑40 (MQ=255)                                       tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            >  1:552886/1‑38 (MQ=255)                                       tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                       <  1:957234/43‑1 (MQ=255)                                           gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                    >  1:1043863/1‑42 (MQ=255)                                            tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      >  1:2322824/1‑71 (MQ=255)                                            tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      <  1:1184875/71‑1 (MQ=255)                                             aTGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      >  1:2728984/1‑70 (MQ=255)                                              tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCt                               <  1:2553512/60‑1 (MQ=255)                                              tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                     >  1:867376/1‑70 (MQ=255)                                               gTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc                                           >  1:1868462/1‑47 (MQ=255)                                                tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:3076080/1‑71 (MQ=255)                                                tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  <  1:2958293/71‑1 (MQ=255)                                                 tgCCTGATAGGTCGCGGGATCGATCTGGTACAGGGACTGGCCTGCTTGCACATCGCTGCCTTCAGTGaaa                  <  1:2195655/69‑1 (MQ=255)                                                 tgCCTGATAGGTCGCGGGATCGATCTGGTACAGGGACTGGCCTGCTTGCACATCGCTGCCTTCAGTGaaa                  <  1:3384932/69‑1 (MQ=255)                                                  gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                 <  1:3660625/70‑1 (MQ=255)                                                   ccTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTa               >  1:2109699/1‑71 (MQ=255)                                                    cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc              >  1:6376/1‑71 (MQ=255)                                                     tGATAGGTCGCAGGATCTATCTGATAGAGAGAGACACCTGCTTc                                        <  1:1836590/44‑1 (MQ=255)                                                     tGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc              <  1:1833739/70‑1 (MQ=255)                                                     tGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACg             >  1:2403095/1‑71 (MQ=255)                                                      gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAc                             >  1:2280155/1‑54 (MQ=255)                                                      gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc            >  1:669171/1‑71 (MQ=255)                                                        tAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                    <  1:1060468/61‑1 (MQ=255)                                                        tAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCtt          <  1:2618231/71‑1 (MQ=255)                                                         aGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCtt          >  1:2440831/1‑70 (MQ=255)                                                           gTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                 <  1:2459011/61‑1 (MQ=255)                                                           gTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAg       <  1:1009595/71‑1 (MQ=255)                                                            tCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAtt                <  1:250702/61‑1 (MQ=255)                                                            tCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAgg      <  1:2127143/71‑1 (MQ=255)                                                            tCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATCTCGCTACCTTCTTTGAAATTACGCTTCAgg      <  1:1113183/71‑1 (MQ=255)                                                             cGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      <  1:2127964/54‑1 (MQ=255)                                                              gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGa     <  1:3420045/70‑1 (MQ=255)                                                                aGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa  >  1:1461189/1‑71 (MQ=255)                                                                      |                                                                 GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  minE/313732‑313865

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑