Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,824,288	G→A	D452D (GAC→GAT)	hypF ←	carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,824,288	0	G	A	100.0%	204.2 / NA	60	D452D (GAC→GAT)	hypF	carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/34);  total (26/34)

AGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGC  >  minE/1824223‑1824351                                                                  |                                                                aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                            >  1:1062295/1‑71 (MQ=255) aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                            <  1:2585505/71‑1 (MQ=255) aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                            >  1:402257/1‑71 (MQ=255) aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                            <  1:245820/71‑1 (MQ=255)  gTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                            >  1:3048270/1‑70 (MQ=255)   tACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGca                                                          >  1:811708/1‑71 (MQ=255)   tACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGc                                                           <  1:2846733/70‑1 (MQ=255)     ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGcac                                                         >  1:461124/1‑70 (MQ=255)     ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGcac                                                         >  1:2239652/1‑70 (MQ=255)     ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGcac                                                         >  1:615374/1‑70 (MQ=255)     ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGca                                                          >  1:2104696/1‑69 (MQ=255)        tttGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACg                                                     <  1:2751211/71‑1 (MQ=255)           gCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACTTa                                                   <  1:156009/70‑1 (MQ=255)           gCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTa                                                   <  1:3028398/70‑1 (MQ=255)            cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGc                                                           <  1:725730/61‑1 (MQ=255)            cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGc                                                           <  1:1252057/61‑1 (MQ=255)            cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGc                                                           <  1:3190005/61‑1 (MQ=255)            cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTAtt                                                 <  1:844649/71‑1 (MQ=255)            cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTAt                                                  <  1:1301882/70‑1 (MQ=255)                     aTTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTa                                         >  1:1287645/1‑70 (MQ=255)                     aTTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTa                                         >  1:1457573/1‑70 (MQ=255)                     aTTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTa                                         >  1:43738/1‑70 (MQ=255)                        cATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAAt                                     >  1:751100/1‑71 (MQ=255)                         aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAAt                                     <  1:581394/70‑1 (MQ=255)                         aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATc                                    <  1:2797837/71‑1 (MQ=255)                         aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATc                                    <  1:421184/71‑1 (MQ=255)                         aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATc                                    <  1:1708042/71‑1 (MQ=255)                          tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCg                                   <  1:770050/71‑1 (MQ=255)                           ttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTg                                                >  1:2062758/1‑57 (MQ=255)                           ttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTg                                                >  1:805764/1‑57 (MQ=255)                               cgcgCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                              <  1:2523461/71‑1 (MQ=255)                                  gcCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtt                           <  1:1682865/71‑1 (MQ=255)                                  gcCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtt                           <  1:664767/71‑1 (MQ=255)                                  gcCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGt                            <  1:461807/70‑1 (MQ=255)                                   cccACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtt                           <  1:2875397/70‑1 (MQ=255)                                   cccACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtt                           <  1:2655207/70‑1 (MQ=255)                                    ccACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                         <  1:1220867/71‑1 (MQ=255)                                      aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                        >  1:1078209/1‑70 (MQ=255)                                       ctggctggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCa                      <  1:2641141/71‑1 (MQ=255)                                         ggctggAGACATAGCCCGGATGCGCATCATGCACaa                                                       >  1:1835201/1‑36 (MQ=255)                                          gctggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAACTCGTAGATGTTTTGCATCa                   <  1:2016764/71‑1 (MQ=255)                                          gctggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCa                   <  1:2644274/71‑1 (MQ=255)                                            tggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAgg                 >  1:1747159/1‑71 (MQ=255)                                              gagaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCg                                   <  1:2396443/51‑1 (MQ=255)                                              gagaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGc                <  1:2418227/70‑1 (MQ=255)                                              gagaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCg               <  1:2308999/71‑1 (MQ=255)                                              gagaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCg               <  1:1932015/71‑1 (MQ=255)                                                      gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAgg                 >  1:888719/1‑61 (MQ=255)                                                      gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt       <  1:1982523/71‑1 (MQ=255)                                                      gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt       <  1:1978708/71‑1 (MQ=255)                                                       cccGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGa                               <  1:2721603/46‑1 (MQ=255)                                                       cccGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCa                      <  1:2290429/55‑1 (MQ=255)                                                         cGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                    >  1:2312262/1‑55 (MQ=255)                                                         cGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                    >  1:2936642/1‑55 (MQ=255)                                                           gATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCgcgc  >  1:1931022/1‑71 (MQ=255)                                                             tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTa                                 >  1:2536136/1‑38 (MQ=255)                                                             tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTa                                 >  1:2537683/1‑38 (MQ=255)                                                             tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTa                                 >  1:3108251/1‑38 (MQ=255)                                                             tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTa                                 >  1:1853866/1‑38 (MQ=255)                                                             tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTa                                 >  1:2597394/1‑38 (MQ=255)                                                                  |                                                                AGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGC  >  minE/1824223‑1824351

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