Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,141,584	A→G	F469S (TTC→TCC)	hybC ←	hydrogenase 2 large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,141,584	0	A	G	100.0%	116.3 / NA	34	F469S (TTC→TCC)	hybC	hydrogenase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

GCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGTGCACTG  >  NC_000913/3141508‑3141661                                                                             |                                                                              gCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGAtt                                                                            >  2:521334/1‑80 (MQ=255)     atgatACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCAccc                                                                       <  1:521301/81‑1 (MQ=255)     atgatACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCAcc                                                                        <  2:521299/80‑1 (MQ=255)           ccGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGcc                                                                 >  2:521345/1‑81 (MQ=255)            cGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCg                                                                >  1:521344/1‑81 (MQ=255)                       aCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCgg                                                      <  1:521303/80‑1 (MQ=255)                       aCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCgg                                                      <  2:521301/80‑1 (MQ=255)                         cATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCtt                                                   >  1:521353/1‑81 (MQ=255)                          aTCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTc                                                  <  1:521304/81‑1 (MQ=255)                                  agagagCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGTCGGAATGTTCGGCTTCACAAAggt                                          >  2:521338/1‑81 (MQ=255)                                  agagagCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAggt                                          >  2:521344/1‑81 (MQ=255)                                  agagagCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAggt                                          <  1:521305/81‑1 (MQ=255)                                     gagCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAggtggt                                       >  2:521353/1‑81 (MQ=255)                                        cATAACGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGtgtg                                    <  2:521305/81‑1 (MQ=255)                                            ccGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATc                                 <  2:521306/80‑1 (MQ=255)                                            ccGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATc                                 <  2:521308/80‑1 (MQ=255)                                            ccGCGCGGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCg                                <  1:521308/81‑1 (MQ=255)                                               cgcgGCGCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCt                             <  1:521310/81‑1 (MQ=255)                                                gcggcgCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCtt                            >  1:521364/1‑81 (MQ=255)                                                 cggcgCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCttt                           <  2:521310/81‑1 (MQ=255)                                                  ggcgCTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTg                          <  1:521312/81‑1 (MQ=255)                                                    cgcTTCGAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGcc                        <  1:521313/81‑1 (MQ=255)                                                          gAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATAtt                  <  1:521314/81‑1 (MQ=255)                                                          gAGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATAtt                  >  1:521349/1‑81 (MQ=255)                                                           aGGAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTg                 >  2:521349/1‑81 (MQ=255)                                                             gAAGCCAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGt               >  1:521343/1‑81 (MQ=255)                                                                 ccAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGAt            >  1:521358/1‑80 (MQ=255)                                                                 ccAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATc           <  2:521315/81‑1 (MQ=255)                                                                 ccAACACCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATc           <  1:521315/81‑1 (MQ=255)                                                                    acacCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGt        <  1:521317/81‑1 (MQ=255)                                                                    acacCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGt        >  2:521354/1‑81 (MQ=255)                                                                      acCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGTg       >  2:521355/1‑80 (MQ=255)                                                                      acCTTTGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGTGc      <  2:521317/81‑1 (MQ=255)                                                                          ttGGATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGTGCACTg  >  2:521347/1‑81 (MQ=255)                                                                             |                                                                              GCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTCGGCTTCACAAAGGTGGTGTGATCGCCTTTGCCGATATTGGTGATCAGTGCACTG  >  NC_000913/3141508‑3141661

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