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%	136.1 / NA	35	R12H (CGT→CAT)	hycH	formate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/16);  total (19/16)

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACG  >  NC_000913/2843365‑2843454                                                |                                           ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                          <  1:1293335/50‑1 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                        >  1:980402/1‑49 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                        >  1:491187/1‑49 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                    <  1:1907062/49‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                    <  1:182807/49‑1 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                  <  1:1702717/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                              <  1:1341130/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                             <  1:1536529/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                             <  1:1491973/50‑1 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                     <  1:325635/50‑1 (MQ=255)                       ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                    >  1:978373/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                   <  1:1312427/50‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                  >  1:1667900/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                 <  1:797433/50‑1 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                 >  1:2184891/1‑50 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa               >  1:262251/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc              >  1:1473340/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc              >  1:1696355/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc              <  1:2050629/50‑1 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt             >  1:373734/1‑49 (MQ=255)                                  catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc         >  1:1414601/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc         <  1:689443/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg        <  1:924495/50‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg       <  1:1207261/50‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg        >  1:1887093/1‑49 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg        >  1:1141408/1‑49 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa      >  1:2185362/1‑50 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa      >  1:449644/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa      >  1:1494457/1‑49 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt     <  1:1325679/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa    >  1:1840071/1‑50 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACAAAACCTTTTCACTCATAGGAt     >  1:1210793/1‑49 (MQ=255)                                       aTAAATTGATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa    <  1:1454592/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc   >  1:640045/1‑50 (MQ=255)                                         aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg  >  1:327034/1‑50 (MQ=255)                                                |                                           CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACG  >  NC_000913/2843365‑2843454

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