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%	62.4 / NA	21	*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 (12/9);  total (12/9)

CATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGG  >  NC_000913/2851759‑2852005                                                                                                                   |                                                                                                                                     cATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCa                                                                                                              >  1:958592/1‑139 (MQ=255)      ccTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGAc                                                                                                         >  2:138470/1‑139 (MQ=255)                cgAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGc                                                                                                                  >  1:712551/1‑120 (MQ=255)                cgAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGc                                                                                                                  <  2:712551/120‑1 (MQ=255)                                           ccTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGGACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAgc                                                                    >  2:191582/1‑139 (MQ=255)                                                    cGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCg                                                           >  1:145449/1‑139 (MQ=255)                                                     gCCACCAGCGGCGAAGGGATGGACTAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                                >  1:790414/1‑69 (MQ=255)                                                     gCCACCAGCGGCGAAGGGATGGACAAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                                <  2:790414/69‑1 (MQ=255)                                                      ccaccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCaa                                                                                                    >  1:696190/1‑96 (MQ=255)                                                      ccaccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCaa                                                                                                    <  2:696190/96‑1 (MQ=255)                                                                       aTGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAg                                        <  2:958592/139‑1 (MQ=255)                                                                         ggACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCt                                      <  2:22831/139‑1 (MQ=255)                                                                               gTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATg                                >  1:246669/1‑139 (MQ=255)                                                                               gTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATg                                >  1:290899/1‑139 (MQ=255)                                                                                tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGa                               >  2:976891/1‑139 (MQ=255)                                                                                          ggctggAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAg                     >  1:980542/1‑139 (MQ=255)                                                                                                      aGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGc         <  2:778004/139‑1 (MQ=255)                                                                                                          acaggCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCATGCGAACGTCCACGTCTGCGGCATTCAGCTCGATGTCGATTTAACGTTAGTCCGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGt     <  2:532097/135‑1 (MQ=255)                                                                                                          aTGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      >  2:832213/1‑58 (MQ=255)                                                                                                          aTGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                      <  1:832213/58‑1 (MQ=255)                                                                                                             tgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTggg  <  2:145449/139‑1 (MQ=255)                                                                                                                   |                                                                                                                                     CATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGG  >  NC_000913/2851759‑2852005

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