Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,850,651	A→G	H2R (CAC→CGC)	hypA →	hydrogenase 3 nickel incorporation protein HypA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,850,651	0	A	G	100.0%	59.2 / NA	18	H2R (CAC→CGC)	hypA	hydrogenase 3 nickel incorporation protein HypA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/9);  total (9/9)

AAAGATGAATTTCGTCGCCGTGTCGACGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCACGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTC  >  NC_000913/2850525‑2850774                                                                                                                               |                                                                                                                            aaaGATGAATTTCGTCGCCGTGTCGACGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAAt                                                                                                                       >  1:419546/1‑133 (MQ=255)               tcgCCGTGTCGACGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAAc                                                                                                         >  1:75565/1‑133 (MQ=255)                  ccGTGTCGACGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACggg                                                                                                      <  2:304371/133‑1 (MQ=255)                           cGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAAt                                                                                             <  1:442735/133‑1 (MQ=255)                                                            gCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCaa                                                            >  1:82063/1‑133 (MQ=255)                                                                ggATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAAcg                                                        <  2:82063/133‑1 (MQ=255)                                                                ggATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAAcg                                                        <  2:311007/133‑1 (MQ=255)                                                                  aTAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAAcgcg                                                      <  1:85672/133‑1 (MQ=255)                                                                           acTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACAGGGCGCAAAACGCGGAACTgggg                                             >  1:143420/1‑133 (MQ=255)                                                                            cTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGt                                            >  1:330382/1‑133 (MQ=255)                                                                               tcacaATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTc                                           <  1:417752/130‑1 (MQ=255)                                                                               gCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTg                                         >  1:64066/1‑133 (MQ=255)                                                                               gCACAATTATTGCTTGTAGCTGGCAATAGATAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTc                                           >  2:417752/1‑131 (MQ=255)                                                                                            ttGTAGCTAGCAATACTTCATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGgc                            <  1:243980/133‑1 (MQ=255)                                                                                               tAGCTGGCAATAGTTAATGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCa                         >  1:436008/1‑133 (MQ=255)                                                                                                      cAATAGTTGATGGGACGCGATACGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCtt                  <  1:188508/133‑1 (MQ=255)                                                                                                      cAATAGTTAATGGGATGCGATATGTGCGTAATAACCCTCTGACAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCtt                  <  2:143420/133‑1 (MQ=255)                                                                                                               aTGGGAGGCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAATTGGGGTCTCGCTCAAAATTGGCGCACTTTCATGTGTCGaaa         >  2:188508/1‑133 (MQ=255)                                                                                                                      gCGATATGCGCGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGctc  >  2:509486/1‑133 (MQ=255)                                                                                                                               |                                                                                                                            AAAGATGAATTTCGTCGCCGTGTCGACGTGTCATTTCGACATCATCGACATTATTCACCGCAGGGATAATCAACACTGGCACAATTATTGCTTGTAGCTGGCAATAGTTAATGGGAGGCGATATGCACGAAATAACCCTCTGCCAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTC  >  NC_000913/2850525‑2850774

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