Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,270	A→G	7.1%	S273S (TCT→TCC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,270	0	A	G	7.1%	196.7 / 5.2	70	S273S (TCT→TCC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base G (3/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGT  >  minE/313202‑313340                                                                     |                                                                       cGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGa                                                                      >  1:2506084/1‑71 (MQ=255) cGTGCACGCACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGa                                                                      <  1:2558874/71‑1 (MQ=255)  gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                                     >  1:963958/1‑71 (MQ=255)  gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                                     <  1:2972675/71‑1 (MQ=255)  gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                                     >  1:2161671/1‑71 (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: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)             aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTg                                                          <  1:3167200/71‑1 (MQ=255)              aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                                     >  1:3023332/1‑59 (MQ=255)              aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                                         >  1:1242164/1‑71 (MQ=255)              aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                                         >  1:602593/1‑71 (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:3886428/52‑1 (MQ=255)                                    aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCa                                                      <  1:3697904/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:1996961/1‑71 (MQ=255)                                             cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc                          >  1:2428621/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:1565714/1‑50 (MQ=255)                                                      aGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAgaga                                      >  1:1962884/1‑50 (MQ=255)                                                       gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCgg                >  1:3242416/1‑71 (MQ=255)                                                         gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg               >  1:2847050/1‑70 (MQ=255)                                                         gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg               >  1:2252918/1‑70 (MQ=255)                                                         gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGa              <  1:382082/71‑1 (MQ=255)                                                          cgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAg                               <  1:2906554/53‑1 (MQ=255)                                                             aGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCa                                <  1:1542628/49‑1 (MQ=255)                                                              gtgTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTa         <  1:2061526/69‑1 (MQ=255)                                                              gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCg                 >  1:645354/1‑63 (MQ=255)                                                              gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt          <  1:541315/70‑1 (MQ=255)                                                              gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt          <  1:2527133/70‑1 (MQ=255)                                                              gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt          >  1:2741882/1‑70 (MQ=255)                                                               ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGc                              >  1:228888/1‑49 (MQ=255)                                                               ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTaa        <  1:2448983/71‑1 (MQ=255)                                                                gTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAAt       <  1:2996727/71‑1 (MQ=255)                                                                    tAGACCCAGTTGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGcc    <  1:2837161/70‑1 (MQ=255)                                                                     aGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGt  >  1:2093092/1‑71 (MQ=255)                                                                     aGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGt  >  1:920101/1‑71 (MQ=255)                                                                     |                                                                       CGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGT  >  minE/313202‑313340

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