Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,550,360	A→G	K77K (AAA→AAG)	tdh →	threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,550,360	0	A	G	100.0%	291.6 / NA	80	K77K (AAA→AAG)	tdh	threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (45/35);  total (45/35)

TCCCGGTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAAGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTG  >  minE/2550293‑2550425                                                                    |                                                                  tCCCGGTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGc                                                                <  1:1885024/71‑1 (MQ=255)  cccGGTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCt                                                               >  1:226873/1‑71 (MQ=255)  cccGGTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCt                                                               >  1:1506005/1‑71 (MQ=255)     ggTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTc                                                             >  1:605481/1‑70 (MQ=255)         ccGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTc                                                             >  1:1448804/1‑66 (MQ=255)           gATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTc                                                             >  1:1733948/1‑64 (MQ=255)           gATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCg                                                      <  1:1493880/71‑1 (MQ=255)               gTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTc                                                             <  1:447848/60‑1 (MQ=255)               gTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCg                                                   <  1:1996216/70‑1 (MQ=255)                  gTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcg                                               <  1:3013965/71‑1 (MQ=255)                   tGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcgc                                              >  1:2527848/1‑71 (MQ=255)                   tGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATc                                                <  1:3004225/69‑1 (MQ=255)                   tGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATc                                                <  1:1345151/69‑1 (MQ=255)                    gggCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcgcg                                             >  1:856416/1‑71 (MQ=255)                     ggCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            >  1:2051071/1‑71 (MQ=255)                       ccATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGtt                                           >  1:764536/1‑70 (MQ=255)                         aTGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCt                                        >  1:2655288/1‑71 (MQ=255)                         aTGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCt                                        >  1:466533/1‑71 (MQ=255)                         aTGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCt                                        >  1:1087549/1‑71 (MQ=255)                         aTGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCt                                        >  1:683834/1‑71 (MQ=255)                            aaTATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGc                                     >  1:1139335/1‑71 (MQ=255)                             atatGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCg                                    <  1:1769957/71‑1 (MQ=255)                              tatGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCg                                    >  1:1346136/1‑70 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGTAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:2478873/58‑1 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:3143873/58‑1 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:1422080/58‑1 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:3149757/58‑1 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:737628/1‑71 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:126807/1‑71 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:123093/1‑71 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:1591180/1‑71 (MQ=255)                                  tCGGTGAAGTGGTAGGTATTGGACAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:2312475/58‑1 (MQ=255)                                   cGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCgg                                                     >  1:819207/1‑48 (MQ=255)                                   cGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:931592/1‑70 (MQ=255)                                   cGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:245899/1‑70 (MQ=255)                                   cGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGc                               >  1:137647/1‑70 (MQ=255)                                    ggTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAg                                 >  1:2221823/1‑67 (MQ=255)                                     gTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGAt                                                 >  1:1404714/1‑50 (MQ=255)                                     gTGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGAt                                                 >  1:213330/1‑50 (MQ=255)                                      tGAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCata                           >  1:325089/1‑71 (MQ=255)                                       gAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCatat                          >  1:2967819/1‑71 (MQ=255)                                       gAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCata                           >  1:1562854/1‑70 (MQ=255)                                       gAAGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCata                           >  1:2764990/1‑70 (MQ=255)                                         aGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCa                        >  1:2965522/1‑71 (MQ=255)                                         aGTGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCa                        >  1:2551320/1‑71 (MQ=255)                                           tGGTAGGTATTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcg                                               <  1:2174901/46‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcgcg                                             >  1:239776/1‑46 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGa                                   <  1:1074905/56‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGa                                   <  1:1981318/56‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCAc                       <  1:1228703/68‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCAc                       <  1:826857/68‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCAc                       <  1:1899718/68‑1 (MQ=255)                                             gtaggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCt                     <  1:2611607/70‑1 (MQ=255)                                              taggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCtgt                   <  1:2134579/71‑1 (MQ=255)                                                ggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTgg                 >  1:2646213/1‑71 (MQ=255)                                                ggtaTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTgg                 <  1:532456/71‑1 (MQ=255)                                                  taTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcgcg                                             <  1:1327421/41‑1 (MQ=255)                                                  taTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATcgcg                                             <  1:2724669/41‑1 (MQ=255)                                                  taTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTc               <  1:2632860/71‑1 (MQ=255)                                                   aTTGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCa              >  1:2355171/1‑71 (MQ=255)                                                    ttGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGtt                                           <  1:2378132/41‑1 (MQ=255)                                                     tGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGtt                                           <  1:2487274/40‑1 (MQ=255)                                                     tGGTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCAtt            >  1:273203/1‑71 (MQ=255)                                                      ggTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:1485780/38‑1 (MQ=255)                                                      ggTCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGt                                            <  1:1985653/38‑1 (MQ=255)                                                        tCAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGcc         <  1:3107451/71‑1 (MQ=255)                                                         cAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTc               >  1:1575779/1‑64 (MQ=255)                                                         cAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTc               >  1:2351014/1‑64 (MQ=255)                                                         cAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTc               >  1:2173125/1‑64 (MQ=255)                                                         cAGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTc               >  1:2630342/1‑64 (MQ=255)                                                          aGGAAGTGAAGGGCTTCAAGATCGGCGCTCGCGTTTCTGGCGAAGGCCATATCa                        >  1:3115831/1‑54 (MQ=255)                                                          aGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCatat                          >  1:199394/1‑52 (MQ=255)                                                          aGGAAGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTgg                 <  1:735684/61‑1 (MQ=255)                                                             aaGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTgg                                      >  1:2530115/1‑37 (MQ=255)                                                             aaGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTgg                                      >  1:3081806/1‑37 (MQ=255)                                                             aaGTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTgg                                      >  1:2582611/1‑37 (MQ=255)                                                               gTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTg  <  1:1632896/71‑1 (MQ=255)                                                               gTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTg  <  1:1847775/71‑1 (MQ=255)                                                               gTGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTg  <  1:1373209/71‑1 (MQ=255)                                                                tGAAGGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTg  <  1:2648861/70‑1 (MQ=255)                                                                    |                                                                  TCCCGGTGCCGATGGTCGTGGGCCATGAATATGTCGGTGAAGTGGTAGGTATTGGTCAGGAAGTGAAAGGCTTCAAGATCGGCGATCGCGTTTCTGGCGAAGGCCATATCACCTGTGGTCATTGCCGCAACTG  >  minE/2550293‑2550425

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