Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,611,639	T→C	D61D (GAT→GAC)	hyfJ →	putative hydrogenase 4 assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,611,639	0	T	C	100.0%	192.9 / NA	52	D61D (GAT→GAC)	hyfJ	putative hydrogenase 4 assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/23);  total (29/23)

GTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACG  >  NC_000913/2611508‑2611774                                                                                                                                    |                                                                                                                                        gTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGc                                                                                                                                 <  1:1076583/140‑1 (MQ=255)  tCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCa                                                                                                                                <  1:1006771/140‑1 (MQ=255)  tCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCa                                                                                                                                >  1:577750/1‑140 (MQ=255)      cAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTgg                                                                                                                            <  1:979029/140‑1 (MQ=255)       aGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGt                                                                                                                           <  2:465383/140‑1 (MQ=255)                               aGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTaa                                                                                                   <  2:577750/140‑1 (MQ=255)                               aGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAACAAGCCCGACGTaa                                                                                                   <  2:145342/140‑1 (MQ=255)                                gCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAg                                                                                                  <  1:286933/140‑1 (MQ=255)                                 ctctCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGa                                                                                                 >  1:277958/1‑140 (MQ=255)                                 ctctCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGa                                                                                                 >  1:308732/1‑140 (MQ=255)                                     cAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCt                                                                                             <  2:308732/140‑1 (MQ=255)                                      aGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTg                                                                                            >  1:1158930/1‑140 (MQ=255)                                      aGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTg                                                                                            >  2:805791/1‑140 (MQ=255)                                                ttACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGtgatga                                                                                  >  1:957856/1‑140 (MQ=255)                                                  aCTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACt                                                                                <  1:215279/140‑1 (MQ=255)                                                  aCTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACt                                                                                >  1:591818/1‑140 (MQ=255)                                                  aCTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGGGATGACt                                                                                >  2:1047138/1‑140 (MQ=255)                                                   cTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACtt                                                                               >  2:605736/1‑140 (MQ=255)                                                   cTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACtt                                                                               >  2:746553/1‑140 (MQ=255)                                                      gctagctaTCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTg                                                                            >  2:74320/1‑140 (MQ=255)                                                       ctagctaTCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTgg                                                                           <  1:605736/140‑1 (MQ=255)                                                        tagctaTCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGt                                                                          <  2:957856/140‑1 (MQ=255)                                                             aTCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGAt                                                                     >  2:1153899/1‑140 (MQ=255)                                                                  ccATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTa                                                                >  1:808281/1‑140 (MQ=255)                                                                   cATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTAt                                                               <  2:808281/140‑1 (MQ=255)                                                                    aTCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTAtt                                                              >  2:555524/1‑140 (MQ=255)                                                                      cACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGa                                                            >  1:1189868/1‑140 (MQ=255)                                                                       aCGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGAc                                                           <  1:74320/140‑1 (MQ=255)                                                                            ggCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATACTGGGGGTGATGACTTTTGGTGAGATTGTTATTGGCgccag                                                      >  1:1086537/1‑140 (MQ=255)                                                                                tgattgatTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCcac                                                  >  1:214166/1‑140 (MQ=255)                                                                                tgattgatTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCcac                                                  >  2:989735/1‑140 (MQ=255)                                                                                tgattgatTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCcac                                                  >  1:153033/1‑140 (MQ=255)                                                                                   ttgattgTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACAcc                                               >  1:133936/1‑140 (MQ=255)                                                                                    tgattgTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATACTGGGGGTGATGACTTTTGGTGAGATTGTTATTGGCGCCAGCCACACCg                                              <  2:1086537/140‑1 (MQ=255)                                                                                              aaTGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGAc                                    <  1:1153899/140‑1 (MQ=255)                                                                                                    gCCTTCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGc                              >  1:863152/1‑140 (MQ=255)                                                                                                       ttCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCAtt                           <  2:1189868/140‑1 (MQ=255)                                                                                                        tCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTc                          >  1:787803/1‑140 (MQ=255)                                                                                                        tCCGCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTc                          >  2:991434/1‑140 (MQ=255)                                                                                                           gCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTcgcg                       <  1:555524/140‑1 (MQ=255)                                                                                                           gCTGCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTcgcg                       <  2:1158930/140‑1 (MQ=255)                                                                                                              gCCCACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCa                    <  1:574046/140‑1 (MQ=255)                                                                                                               cccACTGACGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCAc                   >  2:153963/1‑140 (MQ=255)                                                                                                                      aCGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGa            <  1:879386/140‑1 (MQ=255)                                                                                                                       cGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTcgcg                       >  2:737080/1‑128 (MQ=255)                                                                                                                       cGGAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTcgcg                       <  1:737080/128‑1 (MQ=255)                                                                                                                         gAATACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGa         <  1:805791/140‑1 (MQ=255)                                                                                                                            tACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGAcgc      >  1:512188/1‑140 (MQ=255)                                                                                                                            tACGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGAcgc      <  1:989735/140‑1 (MQ=255)                                                                                                                              cGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGa    >  1:532362/1‑140 (MQ=255)                                                                                                                              cGAAGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGa    >  1:1001189/1‑140 (MQ=255)                                                                                                                                aaGACTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACg  >  1:364734/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                        GTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACG  >  NC_000913/2611508‑2611774

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