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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458                                                |                                               ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              >  1:1476137/1‑50 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:655820/1‑50 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:2568433/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:690332/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:1768885/1‑50 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:747741/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:709372/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  >  1:2016539/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:2621164/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:1771095/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:73671/1‑50 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:1966670/1‑50 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:2500059/50‑1 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  1:1618115/1‑50 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  1:279163/1‑50 (MQ=255)                       ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        >  1:2659550/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  1:2190253/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:1311787/49‑1 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:2216842/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:2166433/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:1958423/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:1420836/1‑50 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:2243299/49‑1 (MQ=255)                                  catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             >  1:1128416/1‑50 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:632691/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:1806400/1‑49 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        >  1:2453214/1‑50 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        >  1:1225592/1‑50 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    >  1:1092877/1‑50 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  1:1783125/50‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:1670891/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:2084711/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:582933/1‑50 (MQ=255)                                                |                                               CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458

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: ‑