Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,509,289	T→C	intergenic (+34/+19)	yicI → / ← yicH	predicted alpha‑glucosidase/conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,509,289	0	T	C	100.0%	240.1 / NA	65	intergenic (+34/+19)	yicI/yicH	predicted alpha‑glucosidase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/30);  total (35/30)

AGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGC  >  minE/2509221‑2509354                                                                     |                                                                  aGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCAc                                                                 >  1:3383837/1‑71 (MQ=255) aGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCAc                                                                 >  1:244316/1‑71 (MQ=255)   ccTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCAccc                                                               <  1:326162/71‑1 (MQ=255)    cTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCg                                                              >  1:1669188/1‑71 (MQ=255)      cAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGc                                                            >  1:47293/1‑71 (MQ=255)       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCAc                                                                 >  1:1526523/1‑65 (MQ=255)       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGc                                                            >  1:2978894/1‑70 (MQ=255)        aGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCga                                                          <  1:2300462/71‑1 (MQ=255)             aTGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTg                                                     <  1:790871/71‑1 (MQ=255)              tGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGa                                                    <  1:1541893/71‑1 (MQ=255)              tGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGa                                                    <  1:2193764/71‑1 (MQ=255)              tGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGa                                                    <  1:642132/71‑1 (MQ=255)              tGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGa                                                    <  1:151254/71‑1 (MQ=255)              tGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGa                                                    <  1:3157094/71‑1 (MQ=255)               gcgcTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGAt                                                   >  1:1294006/1‑71 (MQ=255)                 gcTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGt                                                 <  1:3824813/71‑1 (MQ=255)                  cTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCgag                                                         <  1:147034/62‑1 (MQ=255)                    gACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTa                                              >  1:2257130/1‑71 (MQ=255)                     aCAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTa                                              <  1:1981785/70‑1 (MQ=255)                     aCAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTAc                                             <  1:838410/71‑1 (MQ=255)                     aCAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTAc                                             >  1:402869/1‑71 (MQ=255)                     aCAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTAc                                             >  1:1570358/1‑71 (MQ=255)                     aCAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGACGCACCCGGCGAGATTGATGTTTa                                              <  1:2021246/70‑1 (MQ=255)                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCg                                                              >  1:997903/1‑52 (MQ=255)                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCg                                                              >  1:1009601/1‑52 (MQ=255)                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCg                                           >  1:2093965/1‑71 (MQ=255)                        aTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGc                                          <  1:4083179/71‑1 (MQ=255)                         ttACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCa                                         >  1:1661516/1‑71 (MQ=255)                         ttACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCa                                         >  1:138625/1‑71 (MQ=255)                          tACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTg                                                     >  1:3066003/1‑58 (MQ=255)                          tACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCAc                                        >  1:617249/1‑71 (MQ=255)                             gTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGAt                                                       <  1:283594/53‑1 (MQ=255)                                 tAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTc                                  <  1:2603140/70‑1 (MQ=255)                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCg                                <  1:1227729/70‑1 (MQ=255)                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGt                               <  1:1861320/71‑1 (MQ=255)                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGt                               <  1:3443522/71‑1 (MQ=255)                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGt                               >  1:857056/1‑71 (MQ=255)                                     cTATTGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTc                              >  1:1895560/1‑70 (MQ=255)                                        ttGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGAt                                                   >  1:1863035/1‑46 (MQ=255)                                        ttGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTc                              <  1:2128378/67‑1 (MQ=255)                                        ttGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTc                              <  1:938962/67‑1 (MQ=255)                                        ttGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATg                          >  1:1530276/1‑71 (MQ=255)                                        ttGATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATg                          >  1:4097470/1‑71 (MQ=255)                                          gATGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGt                         <  1:2555280/70‑1 (MQ=255)                                            tGGAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAAc                      <  1:790660/71‑1 (MQ=255)                                              gAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACtt                    >  1:1957003/1‑71 (MQ=255)                                              gAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACt                     <  1:2308700/70‑1 (MQ=255)                                              gAATGCCGGGTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACt                     <  1:1609827/70‑1 (MQ=255)                                                      ggTACGTCTGAGTCGCGCCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTg            >  1:2190311/1‑71 (MQ=255)                                                      ggTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGa                            >  1:2636073/1‑55 (MQ=255)                                                       gTACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTg            >  1:2654993/1‑70 (MQ=255)                                                        tACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTgg           >  1:3624239/1‑70 (MQ=255)                                                        tACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          <  1:563763/71‑1 (MQ=255)                                                        tACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          <  1:2694152/71‑1 (MQ=255)                                                        tACGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          <  1:4026304/71‑1 (MQ=255)                                                        tACGTCTGAGTCGCACCCGGCGAGATTGACGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          <  1:3651450/71‑1 (MQ=255)                                                         aCGTCTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          >  1:1409584/1‑70 (MQ=255)                                                             cTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGa     >  1:3523011/1‑71 (MQ=255)                                                             cTGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGa     >  1:2394054/1‑71 (MQ=255)                                                              tGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTc                   <  1:3248502/56‑1 (MQ=255)                                                              tGAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGa     >  1:3200515/1‑70 (MQ=255)                                                               gAGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATg   >  1:4153370/1‑71 (MQ=255)                                                                aGTCGCCCCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGc  >  1:4122243/1‑71 (MQ=255)                                                                aGTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGa     >  1:2684947/1‑68 (MQ=255)                                                                 gTCGCACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGc  >  1:1595436/1‑70 (MQ=255)                                                                     |                                                                  AGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGC  >  minE/2509221‑2509354

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