Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,258	G→A	9.2%	R277R (CGC→CGT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,258	0	G	A	9.2%	188.9 / 6.8	65	R277R (CGC→CGT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base A (4/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGA  >  minE/313188‑313328                                                                       |                                                                       gCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAgcg                                                                        >  1:1054775/1‑71 (MQ=255)   ccTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                                                      >  1:2621553/1‑71 (MQ=255)    cttcttCCAGACGTGCGCGCACGAACATACCCGGAAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                                                      >  1:2840704/1‑70 (MQ=255)       cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAg                                                                   <  1:3262776/70‑1 (MQ=255)        ttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                                                 >  1:38203/1‑71 (MQ=255)         tcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                                                >  1:1407321/1‑71 (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)                   cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                                                       <  1:347168/70‑1 (MQ=255)                   cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                                                       <  1:342478/70‑1 (MQ=255)                   cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa                                                      <  1:968922/71‑1 (MQ=255)                    gcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAg                                                     <  1:2866077/71‑1 (MQ=255)                      gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                                   >  1:2547432/1‑71 (MQ=255)                      gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                                   >  1:2524937/1‑71 (MQ=255)                      gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                                   >  1:2520459/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)                            aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                         >  1:3023332/1‑59 (MQ=255)                            aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                             >  1:602593/1‑71 (MQ=255)                            aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                             >  1:1242164/1‑71 (MQ=255)                               tACCCGTCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGt                                                                     >  1:3479220/1‑44 (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)                                                tgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGa                                                          >  1:1712043/1‑38 (MQ=255)                                                tgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAcc                                                        <  1:2816754/40‑1 (MQ=255)                                                 gATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAt                        >  1:1660834/1‑71 (MQ=255)                                                  aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCa                                          <  1:3697904/52‑1 (MQ=255)                                                  aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCa                                          <  1:3886428/52‑1 (MQ=255)                                                  aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAtt                       >  1:2336518/1‑71 (MQ=255)                                                  aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAt                        <  1:1027793/70‑1 (MQ=255)                                                    ccGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTc                      >  1:968542/1‑70 (MQ=255)                                                          tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGAt                                            >  1:3545643/1‑42 (MQ=255)                                                          tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGAt                                            >  1:3796927/1‑42 (MQ=255)                                                          tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTa               >  1:3653820/1‑71 (MQ=255)                                                           cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTa               <  1:394430/70‑1 (MQ=255)                                                           cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc              >  1:2428621/1‑71 (MQ=255)                                                           cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc              >  1:1996961/1‑71 (MQ=255)                                                             ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTg                                              <  1:2761731/37‑1 (MQ=255)                                                             ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCg            <  1:724886/71‑1 (MQ=255)                                                              gAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGt           >  1:2879162/1‑71 (MQ=255)                                                                aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc              >  1:2847490/1‑66 (MQ=255)                                                                aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTcc         >  1:123678/1‑71 (MQ=255)                                                                aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTc          >  1:2122436/1‑70 (MQ=255)                                                                 gATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCt        <  1:1398938/71‑1 (MQ=255)                                                                  aTAGCGCGTAGGGTGATAGACCCAGTGGTCTGATc                                           <  1:415541/35‑1 (MQ=255)                                                                  aTAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTg       <  1:661747/71‑1 (MQ=255)                                                                    aGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAgaga                          >  1:1962884/1‑50 (MQ=255)                                                                    aGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAgaga                          >  1:1565714/1‑50 (MQ=255)                                                                     gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCgg    >  1:3242416/1‑71 (MQ=255)                                                                       gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg   >  1:2252918/1‑70 (MQ=255)                                                                       gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg   >  1:2847050/1‑70 (MQ=255)                                                                       gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGa  <  1:382082/71‑1 (MQ=255)                                                                       |                                                                       GCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGA  >  minE/313188‑313328

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