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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843365‑2843457                                                |                                              ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                             <  1:6546011/50‑1 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           >  1:3629915/1‑50 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGcc                                           >  1:356344/1‑48 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGcc                                           >  1:6803365/1‑48 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           <  1:4075283/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                          <  1:310076/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                          <  1:5342831/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                          <  1:2229956/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:2505004/50‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:4582285/50‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:6018512/50‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:63568/50‑1 (MQ=255)      tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                        >  1:4570444/1‑50 (MQ=255)       gcgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                       >  1:5323808/1‑50 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                     <  1:6587991/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                   <  1:312778/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                   <  1:4532495/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 <  1:4379361/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 <  1:659937/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 <  1:248798/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 >  1:2890065/1‑49 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               <  1:2580774/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               <  1:2691108/50‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:690315/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:3380643/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:5783680/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:2028651/1‑49 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                             >  1:3753773/1‑50 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            >  1:3460609/1‑49 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            >  1:2319935/1‑49 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          >  1:837972/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         <  1:863068/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         <  1:6216452/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         >  1:4673683/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         >  1:1511627/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         <  1:6241166/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                         <  1:1790849/50‑1 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                        >  1:3206403/1‑50 (MQ=255)                       ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                       >  1:3448461/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                      <  1:3492053/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                      <  1:5107877/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                      <  1:6803583/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                       <  1:5809988/49‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                     >  1:4076168/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:1138560/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:828935/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:4243572/1‑50 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                   <  1:5917564/50‑1 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  >  1:5677176/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  <  1:1782891/49‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:2205064/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 <  1:1670616/50‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 <  1:2367863/50‑1 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                <  1:6078988/50‑1 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:4729057/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                <  1:3121651/50‑1 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                <  1:85526/50‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc               <  1:2340640/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              >  1:288758/1‑50 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              >  1:412629/1‑50 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              >  1:863671/1‑50 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt             >  1:5576701/1‑50 (MQ=255)                                 ccatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt             >  1:2921278/2‑50 (MQ=255)                                  catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc            >  1:5885630/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc            <  1:3647923/49‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCGCCTTTTCACTCATCgg          <  1:2104830/50‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg           >  1:4010039/1‑49 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:3425877/1‑50 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:2515153/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         <  1:2810490/49‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         <  1:4530234/49‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt        <  1:1767863/50‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt        <  1:1857352/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa       >  1:6323473/1‑50 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa       <  1:3639411/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa       <  1:2843182/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc      <  1:4313996/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc      <  1:4870898/50‑1 (MQ=255)                                         aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg     >  1:6194296/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc    <  1:2136455/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc    <  1:5947847/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc    >  1:4948803/1‑50 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCGCCTTTTCACTCATCGGATAcgcg   <  1:2007898/50‑1 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg   <  1:5016339/50‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  <  1:3568957/50‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  >  1:2008753/1‑50 (MQ=255)                                                |                                              CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843365‑2843457

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