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_0009133,145,8710TC100.0% 152.1 / NA 42M131V (ATG→GTG) hybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/22);  total (20/22)

ACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCATGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATA  >  NC_000913/3145741‑3146006
                                                                                                                                  |                                                                                                                                       
aCTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaataaa                                                                                                                                 <  2:877626/139‑1 (MQ=255)
aCTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCaataaa                                                                                                                                 >  1:877626/1‑139 (MQ=255)
       cAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                          <  2:513660/139‑1 (MQ=255)
       cAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                          <  1:506516/139‑1 (MQ=255)
          ttGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATc                                                                                                                       >  2:612964/1‑139 (MQ=255)
                       aGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGc                                                                                                                                       <  2:887398/110‑1 (MQ=255)
                       aGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGc                                                                                                                                       >  1:887398/1‑110 (MQ=255)
                              gCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCAt                                                                                                   >  2:751776/1‑139 (MQ=255)
                               cAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATc                                                                                                  <  1:612964/139‑1 (MQ=255)
                                    ccGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                             <  1:516825/139‑1 (MQ=255)
                                    ccGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccacca                                                                                             <  1:818653/139‑1 (MQ=255)
                                         ccACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAac                                                                                        <  2:886705/139‑1 (MQ=255)
                                            cGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                          >  1:690959/1‑102 (MQ=255)
                                            cGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGtt                                                                                                                          <  2:690959/102‑1 (MQ=255)
                                                      gAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCtt                                                                           >  1:847842/1‑139 (MQ=255)
                                                        aCCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTg                                                                         >  1:500422/1‑139 (MQ=255)
                                                          cGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTa                                                                       <  1:751776/139‑1 (MQ=255)
                                                           gATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTAc                                                                      >  2:503875/1‑139 (MQ=255)
                                                                  aTAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATg                                                                                                          <  1:601231/96‑1 (MQ=255)
                                                                  aTAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATg                                                                                                          >  2:601231/1‑96 (MQ=255)
                                                                                    tCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTg                                                                         <  2:676002/111‑1 (MQ=255)
                                                                                    tCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTg                                                                         >  1:676002/1‑111 (MQ=255)
                                                                                     ccgccgCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                 <  1:803019/86‑1 (MQ=255)
                                                                                     ccgccgCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                 >  2:803019/1‑86 (MQ=255)
                                                                                      cgccgcTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCg                                                                                                                            <  1:240429/58‑1 (MQ=255)
                                                                                      cgccgcTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCg                                                                                                                            >  2:240429/1‑58 (MQ=255)
                                                                                        ccgcTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTc                                         <  1:503875/139‑1 (MQ=255)
                                                                                                cGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                 >  1:853822/1‑75 (MQ=255)
                                                                                                cGGATGTGATCCACAATCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATcc                                                                                                 <  2:853822/75‑1 (MQ=255)
                                                                                                               aTCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccac                                                                                               >  2:594828/1‑62 (MQ=255)
                                                                                                               aTCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATccac                                                                                               <  1:594828/62‑1 (MQ=255)
                                                                                                                tCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTaa                                                                                                                 >  1:235989/1‑43 (MQ=255)
                                                                                                                tCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTaa                                                                                                                 <  2:235989/43‑1 (MQ=255)
                                                                                                                tCGGCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCAAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGa                 >  1:727331/1‑139 (MQ=255)
                                                                                                                   gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAg              <  2:177565/139‑1 (MQ=255)
                                                                                                                   gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAg              >  1:498957/1‑139 (MQ=255)
                                                                                                                   gCTCACCGGCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAg              >  2:826564/1‑139 (MQ=255)
                                                                                                                          ggCAACCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCg       <  2:498957/139‑1 (MQ=255)
                                                                                                                              aCCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACtt         >  1:437962/1‑133 (MQ=255)
                                                                                                                              aCCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACtt         <  2:437962/133‑1 (MQ=255)
                                                                                                                              aCCACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGAt   >  2:889777/1‑139 (MQ=255)
                                                                                                                               ccACGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATa  <  1:826564/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
ACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACCGCCCCACGCAGAGCAGGAACCGATAGCGATAATGGCTGCTGCGCCTTCCGCCGCTTTGCGGATGTGATCCACAATCGGCTCACCGGCAACCATGCAATAAATACCGTTATCTTTTAATGGGATGGAACCATCCACCACTAACACATACTGCCCTTTGTACTTCTCGAGAGCGTTATGTTTGTTCTCTTCGACCTGATGACCGAAGGCGGCGGAAAGCACTTCGTGATA  >  NC_000913/3145741‑3146006

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