Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,602,161	Δ1 bp	coding (344/2019 nt)	hyfB →	hydrogenase 4, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,602,159	0	G	.	100.0%	114.4 / NA	27	coding (342/2019 nt)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/12);  total (15/12)

ATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGT  >  NC_000913/2602024‑2602272                                                                                                                                        |                                                                                                                  aTTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggct                                                                                                                <  2:205606/139‑3 (MQ=255)   ttGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggcttc                                                                                                              >  2:328620/1‑135 (MQ=255)   ttGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggcttc                                                                                                              >  1:20117/1‑135 (MQ=255)          tttAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGa                                                                                                      <  2:20117/139‑1 (MQ=255)                    cGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTcat                                                                                            >  1:150049/1‑139 (MQ=255)                     gAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggct                                                                                                                >  2:7518/1‑117 (MQ=255)                     gAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggct                                                                                                                <  1:7518/119‑3 (MQ=255)                     gAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTcatc                                                                                           <  1:325769/139‑1 (MQ=255)                               ccTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATgg                                                                                 <  1:109436/139‑1 (MQ=255)                                    cGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGcc                                                                            <  1:328620/139‑1 (MQ=255)                                                  cTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGa                                                              >  1:400666/1‑139 (MQ=255)                                                  cTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGa                                                              >  1:96140/1‑139 (MQ=255)                                                    ggTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACa                                                            <  2:96140/139‑1 (MQ=255)                                                               tGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTcatcgcatc                                                                                      >  1:664240/1‑102 (MQ=255)                                                               tGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTcatcgcatc                                                                                      <  2:664240/102‑1 (MQ=255)                                                               tGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTgg                                                 >  1:558102/1‑139 (MQ=255)                                                                      ggtggtTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCg                                          >  1:135269/1‑139 (MQ=255)                                                                        tggtTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTg                                        <  2:150049/139‑1 (MQ=255)                                                                         ggtTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTgg                                                                     <  1:742509/109‑1 (MQ=255)                                                                         ggtTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTgg                                                                     >  2:742509/1‑109 (MQ=255)                                                                                              gACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTg                                        <  1:80806/117‑1 (MQ=255)                                                                                              gACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTg                                        >  2:80806/1‑117 (MQ=255)                                                                                              gACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGt                  >  1:182599/1‑139 (MQ=255)                                                                                              gACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGt                  >  1:671676/1‑139 (MQ=255)                                                                                              gACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGt                  >  1:708446/1‑139 (MQ=255)                                                                                                    tatGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCt            >  2:564788/1‑139 (MQ=255)                                                                                                         gcgAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAAc                                                          <  2:649543/88‑1 (MQ=255)                                                                                                         gcgAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAAc                                                          >  1:649543/1‑88 (MQ=255)                                                                                                         gcgAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGttt       >  1:148988/1‑139 (MQ=255)                                                                                                           gAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGAt                                                                                   >  2:147906/1‑61 (MQ=255)                                                                                                           gAATACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGAt                                                                                   <  1:147906/61‑1 (MQ=255)                                                                                                              tACGAGGGCAAAGGCGCGGCGGCGAT‑GGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGt  <  2:671676/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                  ATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGT  >  NC_000913/2602024‑2602272

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