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

CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843366‑2843457                                               |                                              ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgga                                            >  1:3122236/1‑49 (MQ=255) ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                            <  1:2658595/50‑1 (MQ=255)   gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           >  1:3548945/1‑49 (MQ=255)    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           <  1:3808563/48‑1 (MQ=255)    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           <  1:3467116/48‑1 (MQ=255)    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:116588/50‑1 (MQ=255)       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                        <  1:398236/48‑1 (MQ=255)       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                        <  1:1219309/48‑1 (MQ=255)       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                       <  1:854878/49‑1 (MQ=255)        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                     <  1:1872475/50‑1 (MQ=255)        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                     <  1:2861522/50‑1 (MQ=255)           tCGGCGGGCGTGGCATCGTTCTCATCCATAAATTTATGGCTCAGTTGACt                                  >  1:2299578/1‑50 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                <  1:1941606/50‑1 (MQ=255)              gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:1568211/1‑50 (MQ=255)                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            <  1:3841718/50‑1 (MQ=255)                   cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          <  1:1706366/50‑1 (MQ=255)                    gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         >  1:140757/1‑50 (MQ=255)                      ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                       >  1:1189569/1‑50 (MQ=255)                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                     <  1:1152704/50‑1 (MQ=255)                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                     >  1:71734/1‑50 (MQ=255)                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    <  1:1831854/50‑1 (MQ=255)                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:1811764/1‑50 (MQ=255)                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:335273/1‑49 (MQ=255)                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                   <  1:3677797/50‑1 (MQ=255)                           cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  >  1:2416333/1‑50 (MQ=255)                           cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  >  1:403877/1‑50 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 <  1:411057/50‑1 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:3529312/1‑50 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:1962697/1‑50 (MQ=255)                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                <  1:4027869/50‑1 (MQ=255)                              tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:3554705/1‑49 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              <  1:1289780/47‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg           <  1:3374015/50‑1 (MQ=255)                                     aaTAAATTTATGGTTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:2363678/1‑49 (MQ=255)                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:1972131/1‑49 (MQ=255)                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:4217713/1‑49 (MQ=255)                                         aaTTTATGGCTCAGTTGACTGAACACCGCCTTTTCACTCATCGGATAcgc    <  1:1875617/50‑1 (MQ=255)                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  <  1:2567691/50‑1 (MQ=255)                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  >  1:1662633/1‑50 (MQ=255)                                               |                                              CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843366‑2843457

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