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

GCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843368‑2843457                                             |                                              gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           >  1:1891877/1‑49 (MQ=255)  ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           <  1:642830/48‑1 (MQ=255)     cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                       <  1:2781755/49‑1 (MQ=255)     cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                       <  1:1568197/49‑1 (MQ=255)      gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                     <  1:168904/50‑1 (MQ=255)           ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 >  1:2404883/1‑49 (MQ=255)           ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 <  1:1409038/49‑1 (MQ=255)           ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 <  1:2172260/49‑1 (MQ=255)              gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                             <  1:698995/50‑1 (MQ=255)                 cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          <  1:1071807/50‑1 (MQ=255)                 cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          >  1:2366700/1‑50 (MQ=255)                 cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          <  1:2324083/50‑1 (MQ=255)                  gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         >  1:549465/1‑50 (MQ=255)                  gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         >  1:259349/1‑50 (MQ=255)                   tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                        >  1:1035154/1‑50 (MQ=255)                   tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                        >  1:1394805/1‑50 (MQ=255)                   tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                        <  1:2685204/50‑1 (MQ=255)                      cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                     >  1:2056151/1‑50 (MQ=255)                      cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                     >  1:1501117/1‑50 (MQ=255)                       aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    <  1:2250043/50‑1 (MQ=255)                          gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 <  1:1144006/50‑1 (MQ=255)                          gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:573459/1‑50 (MQ=255)                           ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:1955333/1‑50 (MQ=255)                           ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                <  1:1178688/50‑1 (MQ=255)                           ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:2765557/1‑50 (MQ=255)                             ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              >  1:2738572/1‑50 (MQ=255)                               catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc            >  1:704774/1‑50 (MQ=255)                                atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt             <  1:2333759/48‑1 (MQ=255)                                atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg           <  1:2740103/50‑1 (MQ=255)                                 tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg          <  1:1920601/50‑1 (MQ=255)                                 tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg           >  1:939388/1‑49 (MQ=255)                                  caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:1993894/1‑50 (MQ=255)                                   aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:543581/1‑49 (MQ=255)                                   aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt        <  1:1920985/50‑1 (MQ=255)                                    aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa       <  1:1661957/50‑1 (MQ=255)                                     tAAAGTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc      >  1:1926627/1‑50 (MQ=255)                                         tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  >  1:1774295/1‑50 (MQ=255)                                             |                                              GCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843368‑2843457

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