Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,774	G→A	6.4%	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	6.4%	240.5 / 3.8	77	Y105Y (TAC→TAT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base A (3/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGT  >  minE/313705‑313832                                                                      |                                                           aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            <  1:750416/70‑1 (MQ=255) aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            >  1:3154363/1‑70 (MQ=255) aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                            >  1:488898/1‑70 (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:3297553/1‑71 (MQ=255)       cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                     >  1:3302138/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)                         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: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: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: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:3695322/1‑71 (MQ=255)                                                          tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt  >  1:3272140/1‑71 (MQ=255)                                                           ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:920251/1‑41 (MQ=255)                                                           ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                               >  1:2717618/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: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)                                                                      |                                                           ACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGT  >  minE/313705‑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: ‑