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

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458                                                |                                               ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              <  1:769248/50‑1 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             >  1:1882250/1‑50 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             >  1:5204045/1‑50 (MQ=255)  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             >  1:1509934/1‑50 (MQ=255)   tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:1720025/1‑50 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:1369429/1‑49 (MQ=255)    gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:4246947/1‑49 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:32118/1‑48 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:3439364/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  1:1168710/49‑1 (MQ=255)     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:25246/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         <  1:2897589/48‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         <  1:1271330/48‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:5287718/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:5714801/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  1:695061/1‑50 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:702589/50‑1 (MQ=255)        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        <  1:543090/49‑1 (MQ=255)           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  1:2559628/50‑1 (MQ=255)              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:4326464/49‑1 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:2983912/1‑49 (MQ=255)                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:107191/50‑1 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:3175941/50‑1 (MQ=255)                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:4728701/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           <  1:3117817/50‑1 (MQ=255)                    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           <  1:3674464/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:2637829/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  1:1009326/50‑1 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:2821765/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:246572/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  1:4818792/1‑50 (MQ=255)                     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  1:574748/50‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:3314402/1‑50 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:229171/50‑1 (MQ=255)                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCgtt                      <  1:5007807/50‑1 (MQ=255)                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:1964937/1‑49 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCGCCTTTTCa                   <  1:5428917/49‑1 (MQ=39)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:2910996/49‑1 (MQ=255)                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  <  1:364701/50‑1 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4590702/1‑50 (MQ=255)                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4029183/1‑50 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:2595851/1‑49 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:2116644/1‑49 (MQ=255)                               tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:1368097/1‑49 (MQ=255)                                ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  1:157195/49‑1 (MQ=255)                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:3290207/1‑50 (MQ=255)                                  catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             >  1:4201023/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:2049778/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:4125655/48‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:2314769/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:4302783/49‑1 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:2168999/1‑50 (MQ=255)                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:4940892/50‑1 (MQ=255)                                   atcaATAAAGTTATGACTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  1:5954375/50‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  1:186989/50‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  1:5000793/50‑1 (MQ=255)                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  1:2814924/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:4410764/50‑1 (MQ=255)                                       aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  1:1056974/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:550076/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:758669/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:287872/50‑1 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:4234080/1‑50 (MQ=255)                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       >  1:401956/1‑50 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  1:4629194/50‑1 (MQ=255)                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:1565650/1‑50 (MQ=255)                                          aaTCTATGGCTCAGCTGACTGAACACCACCTTTTCACTCATCGGATTcgc     >  1:4450205/1‑50 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    <  1:2041000/50‑1 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    <  1:400807/50‑1 (MQ=255)                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    <  1:2507178/50‑1 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:3992324/1‑50 (MQ=255)                                            tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  1:3670513/1‑50 (MQ=255)                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:4424386/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 < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑