Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,931,105	G→A	19.4%	L115L (CTG→TTG)	yjgR ←	predicted ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,931,105	0	G	A	19.4%	144.0 / 27.1	67	L115L (CTG→TTG)	yjgR	predicted ATPase
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base A (6/7);  total (31/36)
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.31e-01

GATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGG  >  minE/2931037‑2931173                                                                     |                                                                     gatgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGAt                                                                    <  1:3518850/71‑1 (MQ=255)   tgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGa                                                                     <  1:2951095/68‑1 (MQ=255)   tgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATcc                                                                  <  1:3366484/71‑1 (MQ=255)   tgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATc                                                                   >  1:2278701/1‑70 (MQ=255)     atatTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCg                                                                 >  1:3344111/1‑70 (MQ=255)      tatTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGaa                                                               <  1:3304566/71‑1 (MQ=255)      tatTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGGCGTGCCAGCAACAGCGGCCCCAGATCCGaa                                                               <  1:1755530/71‑1 (MQ=255)       atTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGaaa                                                              >  1:1770613/1‑71 (MQ=255)       atTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGaa                                                               <  1:3551583/70‑1 (MQ=255)        ttCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAg                                                                      <  1:1554784/62‑1 (MQ=255)          cAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAAc                                                             <  1:3019472/69‑1 (MQ=255)           aGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATACGAAACCGt                                                          >  1:805883/1‑71 (MQ=255)             cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATc                                                                   >  1:2635046/1‑60 (MQ=255)             cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATc                                                                   >  1:600843/1‑60 (MQ=255)             cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATc                                                                   >  1:2000186/1‑60 (MQ=255)             cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCg                                                        >  1:3379484/1‑71 (MQ=255)                gCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGccc                                                     <  1:3491927/71‑1 (MQ=255)                   aGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCa                                                  <  1:550458/71‑1 (MQ=255)                   aGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACAGTCg                                                        <  1:3191248/65‑1 (MQ=255)                    gATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCAc                                                 <  1:3717088/71‑1 (MQ=255)                     aTTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAAcc                                                            <  1:87660/59‑1 (MQ=255)                       tGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCg                                               >  1:2707966/1‑70 (MQ=255)                       tGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCg                                               >  1:2747651/1‑70 (MQ=255)                        gCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCgg                                              <  1:2562077/70‑1 (MQ=255)                        gCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGa                                             >  1:1625304/1‑71 (MQ=255)                         cacaTCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGa                                             <  1:3189984/70‑1 (MQ=255)                             tCGTTGAGATTCAACAGTCGGGCCAGCAACAGCGGCCCCAGATCCGAAAcc                                                            <  1:765645/51‑1 (MQ=255)                               gTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTc                                                         <  1:773134/52‑1 (MQ=255)                               gTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCtt                                      >  1:766447/1‑71 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGa                                                                     >  1:3216027/1‑37 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtct                                   >  1:3773398/1‑71 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtct                                   >  1:2204403/1‑71 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtct                                   >  1:2405791/1‑71 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                    <  1:2800531/70‑1 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                    <  1:290944/70‑1 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                    <  1:648991/70‑1 (MQ=255)                                  gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                    <  1:13691/70‑1 (MQ=255)                                   agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtctc                                  >  1:1467180/1‑71 (MQ=255)                                   agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtct                                   <  1:566700/70‑1 (MQ=255)                                   agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtct                                   <  1:38610/70‑1 (MQ=255)                                    gaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAAcc                                                            >  1:977412/1‑44 (MQ=255)                                     aTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGc                                                       <  1:3222500/48‑1 (MQ=255)                                      ttCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGAt                                            <  1:2590919/58‑1 (MQ=255)                                        cAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCa                              <  1:1884849/70‑1 (MQ=255)                                           cAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCa                              >  1:3450867/1‑67 (MQ=255)                                             gTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGata                        <  1:3154488/71‑1 (MQ=255)                                             gTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGata                        <  1:1555878/71‑1 (MQ=255)                                                 tGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATgcc                    <  1:2679857/71‑1 (MQ=255)                                                    cAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCcaca                 >  1:3793672/1‑71 (MQ=255)                                                    cAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCcaca                 >  1:3867798/1‑71 (MQ=255)                                                      gCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATcc                     <  1:2679261/65‑1 (MQ=255)                                                      gCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCacac                >  1:1749236/1‑70 (MQ=255)                                                        aaCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCAcaccac             <  1:2896052/71‑1 (MQ=255)                                                        aaCAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATg                                           <  1:1414572/41‑1 (MQ=255)                                                         aCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGa                                             <  1:1885889/38‑1 (MQ=255)                                                          cAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCAcaccacca           >  1:1229484/1‑71 (MQ=255)                                                            gCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACaccaccacc         >  1:2043381/1‑71 (MQ=255)                                                            gCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAg                           <  1:1673881/53‑1 (MQ=255)                                                             cGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtctc                                  >  1:1044703/1‑45 (MQ=255)                                                               gCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGa      <  1:783302/71‑1 (MQ=255)                                                                ccccAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                    <  1:888465/40‑1 (MQ=255)                                                                ccccAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGAt     <  1:2327835/71‑1 (MQ=255)                                                                  ccAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTg   >  1:1015029/1‑71 (MQ=255)                                                                  ccAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTg   >  1:3176546/1‑71 (MQ=255)                                                                  ccAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTg   >  1:2632900/1‑71 (MQ=255)                                                                  ccAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTg   >  1:3556035/1‑71 (MQ=255)                                                                  ccAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTg   >  1:3345262/1‑71 (MQ=255)                                                                    aGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGGTATCCCACACCACCACCGGATTgg  >  1:269848/1‑70 (MQ=255)                                                                     |                                                                     GATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGG  >  minE/2931037‑2931173

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