Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,851,873	G→A	*291* (TAG→TAA)  G4S (GGC→AGC)	hypB → hypC →	GTP hydrolase involved in nickel liganding into hydrogenases hydrogenase maturation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,851,873	0	G	A	100.0%	158.1 / NA	49	*291* (TAG→TAA)  G4S (GGC→AGC)	hypB hypC	GTP hydrolase involved in nickel liganding into hydrogenases hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/24);  total (25/24)

TCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGG  >  NC_000913/2851737‑2852005                                                                                                                                         |                                                                                                                                     tcAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGc                                                                                                                                    >  2:782662/1‑139 (MQ=255)     cTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                                <  2:38915/139‑1 (MQ=255)        tGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCg                                                                                                                             <  2:384008/139‑1 (MQ=255)                  aaGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCg                                                                                                                   <  1:1505052/139‑1 (MQ=255)                               gcgcCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGc                                                                                                      <  1:782662/139‑1 (MQ=255)                                  cccGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCaa                                                                                                    <  1:1437164/138‑1 (MQ=255)                                  cccGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCaa                                                                                                    >  2:1437164/1‑138 (MQ=255)                                    cgcgAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCa                                                                                                 >  2:287105/1‑139 (MQ=255)                                      cgAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAgg                                                                                               >  1:574586/1‑139 (MQ=255)                                         aGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGa                                                                                            <  1:287105/139‑1 (MQ=255)                                                 cAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGAc                                                                                    <  2:535634/139‑1 (MQ=255)                                                     aaTTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCt                                                                                >  1:1389238/1‑139 (MQ=255)                                                         gAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGc                                                                                                                           >  1:7356/1‑92 (MQ=255)                                                         gAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGc                                                                                                                           <  2:7356/92‑1 (MQ=255)                                                                 ccTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAgc                                                                    >  1:103094/1‑139 (MQ=255)                                                                       ttCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATg                                                              <  2:103094/139‑1 (MQ=255)                                                                           gCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACg                                                                                                        >  2:1410664/1‑93 (MQ=255)                                                                           gCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACg                                                                                                        <  1:1410664/93‑1 (MQ=255)                                                                            ccaccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCg                                                                                                                   <  1:1497252/81‑1 (MQ=255)                                                                            ccaccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCg                                                                                                                   >  2:1497252/1‑81 (MQ=255)                                                                             caccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAtt                                                        >  1:1653856/1‑139 (MQ=255)                                                                                caGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                                >  2:199024/1‑64 (MQ=255)                                                                                caGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                                <  1:199024/64‑1 (MQ=255)                                                                                             aTGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAg                                        <  2:1389238/139‑1 (MQ=255)                                                                                              tGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      >  1:52771/1‑92 (MQ=255)                                                                                              tGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      <  2:52771/92‑1 (MQ=255)                                                                                              tGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGc                                       >  1:428817/1‑139 (MQ=255)                                                                                              tGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGc                                       >  1:928711/1‑139 (MQ=255)                                                                                                      tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  <  1:1275558/88‑1 (MQ=255)                                                                                                      tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  >  2:1275558/1‑88 (MQ=255)                                                                                                      tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGa                               >  1:1465958/1‑139 (MQ=255)                                                                                                       ggCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGaa                              >  1:300330/1‑139 (MQ=255)                                                                                                        gCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGaaa                             <  2:574586/139‑1 (MQ=255)                                                                                                                      agaCACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAt                                                         <  1:748563/97‑1 (MQ=255)                                                                                                                      agaCACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAt                                                         >  2:748563/1‑97 (MQ=255)                                                                                                                       gaCACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGa                                                          <  2:1279472/95‑1 (MQ=255)                                                                                                                       gaCACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGa                                                          >  1:1279472/1‑95 (MQ=255)                                                                                                                       gaCACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCcgcgcg              <  2:428817/139‑1 (MQ=255)                                                                                                                            aGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCa                                                                                                                 >  2:1696476/1‑35 (MQ=255)                                                                                                                            aGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCa                                                                                                                 <  1:1696476/35‑1 (MQ=255)                                                                                                                                aTGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  >  1:1314118/1‑62 (MQ=255)                                                                                                                                aTGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  <  2:1314118/62‑1 (MQ=255)                                                                                                                                 tgtgCATAAGCGTTCCCGGCCAGATCCTCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGt                                                  <  2:224710/93‑1 (MQ=255)                                                                                                                                 tgtgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGt                                                  >  1:224710/1‑93 (MQ=255)                                                                                                                                  gtgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  >  2:1670551/1‑60 (MQ=255)                                                                                                                                  gtgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                  <  1:1670551/60‑1 (MQ=255)                                                                                                                                   tgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTggg  >  1:629970/1‑139 (MQ=255)                                                                                                                                    gCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      <  1:408760/54‑1 (MQ=255)                                                                                                                                    gCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      >  2:408760/1‑54 (MQ=255)                                                                                                                                         |                                                                                                                                     TCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGG  >  NC_000913/2851737‑2852005

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