Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,142,833	C→T	A53T (GCA→ACA)	hybC ←	hydrogenase 2, large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,142,833	0	C	T	100.0%	106.6 / NA	31	A53T (GCA→ACA)	hybC	hydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/11);  total (20/11)

AATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGAT  >  NC_000913/3142705‑3142967                                                                                                                                 |                                                                                                                                       aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGa                                                                                                                              >  2:389883/1‑139 (MQ=255)           cGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGtttt                                                                                                                    >  1:820476/1‑139 (MQ=255)           cGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGtttt                                                                                                                    >  2:206093/1‑139 (MQ=255)              aTGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCa                                                                                                                 <  2:456914/139‑1 (MQ=255)                gTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACg                                                                                                               <  1:389883/139‑1 (MQ=255)                gTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACg                                                                                                               >  2:1025616/1‑139 (MQ=255)                  aTTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGAt                                                                                                             >  1:483540/1‑139 (MQ=255)                  aTTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGAt                                                                                                             >  2:901710/1‑139 (MQ=255)                        cgTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTc                                                                                                       >  2:873403/1‑139 (MQ=255)                           tAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCAt                                                                                                    >  2:670384/1‑139 (MQ=255)                                      tCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCcaca                                                                                         <  2:483540/139‑1 (MQ=255)                                                 cgcACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGaa                                                                              >  2:488257/1‑139 (MQ=255)                                                                    aCGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAacga                                                           >  1:646365/1‑139 (MQ=255)                                                                               cgcgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCg                                                >  1:348461/1‑139 (MQ=255)                                                                               cgcgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCg                                                >  2:615796/1‑139 (MQ=255)                                                                                gcgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGAt                                                                                                             <  1:1030202/77‑1 (MQ=255)                                                                                gcgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGAt                                                                                                             >  2:1030202/1‑77 (MQ=255)                                                                                 cgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGt                                                                                                                       <  1:844547/66‑1 (MQ=255)                                                                                 cgTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGt                                                                                                                       >  2:844547/1‑66 (MQ=255)                                                                                   tGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGAttt                                            <  2:348461/139‑1 (MQ=255)                                                                                            gtaCATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGa                                   <  2:820476/139‑1 (MQ=255)                                                                                               cATACGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGc                                >  2:508868/1‑139 (MQ=255)                                                                                                  aCGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGa                                               >  1:1100650/1‑121 (MQ=255)                                                                                                  aCGCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGa                                               <  2:1100650/121‑1 (MQ=255)                                                                                                    gCCACAGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTaaa                           >  1:643036/1‑139 (MQ=255)                                                                                                       acaGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGc                        >  2:1008358/1‑139 (MQ=255)                                                                                                         aGATACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGccc                      <  1:670384/139‑1 (MQ=255)                                                                                                           aTACGTTGCACAATCATCCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCt                    <  1:488257/139‑1 (MQ=255)                                                                                                                           tcCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCgg    >  2:1039895/1‑139 (MQ=255)                                                                                                                           tcCATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCgg    <  1:873403/139‑1 (MQ=255)                                                                                                                             cATGTATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGAt  >  2:193807/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       AATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGAT  >  NC_000913/3142705‑3142967

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