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

CTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCG  >  NC_000913/2843369‑2843456                                            |                                             ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                           <  1:1604053/47‑1 (MQ=255)  tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:1288866/1‑47 (MQ=255)    cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                        <  1:13433/47‑1 (MQ=255)    cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                      >  1:1880088/1‑49 (MQ=255)           gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                >  1:1485304/1‑48 (MQ=255)            cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                              >  1:713256/1‑49 (MQ=255)            cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                              >  1:1317209/1‑49 (MQ=255)             gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                              >  1:95999/1‑48 (MQ=255)               gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                           <  1:37725/49‑1 (MQ=255)                cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                           >  1:659098/1‑48 (MQ=255)                  tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacaaca                         >  1:212677/1‑48 (MQ=255)                   ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                       >  1:55116/1‑49 (MQ=255)                   ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                       >  1:132579/1‑49 (MQ=255)                    gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                      <  1:450543/49‑1 (MQ=255)                    gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                      <  1:2010095/49‑1 (MQ=255)                      aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                    >  1:1466629/1‑49 (MQ=255)                      aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                    >  1:1605027/1‑49 (MQ=255)                         gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                  <  1:1828220/48‑1 (MQ=255)                         gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:2056244/1‑49 (MQ=255)                         gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:1462716/1‑49 (MQ=255)                          ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                >  1:39365/1‑49 (MQ=255)                           tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt               >  1:1210400/1‑49 (MQ=255)                               atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc           <  1:1686465/49‑1 (MQ=255)                               atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc           <  1:312566/49‑1 (MQ=255)                               atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc           <  1:447636/49‑1 (MQ=255)                                  aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg         <  1:1067861/48‑1 (MQ=255)                                   aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa        >  1:569742/1‑48 (MQ=255)                                    tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa      <  1:1464309/49‑1 (MQ=255)                                       aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc   >  1:1503640/1‑49 (MQ=255)                                       aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc   <  1:1449716/49‑1 (MQ=255)                                        tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg  >  1:1309671/1‑49 (MQ=255)                                        tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg  >  1:742568/1‑49 (MQ=255)                                            |                                             CTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCG  >  NC_000913/2843369‑2843456

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