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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458                                                |                                               ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  2:230475/50‑1 (MQ=255) ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  2:4128070/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  1:5671853/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  2:2380751/50‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  2:3132173/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:4958656/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:324025/50‑1 (MQ=255)      tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         >  1:897289/1‑50 (MQ=255)       gcgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  2:4971574/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:1147237/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  2:3274399/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:2386780/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:104083/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:5644082/49‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:2559366/49‑1 (MQ=255)        cgcCTCGGCGGGCGTAGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:2975935/49‑1 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  2:2820038/1‑49 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  2:3310890/50‑1 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:3770314/50‑1 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  2:32029/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  2:1950706/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  1:4072058/50‑1 (MQ=255)            tCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  2:3178502/1‑50 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  >  2:3659703/1‑49 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:4647913/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:4324352/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  1:5301646/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  2:1619858/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:4920552/1‑50 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  1:3666563/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  1:3527472/50‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  2:3184450/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  2:5694865/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  2:4911798/1‑50 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  2:4812830/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  2:4416380/1‑50 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:742010/50‑1 (MQ=255)                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  1:1609317/1‑50 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  2:1511929/1‑50 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:2369367/1‑50 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  2:5235006/50‑1 (MQ=255)                  ggcgTGGCATCGTTCGCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  2:4142768/50‑1 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  2:500243/50‑1 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  2:2647288/50‑1 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  2:5531488/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  2:3456093/1‑50 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  1:3215776/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  2:608294/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:4940010/1‑50 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  2:5159745/50‑1 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  2:2600219/1‑50 (MQ=255)                       ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        >  2:4718352/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:1981481/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  2:2543610/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  2:295017/49‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  2:1919582/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:1460449/50‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  2:1179124/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  2:598870/1‑50 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCTCCTTTTc                    <  2:1360125/50‑1 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  2:5510361/1‑49 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                    <  1:3946277/50‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  2:1452927/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  2:4506414/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:744022/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  2:3705374/50‑1 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4220896/1‑49 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:1637274/1‑49 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  2:1694947/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:3577856/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               >  1:171995/1‑50 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  2:2016256/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               <  1:4603477/47‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:4164292/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:5170841/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:5035540/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:3616354/50‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  2:4121537/50‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:4648550/50‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:1327300/50‑1 (MQ=255)                                    tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:609685/50‑1 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:4106417/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:518712/1‑49 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:3037400/1‑49 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  1:676285/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:1536729/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  2:4629760/1‑50 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:2003976/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  2:1198230/1‑50 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  2:1346329/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  2:838681/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  1:4155841/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  2:3188615/50‑1 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    >  2:2362631/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:5393776/1‑50 (MQ=255)                                                |                                               CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑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: ‑