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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458                                                |                                               ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              >  1:8171520/1‑50 (MQ=255) ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  1:7649345/50‑1 (MQ=255) ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  1:6089557/50‑1 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:6822332/1‑49 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:7532102/1‑49 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:2441598/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:5855657/1‑48 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:2622500/49‑1 (MQ=255)       gcgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  1:889368/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         <  1:7361077/48‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:362066/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:3699974/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:7328648/1‑50 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCCCAGtt                                       >  1:1448913/1‑49 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAAATTATGGCTCCGTTGAc                                    >  1:7510937/1‑50 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:5699346/48‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:4414263/49‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  1:83408/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:2595491/1‑50 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:1353005/1‑50 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:8265003/1‑50 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  1:5420763/50‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:7053515/50‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:2518826/50‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:1761493/50‑1 (MQ=255)                 gggctTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:4065200/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:5156942/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:4520554/1‑50 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:169572/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:905397/1‑50 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:7178447/1‑50 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:5409590/1‑49 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  1:6513417/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           <  1:976377/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  1:3576451/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  1:1879783/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:8178083/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:2980478/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  1:2795976/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:973125/49‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:430172/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:1511556/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:7274209/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:5259287/50‑1 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:4140351/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:6019491/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  1:8218368/50‑1 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  1:1046273/50‑1 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:8027330/1‑50 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                    <  1:6967282/50‑1 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:2052659/1‑50 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  1:6058539/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:2848819/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:7281331/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:4694537/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:6990849/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4796601/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:6317761/1‑50 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:1303201/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:3928421/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:2314866/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               >  1:1925807/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:536045/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:4726667/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:3739533/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:7471672/49‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:2556036/50‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:2437206/1‑49 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:5239161/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:7564094/1‑49 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:2420112/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:113323/1‑50 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:5207419/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:306267/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:3363396/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:3879248/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      <  1:149025/49‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:3287439/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:4881149/1‑50 (MQ=255)                                                |                                               CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458

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