Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,910,326	(A)6→7	coding (1318/2697 nt)	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,910,320	1	.	A	95.8%	145.9 / 0.8	48	coding (1312/2697 nt)	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (25/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

ACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑AAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGA  >  minE/2910252‑2910381                                                                      |                                                              aCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑aa                                                             <  1:571805/71‑1 (MQ=255) aCCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑a                                                              >  1:375135/1‑70 (MQ=255)   cTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑aaaa                                                           <  1:2050670/71‑1 (MQ=255)           tACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATtgat                                                    <  1:1890619/71‑1 (MQ=255)               gCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATgaga                                                <  1:329431/71‑1 (MQ=255)                       cGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGttt                                        >  1:951650/1‑71 (MQ=255)                       cGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑AAAAAATTGATgag                                                 <  1:711640/61‑1 (MQ=255)                        gAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑aaaaaaa                                                        >  1:757173/1‑52 (MQ=255)                         aaGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTg                                       <  1:1526244/70‑1 (MQ=255)                          aGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGAt                                     <  1:1222783/71‑1 (MQ=255)                             tACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTc                                  <  1:2027701/71‑1 (MQ=255)                              aCGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGAtt                                              >  1:85685/1‑58 (MQ=255)                              aCGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTc                                  >  1:502479/1‑70 (MQ=255)                               cGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTc                                  >  1:1941016/1‑69 (MQ=255)                                ggATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGa                                >  1:1716900/1‑70 (MQ=255)                                 gATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGc                              >  1:2132949/1‑71 (MQ=255)                                     aaGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCg                          >  1:1969649/1‑71 (MQ=255)                                     aaGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCg                          >  1:1089646/1‑71 (MQ=255)                                      agagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGc                         >  1:191199/1‑71 (MQ=255)                                      agagTCAGCGCGCTCGCTGGCCAGTC‑TTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGcc                        >  1:1320515/1‑71 (MQ=255)                                        agTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGAt                                     >  1:1287417/1‑57 (MQ=255)                                        agTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑aaaaa                                                          <  1:1371360/35‑1 (MQ=255)                                         gTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCg                       >  1:776275/1‑70 (MQ=255)                                          tCAGCGCGCTCGCTGGCCAGTCGTTGGC‑aaaaaaa                                                        <  1:1500272/35‑2 (MQ=255)                                            aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATg                   >  1:2358261/1‑71 (MQ=255)                                               gcgcTCGCTGGCCAGTCGTTGGCAAAAAAATtgatga                                                  <  1:1738458/37‑1 (MQ=38)                                               gcgcTCGCTGGCCAGTCGTTGGCAAAAAAATtgatga                                                  <  1:893390/37‑1 (MQ=38)                                                cgctcgctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGAtt                                    <  1:2008695/50‑1 (MQ=255)                                                  ctcgctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTc                                  >  1:393288/1‑50 (MQ=255)                                                   tcgctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtgg            >  1:1399870/1‑71 (MQ=255)                                                     gctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg             >  1:2006556/1‑68 (MQ=255)                                                     gctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCggtg             >  1:1733195/1‑68 (MQ=255)                                                      ctGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAg         <  1:240052/71‑1 (MQ=255)                                                       tGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCg                       <  1:2348237/56‑1 (MQ=255)                                                       tGGCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTa          >  1:407351/1‑69 (MQ=255)                                                         gCCAGTCGTTGGCAAAAAAATTGCTGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt        <  1:75990/69‑1 (MQ=255)                                                         gCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt        <  1:2221510/69‑1 (MQ=255)                                                         gCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt        <  1:1503383/69‑1 (MQ=255)                                                         gCCAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt        <  1:1374856/69‑1 (MQ=255)                                                          ccAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     <  1:2006125/71‑1 (MQ=255)                                                          ccAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     <  1:510057/71‑1 (MQ=255)                                                          ccAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     <  1:932274/71‑1 (MQ=255)                                                          ccAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     <  1:579879/71‑1 (MQ=255)                                                           cAGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTgg      >  1:347168/1‑69 (MQ=255)                                                            aGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     >  1:114966/1‑69 (MQ=255)                                                            aGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGc     >  1:470125/1‑69 (MQ=255)                                                            aGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAg   <  1:329985/71‑1 (MQ=255)                                                            aGTCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAg   <  1:237033/71‑1 (MQ=255)                                                            aGTCGTTGGC‑AAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTc                 >  1:1218205/1‑56 (MQ=255)                                                              tCGTTGGCAAAAAAATTGATGAGATTCCGTTTGAt                                     <  1:2204041/35‑1 (MQ=38)                                                              tCGTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGa  >  1:1191229/1‑70 (MQ=255)                                                                gTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCgg                      >  1:947052/1‑48 (MQ=255)                                                                gTTGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGa                     >  1:427905/1‑49 (MQ=255)                                                                 ttGGCAAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATg                   >  1:1300745/1‑50 (MQ=255)                                                                      |                                                              ACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC‑AAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGA  >  minE/2910252‑2910381

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