Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,795	A→G	7.6%	P98P (CCT→CCC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,795	0	A	G	7.6%	223.3 / 7.4	79	P98P (CCT→CCC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (3/3);  total (39/40)
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.98e-01

TATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  minE/313728‑313865                                                                    |                                                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                      <  1:2881273/70‑1 (MQ=255) tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                                     >  1:3806933/1‑71 (MQ=255)   ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                                   <  1:2134782/71‑1 (MQ=255)     ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                                   >  1:21546/1‑69 (MQ=255)     ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                  <  1:1700155/70‑1 (MQ=255)     ggCTGAAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                                   >  1:556766/1‑69 (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:620324/1‑69 (MQ=255)            gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                            >  1:2569564/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)                    tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                                                       <  1:3934072/50‑1 (MQ=255)                     ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggac                                                   <  1:1106380/69‑2 (MQ=255)                     ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                                    >  1:2659880/1‑52 (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:3695322/1‑71 (MQ=255)                                   tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                                   >  1:3272140/1‑71 (MQ=255)                                    ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:2717618/1‑41 (MQ=255)                                    ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:920251/1‑41 (MQ=255)                                    ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:2754584/1‑41 (MQ=255)                                    ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                >  1:3783257/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:1184875/71‑1 (MQ=255)                                                tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      >  1:2322824/1‑71 (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:3384932/69‑1 (MQ=255)                                                     tgCCTGATAGGTCGCGGGATCGATCTGGTACAGGGACTGGCCTGCTTGCACATCGCTGCCTTCAGTGaaa                  <  1:2195655/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)                                                                    |                                                                       TATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  minE/313728‑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: ‑