Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,778	G→T	6.2%	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	6.2%	266.5 / 2.5	80	T104K (ACA→AAA) ‡	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base T (3/2);  total (42/38)
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.74e-01

TCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTG  >  minE/313710‑313846                                                                     |                                                                     tCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTc                                                                    >  1:185794/1‑71 (MQ=255)  cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   >  1:3297553/1‑71 (MQ=255)  cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   >  1:3302138/1‑71 (MQ=255)  cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   >  1:1357740/1‑71 (MQ=255)   aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGc                                                                  <  1:355407/71‑1 (MQ=255)    ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   <  1:1414649/69‑1 (MQ=255)    ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   <  1:1466063/69‑1 (MQ=255)     cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                                                >  1:2733143/1‑71 (MQ=255)      gTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGAACACCTTTCGCACTGTCGTATGTCGCCTg                                                               <  1:3397587/71‑1 (MQ=255)       tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                              <  1:842687/71‑1 (MQ=255)       tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                              <  1:2672726/71‑1 (MQ=255)         aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGTCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                             <  1:786936/70‑1 (MQ=255)         aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                              >  1:3890704/1‑69 (MQ=255)         aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                             <  1:588666/70‑1 (MQ=255)         aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                             <  1:2664489/70‑1 (MQ=255)         aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                             <  1:2552917/70‑1 (MQ=255)           ttGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAgg                                                          >  1:383311/1‑71 (MQ=255)            tGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAg                                                           >  1:2139751/1‑69 (MQ=255)                   tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                   <  1:2881273/70‑1 (MQ=255)                   tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                  >  1:3806933/1‑71 (MQ=255)                     ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                   <  1:3325655/52‑1 (MQ=255)                     ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                <  1:2134782/71‑1 (MQ=255)                       ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTc                                                                    >  1:699980/1‑49 (MQ=255)                       ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGt                                                         <  1:2684919/60‑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:3621205/70‑1 (MQ=255)                         cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                             <  1:2781175/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)                                     ttttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGt                                                         <  1:2940451/45‑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)                                         cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                                        >  1:137180/1‑43 (MQ=255)                                         cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                                        >  1:2077633/1‑43 (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:3783257/1‑41 (MQ=255)                                                      ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                             >  1:920251/1‑41 (MQ=255)                                                      ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                             >  1:2754584/1‑41 (MQ=255)                                                      ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                             >  1:2717618/1‑41 (MQ=255)                                                      ttCGCAATGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                     >  1:3799638/1‑65 (MQ=255)                                                       tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                             >  1:3475039/1‑40 (MQ=255)                                                       tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                             >  1:2542452/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)                                                                     |                                                                     TCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTG  >  minE/313710‑313846

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