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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458                                                |                                               ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  1:3702156/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  1:4475515/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  2:2664650/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  1:2522146/50‑1 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGGATCGTTCTCATCAATAAATTTATGGCt                                            >  2:4097047/1‑50 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:3661284/1‑50 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  2:1025834/1‑50 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  2:1154738/1‑49 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:4458581/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  2:516333/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  1:38216/48‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  2:340482/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:3823333/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:3712976/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  2:3000601/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         <  2:1175949/48‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  2:2139702/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:2316622/50‑1 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  2:4032270/1‑49 (MQ=255)         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCCCAGtt                                       >  2:3070049/1‑49 (MQ=255)          ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                                     <  2:1506876/50‑1 (MQ=255)          ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                                     <  1:4305641/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  1:1169216/50‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    >  1:218718/1‑50 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:1318989/48‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  >  1:1953268/1‑49 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:3504827/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:3321347/50‑1 (MQ=255)               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:4406422/1‑50 (MQ=255)               gcgggcgAGGCATCGTTCTCATCAATAAATTTATGCCTCAGTTGACTGaa                                >  2:4296320/1‑50 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:3681031/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:1968750/50‑1 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:4089534/50‑1 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  2:22768/1‑50 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:682733/50‑1 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  1:2489453/50‑1 (MQ=255)                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  1:469925/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           >  1:2667918/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:4255574/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  2:3143307/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:3284061/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  2:2859203/1‑50 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  1:1201566/1‑50 (MQ=255)                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  2:3749824/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  2:762276/50‑1 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  2:4175217/1‑50 (MQ=255)                        gCATCGTTCTCATCAATAAATTTATGACTCAGTTGACTGAACACCACCt                        <  2:3759745/49‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:2366898/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  2:497110/50‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:2134644/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:81225/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:2358907/1‑50 (MQ=255)                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  1:4180194/50‑1 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                    <  1:705732/50‑1 (MQ=255)                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  2:244713/1‑50 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:552396/49‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  <  2:3231406/50‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  2:3465380/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  2:2872207/50‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  2:856556/49‑1 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:3010804/49‑1 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:1853011/1‑50 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  2:2867971/1‑50 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:1611757/1‑50 (MQ=255)                                  aatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               <  2:3804411/47‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  2:2698188/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:3476237/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:4027479/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:1930260/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:1578538/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCATTTCACGCag              <  2:2910553/48‑2 (MQ=37)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:2443774/1‑50 (MQ=255)                                     caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  1:2112516/1‑50 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          <  1:1593782/49‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:899841/1‑49 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  2:1035190/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        >  2:340300/1‑50 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:2182014/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:816549/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:3725549/1‑50 (MQ=255)                                        caaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:669373/2‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACGCCACCTTTTCACTCATCGGATAcgc     <  2:2461679/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  2:3718984/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  1:3036327/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  2:4191216/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  2:282817/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      <  2:381838/49‑1 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    >  1:764232/1‑50 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  2:3261323/50‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  2:2696845/50‑1 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  2:1802010/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  2:986619/1‑50 (MQ=255)                                                |                                               CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458

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