Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,234	G→A	11.5%	H285H (CAC→CAT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,234	0	G	A	11.5%	164.5 / 10.4	61	H285H (CAC→CAT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base A (4/3);  total (33/28)
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.34e-01

ACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGA  >  minE/313166‑313302                                                                     |                                                                     aCTAAAATAGCGTTTGGATTAAGCCCTTCTTCAAGACGTGCGCGCACGAACATACCCGGCAGCAGAGtgtg                                                                    <  1:2866310/71‑1 (MQ=255)  cTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGa                                                                   >  1:2653616/1‑71 (MQ=255)     aaaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                >  1:2839355/1‑71 (MQ=255)     aaaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                >  1:1643598/1‑71 (MQ=255)      aaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                <  1:3018348/70‑1 (MQ=255)           cGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTc                                                          <  1:1249372/71‑1 (MQ=255)             tttGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCgg                                                        >  1:2870226/1‑71 (MQ=255)               tGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGtgtg                                                                    >  1:2929418/1‑57 (MQ=255)               tGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCggg                                                       <  1:3586547/70‑1 (MQ=255)                ggATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGa                                                      <  1:1256628/70‑1 (MQ=255)                 gATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGaa                                                     <  1:1484255/70‑1 (MQ=255)                  aTTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCgg                                                        >  1:782222/1‑66 (MQ=255)                  aTTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                   <  1:2252637/71‑1 (MQ=255)                   ttAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                   <  1:3160175/70‑1 (MQ=255)                       gCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAgcg                                              >  1:1054775/1‑71 (MQ=255)                         ccTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAg                                                >  1:2731923/1‑67 (MQ=255)                         ccTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                            >  1:2621553/1‑71 (MQ=255)                          cttcttCCAGACGTGCGCGCACGAACATACCCGGAAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                            >  1:2840704/1‑70 (MQ=255)                             cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                >  1:2993819/1‑47 (MQ=255)                             cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAg                                         <  1:3262776/70‑1 (MQ=255)                              ttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                       >  1:38203/1‑71 (MQ=255)                               tcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                      >  1:1407321/1‑71 (MQ=255)                                ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGaa                                                     <  1:1627026/55‑1 (MQ=255)                                ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                       <  1:1014023/69‑1 (MQ=255)                                ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                      <  1:1319217/70‑1 (MQ=255)                                 cAGACGTGCGCGCACGAACATACCCGGCAGCAGTGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTg                                     <  1:177122/70‑1 (MQ=255)                                 cAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTg                                     <  1:1392199/70‑1 (MQ=255)                                  atccGTGCACGCACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTAt                                   >  1:1440171/4‑71 (MQ=255)                                     cGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGa                                >  1:2506084/1‑71 (MQ=255)                                     cGTGCACGCACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGa                                <  1:2558874/71‑1 (MQ=255)                                      gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                               >  1:963958/1‑71 (MQ=255)                                      gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                               >  1:2161671/1‑71 (MQ=255)                                      gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                               <  1:2972675/71‑1 (MQ=255)                                         cgcgcACGAACATACCCGGCAGCAGAGTGTGATcc                                                                <  1:3899538/35‑1 (MQ=255)                                         cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                             <  1:342478/70‑1 (MQ=255)                                         cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                             <  1:347168/70‑1 (MQ=255)                                         cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa                            <  1:968922/71‑1 (MQ=255)                                          gcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAg                           <  1:2866077/71‑1 (MQ=255)                                            gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                         >  1:2520459/1‑71 (MQ=255)                                            gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                         >  1:2547432/1‑71 (MQ=255)                                            gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                         >  1:2524937/1‑71 (MQ=255)                                             cACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTgg                        <  1:2958123/71‑1 (MQ=255)                                              aCGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGt                       <  1:1805484/71‑1 (MQ=255)                                                gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCt                     >  1:642806/1‑71 (MQ=255)                                                aaaCATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGag                               <  1:1671618/60‑2 (MQ=255)                                                 aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGa                                    >  1:378178/1‑55 (MQ=255)                                                 aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTg                    <  1:3167200/71‑1 (MQ=255)                                                 aaCATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGcacg                                            <  1:1086704/47‑1 (MQ=255)                                                  aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                   <  1:1570178/39‑1 (MQ=255)                                                  aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                               >  1:3023332/1‑59 (MQ=255)                                                  aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                   >  1:1242164/1‑71 (MQ=255)                                                  aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                   >  1:602593/1‑71 (MQ=255)                                                    aTACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATa                                                 >  1:2628368/1‑39 (MQ=255)                                                     tACCCGTCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGt                                           >  1:3479220/1‑44 (MQ=255)                                                      aCCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATa                                                 >  1:1083030/1‑37 (MQ=255)                                                          ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACgg            >  1:1445506/1‑70 (MQ=255)                                                          ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACg             >  1:3948476/1‑69 (MQ=255)                                                          ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGt           <  1:688501/71‑1 (MQ=255)                                                               cagAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt      >  1:3223652/1‑71 (MQ=255)                                                               cagAGTGTGATCCGGGTTCGGGAAGATAGAGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt      >  1:12574/1‑71 (MQ=255)                                                                 gagTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAAc        >  1:954253/1‑67 (MQ=255)                                                                  agTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt      >  1:1175292/1‑68 (MQ=255)                                                                  agTGTGATCCAGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCag   <  1:2877761/71‑1 (MQ=255)                                                                   gtgtGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATc                 <  1:1593406/56‑1 (MQ=255)                                                                     gtgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCaga  >  1:1108409/1‑69 (MQ=255)                                                                     |                                                                     ACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGA  >  minE/313166‑313302

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