Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,843,412	C→T	R12H (CGT→CAT)	hycH ←	formate hydrogenlyase assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,843,412	0	C	T	100.0%	263.4 / NA	66	R12H (CGT→CAT)	hycH	formate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/33);  total (33/33)

TGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843367‑2843458                                              |                                               tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:4251822/1‑50 (MQ=255)   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:568958/48‑1 (MQ=255)   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:1478981/1‑48 (MQ=255)   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:747697/49‑1 (MQ=255)   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:3013185/50‑1 (MQ=255)   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:391606/50‑1 (MQ=255)      cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         <  1:1654482/48‑1 (MQ=255)      cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:312312/50‑1 (MQ=255)      cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:2901417/49‑1 (MQ=255)       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:3300735/1‑49 (MQ=255)       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:3142014/50‑1 (MQ=255)           cggcggGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  1:2951651/1‑49 (MQ=255)            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:365159/48‑1 (MQ=255)            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:1329627/49‑1 (MQ=255)            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:4406036/49‑1 (MQ=255)            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  1:4446540/50‑1 (MQ=255)              cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:1104994/1‑49 (MQ=255)               gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:533078/50‑1 (MQ=255)               gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:516580/50‑1 (MQ=255)                ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:3024221/1‑50 (MQ=255)                ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:73791/50‑1 (MQ=255)                 gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:1108728/1‑49 (MQ=255)                  cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  1:2314856/1‑50 (MQ=255)                   gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  1:3327562/50‑1 (MQ=255)                   gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:3227329/1‑50 (MQ=255)                   gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:4461482/1‑50 (MQ=255)                    tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  1:4364945/50‑1 (MQ=255)                      gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:1078434/1‑50 (MQ=255)                      gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:374812/1‑50 (MQ=255)                      gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:3341285/1‑50 (MQ=255)                      gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:3672695/49‑1 (MQ=255)                      gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:247475/49‑1 (MQ=255)                       cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:4656228/50‑1 (MQ=255)                        aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:4742331/1‑50 (MQ=255)                          cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:1830204/1‑50 (MQ=255)                          cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:2838289/1‑50 (MQ=255)                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:1993589/49‑1 (MQ=255)                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:3401081/1‑50 (MQ=255)                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:278883/1‑50 (MQ=255)                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:2792800/1‑50 (MQ=255)                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:296042/1‑50 (MQ=255)                            ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:1255664/1‑50 (MQ=255)                            ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:3281294/1‑50 (MQ=255)                            ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  1:145438/50‑1 (MQ=255)                            ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  1:1984010/50‑1 (MQ=255)                             tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:3905743/1‑49 (MQ=255)                             tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4441806/1‑49 (MQ=255)                              ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:2951698/49‑1 (MQ=255)                                catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             >  1:3295947/1‑50 (MQ=255)                                catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             >  1:4244651/1‑50 (MQ=255)                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               <  1:2870459/47‑1 (MQ=255)                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:2179374/48‑1 (MQ=255)                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:1076359/49‑1 (MQ=255)                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:3940035/50‑1 (MQ=255)                                  tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:129837/50‑1 (MQ=255)                                  tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:2073766/1‑49 (MQ=255)                                   caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:692095/1‑50 (MQ=255)                                   caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:2899339/1‑50 (MQ=255)                                    aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:1512948/1‑49 (MQ=255)                                    aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:572217/1‑49 (MQ=255)                                    aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  1:1970595/50‑1 (MQ=255)                                     aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:531738/50‑1 (MQ=255)                                      tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:557738/50‑1 (MQ=255)                                       aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  1:4953782/1‑50 (MQ=255)                                          tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  1:4137429/50‑1 (MQ=255)                                           ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:1061332/1‑50 (MQ=255)                                              |                                               TGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843367‑2843458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑