Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,608,141	A→T	T409S (ACC→TCC)	hyfF →	hydrogenase 4, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,608,141	0	A	T	100.0%	90.0 / NA	29	T409S (ACC→TCC)	hyfF	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/15);  total (14/15)

GCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGC  >  NC_000913/2608017‑2608221                                                                                                                             |                                                                                 gcgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGa                                                                    <  2:279646/139‑1 (MQ=255)  cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGAc                                                                   >  1:1003542/1‑139 (MQ=255)                 gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAAt                                                        >  2:442867/1‑135 (MQ=255)                 gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAAt                                                        <  1:442867/135‑1 (MQ=255)                  tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCTCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACct                                                   <  1:805227/139‑1 (MQ=255)                  tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACct                                                   >  2:175702/1‑139 (MQ=255)                   cTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACctg                                                  <  1:175702/139‑1 (MQ=255)                            gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGc                                                               <  1:243292/117‑1 (MQ=255)                            gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGc                                                               >  2:243292/1‑117 (MQ=255)                            gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc                                         >  1:1049804/1‑139 (MQ=255)                            gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc                                         >  2:133039/1‑139 (MQ=255)                             cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCg                                        >  1:1084693/1‑139 (MQ=255)                                aaaaTCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGc                                                               >  1:673112/1‑113 (MQ=255)                                aaaaTCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGc                                                               <  2:673112/113‑1 (MQ=255)                                           aTTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTg                          <  2:1003542/139‑1 (MQ=255)                                                     gTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTcc                                     <  1:796924/118‑1 (MQ=255)                                                     gTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTcc                                     >  2:796924/1‑118 (MQ=255)                                                     gTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGc                <  2:16318/139‑1 (MQ=255)                                                      tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCctgc                                  <  1:458305/120‑1 (MQ=255)                                                      tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCctgc                                  >  2:458305/1‑120 (MQ=255)                                                      tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTg                          <  2:726215/128‑1 (MQ=255)                                                      tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTg                          >  1:726215/1‑128 (MQ=255)                                                                cggcggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGTCGGACTGGc                                                               >  2:1072565/1‑81 (MQ=255)                                                                cggcggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGTCGGACTGGc                                                               <  1:1072565/81‑1 (MQ=255)                                                                cggcggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTc                                      <  2:979205/106‑1 (MQ=255)                                                                cggcggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTc                                      >  1:979205/1‑106 (MQ=255)                                                                   cggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGc  <  1:133039/139‑1 (MQ=255)                                                                                                                 aGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCct                                    <  2:197287/59‑1 (MQ=255)                                                                                                                 aGCGAATTTATGTCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCct                                    >  1:197287/1‑59 (MQ=255)                                                                                                                             |                                                                                 GCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGC  >  NC_000913/2608017‑2608221

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: ‑