Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,590,722	(T)6→5	coding (323/375 nt)	yfeK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,590,717	0	T	.	100.0%	363.3 / NA	82	coding (318/375 nt)	yfeK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (47/35);  total (47/35)

GTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTT  >  minE/1590651‑1590781                                                                   |                                                                 gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAATAGCGATGAGAACGCACAGCCtttt                                                               <  1:2749705/70‑1 (MQ=255) gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttt                                                               <  1:3349246/70‑1 (MQ=255) gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttt                                                               <  1:9007/70‑1 (MQ=255) gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttt                                                               <  1:4706428/70‑1 (MQ=255) gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttt                                                               <  1:4064068/70‑1 (MQ=255) gtcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttt                                                               <  1:2373007/70‑1 (MQ=255)  tcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:2110483/71‑2 (MQ=255)  tcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:840338/71‑2 (MQ=255)  tcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:2087505/71‑2 (MQ=255)  tcATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:1955326/71‑2 (MQ=255)   cATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:1372625/70‑2 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttttaca                                                           >  1:106006/1‑68 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttttaca                                                           >  1:1984375/1‑68 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttttaca                                                           >  1:3656284/1‑68 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttttaca                                                           >  1:3645720/1‑68 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCtttttaca                                                           >  1:1759357/1‑68 (MQ=255)    aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCttttta                                                             <  1:4113604/69‑2 (MQ=255)       gATTACTGGGAAGCCGTGTATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:1758376/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATg                                                        >  1:1719466/1‑70 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:107682/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:1745410/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:699930/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:471404/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:56592/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:1875737/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:4567230/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:3666685/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:661399/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:4559445/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:3523993/1‑71 (MQ=255)       gATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGc                                                       >  1:4556424/1‑71 (MQ=255)           aCTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:461719/70‑1 (MQ=255)           aCTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:2469171/70‑1 (MQ=255)           aCTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:2564880/70‑1 (MQ=255)           aCTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:1462241/70‑1 (MQ=255)            cTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:2357899/69‑1 (MQ=255)                gAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGtt                                                    <  1:4624129/65‑1 (MQ=255)                 aaGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTgcg                                             >  1:252953/1‑71 (MQ=255)                 aaGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TATTTACATGCGTTAATTgcg                                             >  1:1583654/1‑71 (MQ=255)                        atatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACg                                       >  1:2478003/1‑70 (MQ=255)                        atatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACg                                       >  1:1285512/1‑70 (MQ=255)                        atatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACg                                       >  1:23214/1‑70 (MQ=255)                        atatTGTGAAGATCCCAGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACg                                       >  1:1080561/1‑70 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:3196860/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:4803953/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:4126462/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:1160835/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:2735334/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:1071333/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:1368937/65‑1 (MQ=255)                         tatTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:944227/65‑1 (MQ=255)                             gtgAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATaa                                  <  1:4513909/70‑1 (MQ=255)                             gtgAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATaa                                  <  1:3274345/70‑1 (MQ=255)                             gtgAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATaa                                  <  1:760152/70‑1 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCAGAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:2758244/1‑71 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:4910241/1‑71 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:1346117/1‑71 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:172594/1‑71 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:2162204/1‑71 (MQ=255)                                 aGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACgg                             >  1:512866/1‑71 (MQ=255)                                  gaaccccGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:806235/68‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:4217139/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:1624209/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:3775419/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:1335720/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:1341787/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:3230046/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:1856806/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:2017391/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:2462901/71‑1 (MQ=255)                                  gATCCCCGGTAAGAGAGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:3017818/71‑1 (MQ=255)                                  gATCCCCGGTAAGAACGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            <  1:3158403/71‑1 (MQ=255)                                      cccGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCa                                           <  1:3668460/52‑1 (MQ=255)                                        cGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGc                      <  1:3877873/71‑1 (MQ=255)                                          gTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCg                     >  1:3507166/1‑70 (MQ=255)                                          gTAAGAGCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCg                     >  1:1423270/1‑70 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTgcgc                                            >  1:970224/1‑41 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCggaagga               >  1:1088461/1‑70 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCggaagga               >  1:3060220/1‑70 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCggaagga               >  1:4731186/1‑70 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCggaagga               >  1:1755883/1‑70 (MQ=255)                                                gCGATGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCggaagga               >  1:2436534/1‑70 (MQ=255)                                                    tGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTa          <  1:1286082/71‑1 (MQ=255)                                                    tGAGAACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTa          <  1:565874/71‑1 (MQ=255)                                                      agaACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGTATAAAACGGTGCCTGCGGAAGGAAATTAAt        >  1:395993/1‑71 (MQ=255)                                                      agaACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGTATAAAACGGTGCCTGCGGAAGGAAATTAAt        >  1:867208/1‑71 (MQ=255)                                                      agaACGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGTATAAAACGGTGCCTGCGGAAGGAAATTAAt        >  1:1699477/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGt                            >  1:2770239/1‑47 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGc                          >  1:3900020/1‑49 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAAt        >  1:4204737/1‑67 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:2462260/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:2407069/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:3017380/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:3107237/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:3591094/1‑71 (MQ=255)                                                          cGCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc    >  1:936837/1‑71 (MQ=255)                                                           gCACAGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCt   >  1:1191350/1‑71 (MQ=255)                                                            cacaGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCtt  <  1:4692880/71‑1 (MQ=255)                                                            cacaGCC‑TTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCtt  <  1:935786/71‑1 (MQ=255)                                                                   |                                                                 GTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTT  >  minE/1590651‑1590781

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