Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,602,017	(C)7→8	coding (200/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,010	1	.	C	100.0%	72.5 / NA	23	A65P (GCC→CCC)	hyfB	hydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/12);  total (11/12)

ATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCC  >  NC_000913/2601931‑2602084                                                                                 |                                                                           aTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑c                                                                           <  1:464288/81‑1 (MQ=255) aTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑c                                                                           >  2:464313/1‑81 (MQ=255)    gCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑cccc                                                                        >  1:464338/1‑81 (MQ=255)        gtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑cccccc                                                                      <  1:464289/79‑1 (MQ=255)              ggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGa                                                               <  2:464290/81‑1 (MQ=255)                    aTTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTg                                                          <  2:464291/80‑1 (MQ=255)                     ttGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGcc                                                        <  1:464292/81‑1 (MQ=255)                       ggCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGa                                                      >  1:464321/1‑81 (MQ=255)                        gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGAt                                                     >  2:464316/1‑81 (MQ=255)                                 gCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGCTTTAACCCTg                                            >  1:464328/1‑81 (MQ=255)                                   ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                          <  1:464294/81‑1 (MQ=255)                                   ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                          <  2:464293/81‑1 (MQ=255)                                         aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATgg                                     <  2:464295/80‑1 (MQ=255)                                          cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                   <  1:464296/81‑1 (MQ=255)                                          cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                   <  2:464297/81‑1 (MQ=255)                                            cAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAg                                 <  1:464298/81‑1 (MQ=255)                                            cAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAg                                 >  2:464325/1‑81 (MQ=255)                                              aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCGGCGAATGGATAGcc                               <  2:464299/81‑1 (MQ=255)                                               ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCct                              <  1:464300/81‑1 (MQ=255)                                               ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCct                              >  2:464329/1‑81 (MQ=255)                                                             aaCTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATgg                >  2:464318/1‑81 (MQ=255)                                                                     tCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGtt        <  1:464301/81‑1 (MQ=255)                                                                      cGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa       >  1:464326/1‑81 (MQ=255)                                                                      cGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa       >  1:464322/1‑81 (MQ=255)                                                                      cGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa       >  2:464320/1‑81 (MQ=255)                                                                         ccGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAtct    >  1:464335/1‑81 (MQ=255)                                                                           gTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTcc  >  2:464319/1‑81 (MQ=255)                                                                                 |                                                                           ATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCC  >  NC_000913/2601931‑2602084

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