Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	801,414	A→G	S186S (AGT→AGC)	dhaL ←	dihydroxyacetone kinase, C‑terminal domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	801,414	0	A	G	100.0%	202.3 / NA	55	S186S (AGT→AGC)	dhaL	dihydroxyacetone kinase, C‑terminal domain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/24);  total (31/24)

CTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGC  >  minE/801353‑801479                                                              |                                                                  cTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCa                                                          >  1:1796549/1‑71 (MQ=255) cTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCa                                                          <  1:3058868/71‑1 (MQ=255) cTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCa                                                          >  1:1209737/1‑71 (MQ=255) cTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCa                                                          <  1:935047/71‑1 (MQ=255)   aaCGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTg                                                                >  1:1829151/1‑63 (MQ=255)   aaCGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCAcc                                                        >  1:364972/1‑71 (MQ=255)    aCGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCg                                                       <  1:2367435/71‑1 (MQ=255)     cGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCAcc                                                        >  1:1165141/1‑69 (MQ=255)      gCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCAcc                                                        >  1:168441/1‑68 (MQ=255)        cAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGAt                                                   <  1:657642/71‑1 (MQ=255)         aaCATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGAt                                                   >  1:2305713/1‑70 (MQ=255)          aCATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATaa                                                 >  1:1204195/1‑71 (MQ=255)           cATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAAc                                                >  1:3067655/1‑71 (MQ=255)               tGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTgg                                             >  1:710709/1‑70 (MQ=255)               tGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGgc                                            >  1:1606639/1‑71 (MQ=255)                gCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGgc                                            <  1:29773/70‑1 (MQ=255)                  atcatAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGcgcg                                         <  1:1673884/71‑1 (MQ=255)                  atcatAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGcgcg                                         <  1:1783208/71‑1 (MQ=255)                    catAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGc                                       >  1:107393/1‑71 (MQ=255)                        aaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                   <  1:3100840/71‑1 (MQ=255)                        aaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCtt                                    <  1:2003976/70‑1 (MQ=255)                        aaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCtt                                    <  1:1531624/70‑1 (MQ=255)                         aCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                   >  1:1841830/1‑70 (MQ=255)                              aCCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGc                             <  1:1993810/71‑1 (MQ=255)                                cGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCtt                           >  1:674233/1‑71 (MQ=255)                                  aGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTCCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTg                          >  1:2991854/1‑70 (MQ=255)                                  aGGTCGCGCCGGGATCATCGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTg                          >  1:2637767/1‑70 (MQ=255)                                   ggTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCa                        <  1:1531110/71‑1 (MQ=255)                                    gTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCat                       <  1:3130759/71‑1 (MQ=255)                                       gcgcCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCgg                               <  1:2709417/60‑1 (MQ=255)                                         gcCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGc                             >  1:2749795/1‑60 (MQ=255)                                          ccGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaat                 >  1:625776/1‑71 (MQ=255)                                           cGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACt                                               >  1:38933/1‑40 (MQ=255)                                           cGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACt                                               >  1:638110/1‑40 (MQ=255)                                            gggATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGCGCTTGCatcgtaatcg               >  1:1182233/1‑71 (MQ=255)                                             ggATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaatcg               >  1:1183897/1‑70 (MQ=255)                                              gATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaatcgta             <  1:1122063/71‑1 (MQ=255)                                               aTCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGc                                 >  1:579511/1‑50 (MQ=255)                                               aTCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCa                        <  1:689548/59‑1 (MQ=255)                                               aTCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTAc            >  1:1236630/1‑71 (MQ=255)                                                 ccTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaatcgta             >  1:1779047/1‑68 (MQ=255)                                                 ccTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACtt          >  1:1903195/1‑71 (MQ=255)                                                   tGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgta                   <  1:2791362/60‑1 (MQ=255)                                                   tGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgta                   <  1:2439815/60‑1 (MQ=255)                                                    ggTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCat                       <  1:2005814/55‑1 (MQ=255)                                                     gTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTtgtg      >  1:74947/1‑71 (MQ=255)                                                      tGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                   <  1:3171243/41‑1 (MQ=255)                                                      tGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTg                                  >  1:261720/1‑42 (MQ=255)                                                      tGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTg                                  >  1:579721/1‑42 (MQ=255)                                                       gACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                   <  1:1189069/40‑1 (MQ=255)                                                       gACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaat                 >  1:2929015/1‑58 (MQ=255)                                                       gACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACt           <  1:2418803/64‑1 (MQ=255)                                                        aCCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCa    <  1:1476077/70‑1 (MQ=255)                                                        aCCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCa    <  1:2130256/70‑1 (MQ=255)                                                         ccAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGc  >  1:303112/1‑71 (MQ=255)                                                              |                                                                  CTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGC  >  minE/801353‑801479

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