Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,450,298	A→G	W78R (TGG→CGG)	yifE ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,450,298	0	A	G	100.0%	312.4 / NA	86	W78R (TGG→CGG)	yifE	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/42);  total (44/42)

CTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTCC  >  minE/2450230‑2450365                                                                     |                                                                    cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                   >  1:265975/1‑71 (MQ=255)  tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                  >  1:227845/1‑71 (MQ=255)  tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                  >  1:1066822/1‑71 (MQ=255)  tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                  >  1:2456400/1‑71 (MQ=255)  tGCGGTTTACAGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                   >  1:1806350/1‑70 (MQ=255)   gCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcg                                                                 >  1:1375304/1‑71 (MQ=255)     ggTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcgcg                                                               >  1:515362/1‑71 (MQ=255)      gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                              <  1:2684470/71‑1 (MQ=255)        ttACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTc                                                            >  1:1082934/1‑71 (MQ=255)        ttACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTc                                                            >  1:1727675/1‑71 (MQ=255)         tACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCt                                                           >  1:445143/1‑71 (MQ=255)          aCCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCt                                                           <  1:201482/70‑1 (MQ=255)          aCCGCCGGAAAGGGTGTGAAAACG‑TTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGtt                                                             <  1:114325/67‑1 (MQ=255)            cgccgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCt                                                           >  1:3006187/1‑68 (MQ=255)             gccgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCtt                                                       >  1:392972/1‑71 (MQ=255)             gccgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCtt                                                       <  1:1254513/71‑1 (MQ=255)              ccgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCtt                                                       <  1:2921659/70‑1 (MQ=255)               cgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCt                                                     >  1:2747142/1‑71 (MQ=255)               cgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCt                                                     <  1:643573/71‑1 (MQ=255)                 gAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCt                                                     >  1:2812488/1‑69 (MQ=255)                 gAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGt                                                   <  1:1726409/71‑1 (MQ=255)                  aaaGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTc                                                  <  1:1710290/71‑1 (MQ=255)                   aaGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCa                                                 <  1:366844/71‑1 (MQ=255)                   aaGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCa                                                 <  1:341936/71‑1 (MQ=255)                     gggTGTGAAAACGTTTTGGACGCTTAATACGCGTGATATACTTGGACCGCACGCGTTCTGCTTCTGTCACt                                               <  1:555816/71‑1 (MQ=255)                       gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCtt                                                       <  1:1674284/61‑1 (MQ=255)                       gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTgg                                             <  1:204546/71‑1 (MQ=255)                        tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                            <  1:1300387/71‑1 (MQ=255)                        tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                            >  1:333977/1‑71 (MQ=255)                        tctgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTgg                                             <  1:1192546/68‑1 (MQ=255)                          tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGCCCGCACGCGTTCTGCTTCTGTCACTGGCt                                           >  1:2946056/1‑70 (MQ=255)                          tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTc                                                      >  1:2802260/1‑59 (MQ=255)                           gAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCa                                         >  1:648128/1‑71 (MQ=255)                            aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGtt                                                             >  1:1392707/1‑50 (MQ=255)                            aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCa                                         >  1:1807662/1‑70 (MQ=255)                            aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCAc                                        >  1:2553798/1‑71 (MQ=255)                             aaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACg                                       <  1:2091442/71‑1 (MQ=255)                             aaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACg                                       >  1:1519132/1‑71 (MQ=255)                              aaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTc                                                      <  1:1603598/55‑1 (MQ=255)                              aaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGt                                      <  1:994657/71‑1 (MQ=255)                               aCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGtt                                     >  1:1308896/1‑71 (MQ=255)                                cGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                    >  1:3008828/1‑71 (MQ=255)                                cGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                    >  1:2608215/1‑71 (MQ=255)                                 gTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGc                                                         <  1:2817229/49‑1 (MQ=255)                                 gTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCg                                   >  1:2097003/1‑71 (MQ=255)                                   tttGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCa                                         >  1:538441/1‑63 (MQ=255)                                     tGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTc                                          <  1:1798860/60‑1 (MQ=255)                                     tGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCa                                >  1:2085424/1‑70 (MQ=255)                                     tGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCa                                <  1:1602407/70‑1 (MQ=255)                                      ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                               <  1:1370598/70‑1 (MQ=255)                                        aCGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACgg                             >  1:2255062/1‑70 (MQ=255)                                         cGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGc                            >  1:1554372/1‑70 (MQ=255)                                           cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                         <  1:268115/71‑1 (MQ=255)                                           cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                         <  1:900226/71‑1 (MQ=255)                                            ttAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATAc                        >  1:1456122/1‑71 (MQ=255)                                             tAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGt                                      >  1:741279/1‑56 (MQ=255)                                             tAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGc                            >  1:1875888/1‑66 (MQ=255)                                             tAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGc                            >  1:737066/1‑66 (MQ=255)                                             tAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACt                       >  1:903328/1‑71 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACg                                       >  1:2960302/1‑53 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACt                       <  1:1120860/69‑1 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTg                      <  1:1992649/70‑1 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                     <  1:120791/71‑1 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                     <  1:1561129/71‑1 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                     <  1:1255177/71‑1 (MQ=255)                                               aTACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                     >  1:1255337/1‑71 (MQ=255)                                                tACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCt                                                     <  1:271131/38‑1 (MQ=255)                                                tACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTAACTGGCTCACGTTCGCCACGGCATACTGCt                    <  1:2527633/71‑1 (MQ=255)                                                 aCGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                    >  1:581344/1‑54 (MQ=255)                                                 aCGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                    >  1:577900/1‑54 (MQ=255)                                                 aCGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTa                   >  1:538410/1‑71 (MQ=255)                                                  cgcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCa                                         >  1:1647087/1‑48 (MQ=255)                                                  cgcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCa                                         >  1:2614230/1‑48 (MQ=255)                                                  cgcgTCATATACTTGGACCGCACGCGGTCTGCTTCTGTCACt                                               <  1:804968/42‑1 (MQ=255)                                                   gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                               >  1:778985/1‑57 (MQ=255)                                                    cgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGa                <  1:1525709/71‑1 (MQ=255)                                                      tCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                              <  1:2565116/39‑1 (MQ=255)                                                      tCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                 >  1:63479/1‑68 (MQ=255)                                                       cATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGt                                      <  1:1371686/46‑1 (MQ=255)                                                       cATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCa                                <  1:2599164/52‑1 (MQ=255)                                                        atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                              >  1:2325062/1‑37 (MQ=255)                                                        atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAagag            >  1:653217/1‑71 (MQ=255)                                                          ataCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCt                    <  1:1355976/61‑1 (MQ=255)                                                          ataCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt          >  1:2345008/1‑71 (MQ=255)                                                          ataCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGt           <  1:579301/70‑1 (MQ=255)                                                           taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt         <  1:2125168/71‑1 (MQ=255)                                                           taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt         >  1:2034941/1‑71 (MQ=255)                                                            aCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt         <  1:2623518/70‑1 (MQ=255)                                                             cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGaa               >  1:439535/1‑63 (MQ=255)                                                                ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTctct     <  1:327932/70‑1 (MQ=255)                                                                  aCCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTcc  <  1:991622/71‑1 (MQ=255)                                                                     |                                                                    CTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTCC  >  minE/2450230‑2450365

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