Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,849,543	G→A	R103C (CGT→TGT)	hycB ←	hydrogenase 3, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,849,543	0	G	A	100.0%	165.3 / NA	46	R103C (CGT→TGT)	hycB	hydrogenase 3, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/20);  total (26/20)

TCATCAAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCAC  >  NC_000913/2849401‑2849689                                                                                                                                               |                                                                                                                                                   tcatcaAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCCTTCGGGGGATTGGCGTTTGCCGGAATATCCAGCGGACAGCTTCCgg                                                                                                                                            >  2:378214/1‑151 (MQ=255)     caAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGaaa                                                                                                                                         <  2:16083/150‑1 (MQ=255)     caAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGAGACGGCGGCGCTTTCCGGGGAGTGCAGGTTGGCGGAATATCCCGCGGACAGCTGCCgcgaa                                                                                                                                         >  2:488820/1‑146 (MQ=255)      aaaGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAAt                                                                                                                                       <  2:598771/151‑1 (MQ=255)           tACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATAACCAGCGGACAGCTGCCGGAAAATTCaa                                                                                                                                   >  2:801111/1‑150 (MQ=255)            aCAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAAtt                                                                                                                                 <  2:83599/151‑1 (MQ=255)            aCAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGGTTGCCGGGATATCCAGCGGGCAGCTGCCGGAAAATTCAAtt                                                                                                                                 >  2:15488/1‑151 (MQ=255)                    cacaTTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCaa                                                                                                                          >  1:371703/1‑150 (MQ=255)                     acaTTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCaaa                                                                                                                         >  1:209873/1‑150 (MQ=255)                         ttGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGCGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATAACCAGCGGACAGCTGCCGGGAAATTCAATTGCGCCAAACggg                                                                                                                     >  2:469627/1‑150 (MQ=255)                           gACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCTGAAAATTCAATTGCGCCAAACGGGCa                                                                                                                   <  2:763700/150‑1 (MQ=255)                            aCGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAgg                                                                                                                 >  1:415765/1‑151 (MQ=255)                             cGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGc                                                                                                                <  2:138214/151‑1 (MQ=255)                                        cgAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGcac                                                                                                      >  1:789544/1‑150 (MQ=255)                                        cgAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGGGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGca                                                                                                       >  2:874917/1‑149 (MQ=255)                                          aaTACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGc                                                                                                   >  1:809547/1‑151 (MQ=255)                                                tGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGAGATGCCGCACAGCTTGc                                                                                               >  1:544490/1‑149 (MQ=255)                                                 ggCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAg                                                                                             <  1:802163/150‑1 (MQ=255)                                                     cccAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTAc                                                                                        >  1:392874/1‑151 (MQ=255)                                                        aGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGca                                                                                     <  1:509343/151‑1 (MQ=255)                                                                     gCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTttcatt                                                                         <  1:592717/150‑1 (MQ=255)                                                                         acacGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCCTACGCACAGGCCTTCATTCaa                                                                      >  2:54020/1‑149 (MQ=255)                                                                               gcCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACATGCTTTTATTCAACTGCa                                                                 >  2:566915/1‑148 (MQ=255)                                                                                      cAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATc                                                       >  2:650056/1‑151 (MQ=255)                                                                                                  gTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGAt                                           >  1:685591/1‑151 (MQ=255)                                                                                                  gTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGAt                                           >  1:320010/1‑151 (MQ=255)                                                                                                          ttCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAGACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTa                                    <  2:276508/150‑1 (MQ=255)                                                                                                           tCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTaa                                   <  1:805469/150‑1 (MQ=255)                                                                                                            cGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAAcc                                 >  1:25230/1‑151 (MQ=255)                                                                                                            cGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAAcc                                 >  1:897903/1‑151 (MQ=255)                                                                                                               ggTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCggg                              <  1:372489/151‑1 (MQ=255)                                                                                                                gTATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGc                             <  1:622927/151‑1 (MQ=255)                                                                                                                 tATTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGc                             >  2:903330/1‑150 (MQ=255)                                                                                                                  aTTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCa                            <  2:807036/150‑1 (MQ=255)                                                                                                                  aTTGGCGTTTGCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCa                            <  2:421128/150‑1 (MQ=255)                                                                                                                            gCCGGAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACcgcg                  <  2:421079/150‑1 (MQ=255)                                                                                                                               ggAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAgg              <  1:632735/151‑1 (MQ=255)                                                                                                                                gAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAggg             <  1:55668/151‑1 (MQ=255)                                                                                                                                gAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAggg             >  1:464333/1‑151 (MQ=255)                                                                                                                                gAATATCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAggg             <  2:789493/151‑1 (MQ=255)                                                                                                                                   tatCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGc          >  1:444792/1‑151 (MQ=255)                                                                                                                                    atCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTgaa         <  2:588778/151‑3 (MQ=255)                                                                                                                                    atCCAGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGc          >  2:854592/1‑150 (MQ=255)                                                                                                                                        aGCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCtt     >  1:537931/1‑151 (MQ=255)                                                                                                                                         gCGGACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTc    <  1:91950/151‑1 (MQ=255)                                                                                                                                            gACAGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTcac  >  2:425427/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                   TCATCAAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCAC  >  NC_000913/2849401‑2849689

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