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

CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843366‑2843458                                               |                                               ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             >  1:2278685/1‑50 (MQ=255)    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:3630905/48‑1 (MQ=255)    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:7216273/49‑1 (MQ=255)     tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         >  1:4635201/1‑50 (MQ=255)     tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         >  1:3139333/1‑50 (MQ=255)      gcgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  1:6230632/1‑50 (MQ=255)      gcgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  1:5611854/1‑50 (MQ=255)       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:8473044/50‑1 (MQ=255)       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:7049877/49‑1 (MQ=255)        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:4298948/1‑49 (MQ=255)        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:1999318/1‑49 (MQ=255)            cggcggGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  1:1724492/1‑49 (MQ=255)            cggcggGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  1:689932/1‑49 (MQ=255)            cggcggGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  1:5626399/1‑49 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:1072533/48‑1 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:3454534/48‑1 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:2390013/48‑1 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:2125600/48‑1 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  >  1:4796397/1‑49 (MQ=255)             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  1:1664064/50‑1 (MQ=255)               cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:4551503/50‑1 (MQ=255)               cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:7246506/50‑1 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:6009128/50‑1 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:8436095/50‑1 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:6653805/1‑50 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:7869788/50‑1 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:8263560/1‑50 (MQ=255)                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:3356915/1‑50 (MQ=255)                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:2071008/1‑50 (MQ=255)                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  1:772400/50‑1 (MQ=255)                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  1:8398569/50‑1 (MQ=255)                   cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  1:6500487/1‑50 (MQ=255)                   cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           <  1:7183879/50‑1 (MQ=255)                    gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:4389769/1‑50 (MQ=255)                      ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        >  1:5767824/1‑50 (MQ=255)                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:3766359/1‑50 (MQ=255)                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:5516027/1‑50 (MQ=255)                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:5556429/1‑50 (MQ=255)                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:4649918/49‑1 (MQ=255)                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:8382466/50‑1 (MQ=255)                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  1:1381977/50‑1 (MQ=255)                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:1904213/1‑50 (MQ=255)                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:1092693/1‑49 (MQ=255)                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:8268883/1‑49 (MQ=255)                           cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:1084985/1‑50 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:5586014/49‑1 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:8170273/49‑1 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  <  1:6371397/50‑1 (MQ=255)                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:4373571/1‑50 (MQ=255)                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:2438499/1‑50 (MQ=255)                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  1:3074278/50‑1 (MQ=255)                              tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:8258910/1‑49 (MQ=255)                               ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               >  1:722030/1‑50 (MQ=255)                               ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               >  1:8538837/1‑50 (MQ=255)                               ctcATCAATAAATTTATGGCTCAATTGACTGAACACCGCCTTTTCACTc                <  1:2028243/49‑1 (MQ=38)                                tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:3867493/1‑50 (MQ=255)                                tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:4257101/1‑50 (MQ=255)                                 catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             >  1:6500460/1‑50 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               <  1:4031159/47‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:8518849/48‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:8270600/48‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:6077229/49‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:7571968/49‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:7532233/49‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:5990122/49‑1 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:1004572/1‑50 (MQ=255)                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:590764/50‑1 (MQ=255)                                   tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:4787467/1‑49 (MQ=255)                                    caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:3393131/1‑50 (MQ=255)                                     aaTAAATTTATGGCTCGGTTGACTGAACACCACCTTTTCACTCATCGGa          <  1:3567656/49‑1 (MQ=255)                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:6989856/1‑49 (MQ=255)                                      aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:5217155/50‑1 (MQ=255)                                      aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        >  1:1231076/1‑50 (MQ=255)                                       tAAATTTATGGCTCAGTTTACTGAACACCACCTTTTCACTCATCGGATAc       >  1:2300454/1‑50 (MQ=255)                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:6349650/50‑1 (MQ=255)                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:4003899/1‑50 (MQ=255)                                        aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  1:367930/1‑50 (MQ=255)                                        aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  1:7997902/1‑50 (MQ=255)                                        aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  1:7197466/1‑50 (MQ=255)                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:7600762/1‑50 (MQ=255)                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  1:4424055/50‑1 (MQ=255)                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      <  1:3797615/49‑1 (MQ=255)                                          aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    <  1:1363768/50‑1 (MQ=255)                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  1:3437088/50‑1 (MQ=255)                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:1857362/1‑50 (MQ=255)                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:4225982/1‑50 (MQ=255)                                            ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:2543254/1‑50 (MQ=255)                                               |                                               CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843366‑2843458

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