Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,602,223	(T)5→4	coding (406/2019 nt)	hyfB →	hydrogenase 4 component B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,602,219	0	T	.	100.0%	115.4 / NA	28	coding (402/2019 nt)	hyfB	hydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/12);  total (16/12)

GCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAAC  >  NC_000913/2602140‑2602293                                                                                |                                                                           gCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCtt                                                                           >  1:464362/1‑81 (MQ=255) gCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCtt                                                                           >  1:464348/1‑81 (MQ=255)   aaaGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCtttt                                                                         >  2:464361/1‑81 (MQ=255)     aGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCttttgg                                                                       >  1:464349/1‑79 (MQ=255)        ggcggcggcgATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTc                                                                   <  1:464316/80‑1 (MQ=255)          cggcggcgATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCAt                                                                 <  1:464317/81‑1 (MQ=255)              ggcgATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTg                                                             >  2:464376/1‑81 (MQ=255)                   tGGGCTCCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGtt                                                        >  2:464356/1‑81 (MQ=255)                     ggCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTg                                                      <  1:464318/81‑1 (MQ=255)                         tcttTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAat                                                  >  1:464351/1‑81 (MQ=255)                             tATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAatgatg                                              >  1:464371/1‑81 (MQ=255)                               tGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTc                                            <  1:464319/81‑1 (MQ=255)                               tGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTc                                            >  2:464341/1‑81 (MQ=255)                                      tttCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCt                                     >  1:464350/1‑81 (MQ=255)                                        tCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTc                                   >  1:464373/1‑81 (MQ=255)                                         catcgcatcgATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTc                                   <  1:464320/80‑1 (MQ=255)                                          atcgcatcgATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCt                                 <  2:464321/81‑1 (MQ=255)                                             gcatcgATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGt                              >  1:464354/1‑81 (MQ=255)                                             gcatcgATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGt                              >  1:464355/1‑81 (MQ=255)                                                tcgATGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTc                           <  2:464322/81‑1 (MQ=255)                                                   aTGGTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTCGAAATGATGTCGCTGTCTTCCTGGTTTCTgg                        >  2:464353/1‑81 (MQ=255)                                                     ggTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGt                       >  2:464379/1‑80 (MQ=255)                                                     ggTTGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGt                       >  1:464365/1‑80 (MQ=255)                                                       ttGCCCTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCAt                    <  1:464323/81‑1 (MQ=255)                                                           ccTGCTGGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGcc                <  2:464324/81‑1 (MQ=255)                                                             tgctgGTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCcag              >  2:464370/1‑81 (MQ=255)                                                                 ggTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCcaggcag          >  2:464384/1‑81 (MQ=255)                                                                 ggTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCcaggcag          >  1:464366/1‑81 (MQ=255)                                                                  gTGATGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCcaggca           <  1:464325/79‑1 (MQ=255)                                                                     aTGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATa      <  2:464326/81‑1 (MQ=255)                                                                      tGGACAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATa      >  2:464381/1‑80 (MQ=255)                                                                         aCAACGC‑TTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAAc  <  1:464327/81‑1 (MQ=255)                                                                                |                                                                           GCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAAC  >  NC_000913/2602140‑2602293

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