Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,145,871	T→C	M131V (ATG→GTG)	hybO ←	hydrogenase 2, small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,145,871	0	T	C	100.0%	202.6 / NA	56	M131V (ATG→GTG)	hybO	hydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/30);  total (26/30)

GACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCATGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGAT  >  NC_000913/3145740‑3146005                                                                                                                                    |                                                                                                                                       gACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaataa                                                                                                                                 <  2:375095/139‑1 (MQ=255)         cAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                         <  1:744271/139‑1 (MQ=255)         cAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                         <  2:804968/139‑1 (MQ=255)                gTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCtttt                                                                                                                  >  1:674244/1‑139 (MQ=255)                gTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCtttt                                                                                                                  >  2:784489/1‑139 (MQ=255)                          gCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATgg                                                                                                        >  2:215696/1‑139 (MQ=255)                                gCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaataaata                                                                                                                              >  2:1568983/1‑111 (MQ=255)                                gCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaataaata                                                                                                                              <  1:1568983/111‑1 (MQ=255)                                    cacCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccac                                                                                              <  2:832609/139‑1 (MQ=255)                                      ccGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                            <  2:468406/139‑1 (MQ=255)                                                cAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATa                                                                                  <  1:784489/139‑1 (MQ=255)                                                 agagCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaat                                                                                                                                   <  1:1405641/89‑1 (MQ=255)                                                 agagCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaat                                                                                                                                   >  2:1405641/1‑89 (MQ=255)                                                     cAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGcc                                                                             <  2:978082/139‑1 (MQ=255)                                                         aaCCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCttt                                                                         >  1:969766/1‑139 (MQ=255)                                                          aCCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTg                                                                        <  2:1629276/139‑1 (MQ=255)                                                                   gATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGt                                                                                                                          <  2:882043/80‑1 (MQ=255)                                                                   gATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGt                                                                                                                          >  1:882043/1‑80 (MQ=255)                                                                    aTAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                <  2:494353/105‑1 (MQ=255)                                                                    aTAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                >  1:494353/1‑105 (MQ=255)                                                                    aTAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCga                                                              <  2:674244/139‑1 (MQ=255)                                                                       aTGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGagag                                                           <  2:1085343/139‑1 (MQ=255)                                                                            tgctgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGt                                                        <  1:976581/137‑1 (MQ=255)                                                                            tgctgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGt                                                        >  2:976581/1‑137 (MQ=255)                                                                               tgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                            >  2:1498023/1‑98 (MQ=255)                                                                               tgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                            <  1:1498023/98‑1 (MQ=255)                                                                               tgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                   <  2:969766/139‑1 (MQ=255)                                                                               tgcGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgt                                                   >  2:1413862/1‑139 (MQ=255)                                                                                   ccTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTa                                                                                                                        >  1:941239/1‑66 (MQ=255)                                                                                   ccTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTa                                                                                                                        <  2:941239/66‑1 (MQ=255)                                                                                    cTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTAtgtttgtt                                              <  1:1317869/139‑1 (MQ=255)                                                                                      tCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATc                                                                                                 <  2:230027/86‑1 (MQ=255)                                                                                      tCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATc                                                                                                 >  1:230027/1‑86 (MQ=255)                                                                                        cgccgcTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGc                                                                                                                                      >  2:118582/1‑47 (MQ=255)                                                                                        cgccgcTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGc                                                                                                                                      <  1:118582/47‑1 (MQ=255)                                                                                                tGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtg                                  >  2:189005/1‑139 (MQ=255)                                                                                                  cGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccac                                                                                              >  1:775910/1‑77 (MQ=255)                                                                                                  cGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccac                                                                                              <  2:775910/77‑1 (MQ=255)                                                                                                  cGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCtgat                                >  1:1148714/1‑139 (MQ=255)                                                                                                        tgATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATg                                                                                                         <  1:973936/60‑1 (MQ=255)                                                                                                        tgATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATg                                                                                                         >  2:973936/1‑60 (MQ=255)                                                                                                        tgATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCgag                                                             >  2:115744/1‑104 (MQ=255)                                                                                                        tgATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCgag                                                             <  1:115744/104‑1 (MQ=255)                                                                                                          aTCCACAATCGGCTCACCAGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAg                        >  2:1099226/1‑139 (MQ=255)                                                                                                           tCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCGCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAgg                       <  1:1099226/139‑1 (MQ=255)                                                                                                           tCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAgg                       >  2:1315323/1‑139 (MQ=255)                                                                                                                    ggCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                >  2:21101/1‑57 (MQ=255)                                                                                                                    ggCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                <  1:21101/57‑1 (MQ=255)                                                                                                                     gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                            >  2:1648867/1‑60 (MQ=255)                                                                                                                     gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                            <  1:1648867/60‑1 (MQ=255)                                                                                                                     gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAg             >  2:675606/1‑139 (MQ=255)                                                                                                                       tCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCa           <  1:215696/139‑1 (MQ=255)                                                                                                                            ggCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCg      <  1:1413862/139‑1 (MQ=255)                                                                                                                              cAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACt                                                                                >  1:470916/1‑63 (MQ=255)                                                                                                                              cAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACt                                                                                <  2:470916/63‑1 (MQ=255)                                                                                                                                aCCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGAt  >  1:1533928/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       GACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCATGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGAT  >  NC_000913/3145740‑3146005

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