Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,034,326	A→G	N358S (AAC→AGC)	hyaB →	hydrogenase 1, large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,034,326	0	A	G	100.0%	174.2 / NA	49	N358S (AAC→AGC)	hyaB	hydrogenase 1, large subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/20);  total (29/20)

AACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCTG  >  NC_000913/1034235‑1034420                                                                                            |                                                                                               aaCGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAgg                                                                                        >  1:1171812/1‑100 (MQ=255) aaCGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAgg                                                                                        >  1:690091/1‑100 (MQ=255)      cGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCggg                                                                                   >  1:180475/1‑100 (MQ=255)      cGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCggg                                                                                   >  2:1240102/1‑100 (MQ=255)         nttCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGt                                                                                >  2:868392/2‑100 (MQ=255)         cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGt                                                                                <  1:243629/100‑1 (MQ=255)         cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGt                                                                                >  2:930069/1‑100 (MQ=255)          ttCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTc                                                                               <  2:598799/100‑1 (MQ=255)           ncaacaaTGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCa                                                                              >  2:492763/2‑100 (MQ=255)            caacaaTGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCa                                                                              >  1:417615/1‑99 (MQ=255)            caacaaTGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCAt                                                                             >  2:155667/1‑100 (MQ=255)            caacaaTGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCAt                                                                             >  1:1251613/1‑100 (MQ=255)               ccaggggctgaCAGTGTATTTGGTCGTTCCGGAGCAGGTGTATGAGTTTGTCGACCGCGCCCGGTATCGCTATCCCAGCGATCAGGTCGGGCGTCATc                                                                            <  2:599981/92‑1 (MQ=255)                   gtgCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCg                                                                      <  2:1091812/100‑1 (MQ=255)                     gctgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGAt                                                                    >  2:682287/1‑100 (MQ=255)                      ctgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATg                                                                   >  1:587691/1‑100 (MQ=255)                      ctgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATg                                                                   <  1:184028/100‑1 (MQ=255)                        gcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGc                                                                 >  1:216715/1‑100 (MQ=255)                          cAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTATGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCAt                                                               <  1:1240113/100‑1 (MQ=255)                               gATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCAcc                                                           <  1:845814/99‑1 (MQ=255)                                aTTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGa                                                         >  1:890077/1‑100 (MQ=255)                                  ttGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGAcc                                                       <  2:207211/100‑1 (MQ=255)                                   tGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGAccc                                                      >  1:952634/1‑100 (MQ=255)                                      ttGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTg                                                   >  1:69338/1‑100 (MQ=255)                                         aTCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTa                                                >  1:969651/1‑100 (MQ=255)                                         aTCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTa                                                >  1:107400/1‑100 (MQ=255)                                           ccGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACa                                              >  2:660502/1‑100 (MQ=255)                                             gcagcaGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAAc                                            <  1:593968/100‑1 (MQ=255)                                                  agGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCgg                                       <  2:726210/100‑1 (MQ=255)                                                   ggTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGc                                      <  1:751161/100‑1 (MQ=255)                                                    gTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCccgccg                                     >  2:20028/1‑100 (MQ=255)                                                    gTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCg                                     <  1:183577/100‑1 (MQ=255)                                                     tGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGa                                    >  1:1201151/1‑100 (MQ=255)                                                       cAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATg                                  >  1:502385/1‑100 (MQ=255)                                                       cAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATg                                  <  2:41438/100‑1 (MQ=255)                                                          gAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATCTCa                               <  2:874418/100‑1 (MQ=255)                                                              ttGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAgg                           <  2:1171801/100‑1 (MQ=255)                                                                gTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCa                         >  2:676812/1‑100 (MQ=255)                                                                gTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCa                         >  1:349319/1‑100 (MQ=255)                                                                   gACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCg                      >  1:1157533/1‑100 (MQ=255)                                                                   gACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCg                      <  2:417611/100‑1 (MQ=255)                                                                      cACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATa                   <  2:69337/100‑1 (MQ=255)                                                                            tGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACa             <  1:724655/100‑1 (MQ=255)                                                                             ggTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACAt            >  2:424581/1‑100 (MQ=255)                                                                                aTCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTca         <  1:370199/100‑1 (MQ=255)                                                                                  cGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTcagc       >  2:1060652/1‑100 (MQ=255)                                                                                   gATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTcagca      >  1:56872/1‑100 (MQ=255)                                                                                       tCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtg  >  2:459563/1‑100 (MQ=255)                                                                                       tCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtg  <  2:952626/100‑1 (MQ=255)                                                                                            |                                                                                               AACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCTG  >  NC_000913/1034235‑1034420

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