Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,768	A→G	6.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	6.1%	233.4 / 4.8	82	S107S (AGT→AGC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base G (3/2);  total (43/39)
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.94e-01

TCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGT  >  minE/313698‑313832                                                                       |                                                                 tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                  >  1:906018/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                  >  1:2098231/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                  >  1:1958271/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                  >  1:2783148/1‑71 (MQ=255) tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                                  >  1:3100749/1‑71 (MQ=255)     aTAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTg                                                               <  1:2406632/70‑1 (MQ=255)     aTAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGt                                                              <  1:3948668/71‑1 (MQ=255)      tAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTc                                                             >  1:1265410/1‑71 (MQ=255)        aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            >  1:488898/1‑70 (MQ=255)        aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            >  1:3154363/1‑70 (MQ=255)        aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            <  1:750416/70‑1 (MQ=255)        aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGt                                                           <  1:2857252/71‑1 (MQ=255)             tCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTc                                                      >  1:185794/1‑71 (MQ=255)              cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     >  1:1357740/1‑71 (MQ=255)              cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     >  1:3302138/1‑71 (MQ=255)              cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     >  1:3297553/1‑71 (MQ=255)              cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGC‑‑GATCACCTTTCGCa                                                                  >  1:2478015/1‑56 (MQ=38)               aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGc                                                    <  1:355407/71‑1 (MQ=255)                ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     <  1:1466063/69‑1 (MQ=255)                ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     <  1:1414649/69‑1 (MQ=255)                 cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                                  >  1:2733143/1‑71 (MQ=255)                  gTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGAACACCTTTCGCACTGTCGTATGTCGCCTg                                                 <  1:3397587/71‑1 (MQ=255)                   tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                <  1:2672726/71‑1 (MQ=255)                   tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                <  1:842687/71‑1 (MQ=255)                     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGTCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                               <  1:786936/70‑1 (MQ=255)                     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                >  1:3890704/1‑69 (MQ=255)                     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                               <  1:2664489/70‑1 (MQ=255)                     aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                               <  1:588666/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)                             gATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCAc                                                                 >  1:272578/1‑44 (MQ=255)                               tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                     <  1:2881273/70‑1 (MQ=255)                               tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                    >  1:3806933/1‑71 (MQ=255)                                atTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                         >  1:1152723/1‑49 (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: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)                                             tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                          >  1:2001585/1‑35 (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:3551790/70‑1 (MQ=255)                                                    tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                <  1:1487958/70‑1 (MQ=255)                                                     cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                          >  1:2077633/1‑43 (MQ=255)                                                     cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                          >  1:137180/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:3695322/1‑71 (MQ=255)                                                                 tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt  >  1:3272140/1‑71 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:2717618/1‑41 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:2754584/1‑41 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:3783257/1‑41 (MQ=255)                                                                  ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:920251/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)                                                                       |                                                                 TCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGT  >  minE/313698‑313832

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