Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,611,715	G→A	D87N (GAC→AAC)	hyfJ →	putative processing element hydrogenase 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,611,715	0	G	A	100.0%	107.2 / NA	34	D87N (GAC→AAC)	hyfJ	putative processing element hydrogenase 4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/15);  total (19/15)

ATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCATCT  >  NC_000913/2611589‑2611844                                                                                                                               |                                                                                                                                  attgattgTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCcaca                                                                                                                       >  1:292854/1‑139 (MQ=255)            gAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCcaca                                                                                                                       <  1:223518/128‑1 (MQ=255)            gAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCcaca                                                                                                                       >  2:223518/1‑128 (MQ=255)                       tCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCAtt                                                                                                 >  2:456648/1‑139 (MQ=255)                       tCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCAtt                                                                                                 >  1:822877/1‑139 (MQ=255)                        ccacTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCg                                                                                                                    >  2:399457/4‑119 (MQ=255)                        ccacTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCg                                                                                                                    <  1:399457/116‑1 (MQ=255)                                       ggAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGAt                                                                                 <  2:831790/139‑1 (MQ=255)                                        gAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCcaca                                                                                                                       <  1:83301/100‑1 (MQ=255)                                        gAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCcaca                                                                                                                       >  2:83301/1‑100 (MQ=255)                                           tACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCa                                                                                          >  1:233948/1‑126 (MQ=255)                                           tACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCa                                                                                          <  2:233948/126‑1 (MQ=255)                                                     ggCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCtgt                                                                   >  2:200863/1‑139 (MQ=255)                                                         tGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTgg                                                               <  1:456648/139‑1 (MQ=255)                                                           cACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCa                                                             >  2:296612/1‑139 (MQ=255)                                                                gTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCg                                                                                                                    <  1:655919/79‑1 (MQ=255)                                                                gTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCg                                                                                                                    >  2:655919/1‑79 (MQ=255)                                                                gTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGgcgc                                                        >  2:289830/1‑139 (MQ=255)                                                                          agCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGAAGACGCGACGTCTGTGTGGCAGGCGCGTAGCattca                                              >  1:562013/1‑139 (MQ=255)                                                                            cAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCat                                                 >  1:610289/1‑134 (MQ=255)                                                                            cAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCat                                                 <  2:610289/134‑1 (MQ=255)                                                                             aaGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCattcaatt                                           <  1:289830/139‑1 (MQ=255)                                                                                 ccGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAattcatt                                       >  1:20206/1‑139 (MQ=255)                                                                                              cTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTaaca                                                  <  2:380796/115‑4 (MQ=255)                                                                                              cTGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAACa                                                  >  1:380796/1‑115 (MQ=255)                                                                                               tGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGTCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGAcg                                                                             <  2:110799/87‑1 (MQ=255)                                                                                               tGGGGGTGATGACTTTTGGTGAGATTGTTATTAACGCCAGTCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGAcg                                                                             >  1:110799/1‑87 (MQ=255)                                                                                                    gTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCattcaattca                                         <  2:160928/118‑1 (MQ=255)                                                                                                    gTGATGACTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCattcaattca                                         >  1:160928/1‑118 (MQ=255)                                                                                                       atgaCTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGc                                                                                     <  2:38690/71‑1 (MQ=255)                                                                                                       atgaCTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGc                                                                                     >  1:38690/1‑71 (MQ=255)                                                                                                         gaCTTTTGGTGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAg               >  1:497714/1‑139 (MQ=255)                                                                                                                  tGAGATTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGcc      <  2:497714/139‑1 (MQ=255)                                                                                                                      aTTGTTATTAACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCatct  <  1:296612/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                  ATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAGCAAGCCCGACGTAAGATGCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCATCT  >  NC_000913/2611589‑2611844

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