Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,330,586	A→G	Y334H (TAC→CAC)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,330,586	0	A	G	100.0%	187.5 / NA	53	Y334H (TAC→CAC)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/25);  total (28/25)

CACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGC  >  minE/1330521‑1330652                                                                  |                                                                   cacCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggt                                                               <  1:799603/71‑1 (MQ=255)  acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                              >  1:2864469/1‑71 (MQ=255)      gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                          >  1:1193981/1‑71 (MQ=255)         cGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCt                                                        <  1:1152723/70‑1 (MQ=255)          ggCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCa                                                      <  1:993239/71‑1 (MQ=255)             gACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCAttt                                                   >  1:2158525/1‑71 (MQ=255)              aCCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTc                                                  <  1:2660703/71‑1 (MQ=255)                cagcagCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtgg                                                             >  1:476146/1‑58 (MQ=255)                cagcagCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtgg                                                             >  1:414498/1‑58 (MQ=255)                cagcagCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCg                                                 >  1:1764593/1‑70 (MQ=255)                  gcagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                              >  1:2190203/1‑55 (MQ=255)                  gcagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCt                                                        <  1:1098244/61‑1 (MQ=255)                  gcagcCACGGGATAATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGcca                                              >  1:2413402/1‑71 (MQ=255)                      ccACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTc                                                       <  1:931233/58‑1 (MQ=255)                       cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCg                                                 >  1:1895607/1‑63 (MQ=255)                        aCGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTggg                                                            <  1:1821639/51‑1 (MQ=255)                        aCGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGccacca                                           >  1:1359389/1‑68 (MQ=255)                         cGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTc                                                  >  1:735036/1‑60 (MQ=255)                         cGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGa                                       <  1:861307/71‑1 (MQ=255)                          gggATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggt                                                               <  1:1823690/46‑1 (MQ=255)                          gggATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggt                                                               >  1:656461/1‑46 (MQ=255)                          gggATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggt                                                               <  1:2444077/46‑1 (MQ=255)                          gggATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                                         <  1:2703631/68‑1 (MQ=255)                            gATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATg                                     <  1:2370809/70‑1 (MQ=255)                             aTCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGcc                                               <  1:287265/59‑1 (MQ=255)                             aTCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGa                                       <  1:2274188/67‑1 (MQ=255)                             aTCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGa                                       <  1:2277823/67‑1 (MQ=255)                                aTGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggt                                                               >  1:1939817/1‑40 (MQ=255)                                aTGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTc                                >  1:103007/1‑71 (MQ=255)                                    cGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTcc                            >  1:1199849/1‑71 (MQ=255)                                      cTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCtt                          >  1:1913089/1‑71 (MQ=255)                                      cTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCtt                          >  1:1974885/1‑71 (MQ=255)                                       tGTCGCAGTAAGCGAAATCACGGAAGTGGTGGTGGGCGCTCATTTCGCCACCATAGATgg                                    >  1:235280/1‑60 (MQ=255)                                           gCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGAt                                      <  1:1563325/54‑1 (MQ=255)                                           gCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCa                     >  1:2352954/1‑71 (MQ=255)                                              gTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGcca                                              <  1:2259143/43‑1 (MQ=255)                                              gTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATAc                  >  1:1120117/1‑71 (MQ=255)                                               tAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGa                                       >  1:2206223/1‑49 (MQ=255)                                                aaGCGAAATCACGGAAGTGATGGTGGGCGCTCAttt                                                   >  1:2693349/1‑36 (MQ=255)                                                aaGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGt                >  1:1231962/1‑71 (MQ=255)                                                 aGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCg                                  <  1:1866338/52‑1 (MQ=255)                                                 aGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGtt               <  1:2520231/71‑1 (MQ=255)                                                  gCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTc              <  1:2139707/71‑1 (MQ=255)                                                  gCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTc              <  1:1899744/71‑1 (MQ=255)                                                     aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt            >  1:1586661/1‑70 (MQ=255)                                                          aCGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                                         <  1:586964/36‑1 (MQ=255)                                                           cGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATaaa     <  1:2166191/71‑1 (MQ=255)                                                            ggAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATg                                     >  1:1522041/1‑38 (MQ=255)                                                            ggAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCg                                  <  1:3026379/41‑1 (MQ=255)                                                            ggAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAg    >  1:2533392/1‑71 (MQ=255)                                                             gAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt             <  1:1332782/61‑1 (MQ=255)                                                              aaGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGc  >  1:608598/1‑71 (MQ=255)                                                              aaGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGc  >  1:1564363/1‑71 (MQ=255)                                                                  |                                                                   CACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGC  >  minE/1330521‑1330652

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