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%	127.3 / NA	40	*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 (23/17);  total (23/17)

TGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTG  >  NC_000913/2851744‑2852003                                                                                                                                  |                                                                                                                                   tGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCg                                                                                                                           >  1:1113243/1‑139 (MQ=255) tGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCg                                                                                                                           >  2:571800/1‑139 (MQ=255)  gACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCgg                                                                                                                          >  2:711132/1‑139 (MQ=255)        gAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGAt                                                                                                                    <  2:46715/139‑1 (MQ=255)             gTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCa                                                                                                               >  2:516729/1‑139 (MQ=255)                               cgAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAgg                                                                                             >  1:948779/1‑139 (MQ=255)                                                        atcCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCa                                                                    >  2:175303/1‑139 (MQ=255)                                                           catATGTCCGCCGCCAGCGGCGAAGGGATGGACCAGTGGCTAAACTGGCTGGAGACACAGCGATGTGTATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGACGAAAGTCGACGTCTGCGGCATTCAgcg                                                                 <  2:836414/137‑1 (MQ=255)                                                                ttCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATg                                                            <  1:516729/139‑1 (MQ=255)                                                                   cGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGcc                                                                                                                        <  2:694217/76‑1 (MQ=255)                                                                   cGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGcc                                                                                                                        >  1:694217/1‑76 (MQ=255)                                                                    gCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                              <  2:781051/69‑1 (MQ=255)                                                                    gCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTc                                                                                                                              >  1:781051/1‑69 (MQ=255)                                                                    gCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGa                                                        >  2:132415/1‑139 (MQ=255)                                                                        ccaGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTa                                                    >  2:1213066/1‑139 (MQ=255)                                                                            cggcgAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGt                                                <  2:1113243/139‑1 (MQ=255)                                                                               cgAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGa                                                                                          <  1:1101472/94‑1 (MQ=255)                                                                               cgAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCAGTGACGGCAACCAGGCGa                                                                                          >  2:1101472/1‑94 (MQ=255)                                                                                 aaGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGACAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTCGTc                                           >  2:274895/1‑139 (MQ=255)                                                                                           ccAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTc                                                                                     >  1:32166/1‑87 (MQ=255)                                                                                           ccAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTc                                                                                     <  2:32166/87‑1 (MQ=255)                                                                                               tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCa                                                                                                               <  1:48771/57‑1 (MQ=255)                                                                                               tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCa                                                                                                               >  2:48771/1‑57 (MQ=255)                                                                                               tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGa                             >  2:1165050/1‑139 (MQ=255)                                                                                               tGGCTGAACTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGa                             >  2:1342791/1‑139 (MQ=255)                                                                                                      aCTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAgc                                                                  <  1:1504474/95‑1 (MQ=255)                                                                                                      aCTGGCTGGAGACACAGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAgc                                                                  >  2:1504474/1‑95 (MQ=255)                                                                                                         ggctggAGACACAGCGATGTGCGTAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAg                   >  2:566606/1‑139 (MQ=255)                                                                                                                    caGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCg                                                                           >  2:292562/1‑72 (MQ=255)                                                                                                                    caGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCg                                                                           <  1:292562/72‑1 (MQ=255)                                                                                                                    caGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAtt                                                      <  1:443824/93‑1 (MQ=255)                                                                                                                    caGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAtt                                                      >  2:443824/1‑93 (MQ=255)                                                                                                                    caGCGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTggg        <  1:571800/139‑1 (MQ=255)                                                                                                                       cGATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCa     <  1:175303/139‑1 (MQ=255)                                                                                                                        gATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGa                                                        >  2:1282779/1‑87 (MQ=255)                                                                                                                        gATGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGa                                                        <  1:1282779/87‑1 (MQ=255)                                                                                                                         aTGTGCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGt   >  2:604513/1‑139 (MQ=255)                                                                                                                          tgtgCATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTg  <  1:604513/139‑1 (MQ=255)                                                                                                                              cATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCAtt                                                                      <  1:1044928/67‑1 (MQ=255)                                                                                                                              cATAAGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCAtt                                                                      >  2:1044928/1‑67 (MQ=255)                                                                                                                                  |                                                                                                                                   TGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTG  >  NC_000913/2851744‑2852003

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