Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,843,412 C→T R12H (CGT→CAT)  hycH ← formate hydrogenlyase assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,4120CT100.0% 264.2 / NA 65R12H (CGT→CAT) hycHformate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/33);  total (32/33)

CGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCAC  >  NC_000913/2843380‑2843442
                                |                              
cgggcgTGGCATCGTTCTCATCAATAAATTTATgg                              <  1:6297893/35‑1 (MQ=255)
  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCt                            <  1:4405815/35‑1 (MQ=255)
    cgTGGCATCGTTCTCATCAATAAATTTATGGCTCa                          >  1:4510780/1‑35 (MQ=255)
     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                         <  1:5281817/35‑1 (MQ=38)
     gTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                         <  1:7493141/35‑1 (MQ=38)
      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                        >  1:7123316/1‑35 (MQ=38)
      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                        <  1:2391919/35‑1 (MQ=255)
      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                        >  1:6149197/1‑35 (MQ=255)
       ggCATCGTTCTCATCAATAAATTTATGGCTCAGtt                       <  1:1597301/35‑1 (MQ=255)
       ggCATCGTTCTCATCAATAAATTTATGGCTCAGtt                       <  1:266965/35‑1 (MQ=38)
       ggCATCGTTCTCATCAATAAATTTATGGCTCAGtt                       >  1:5170944/1‑35 (MQ=38)
       ggCATCGTTCTCATCAATAAATTTATGGCTCAGtt                       <  1:608049/35‑1 (MQ=255)
        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                      >  1:4440528/1‑35 (MQ=38)
        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                      <  1:2680025/35‑1 (MQ=255)
        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                      >  1:917242/1‑35 (MQ=255)
        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                      >  1:930617/1‑35 (MQ=255)
         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                     >  1:5148602/1‑35 (MQ=255)
          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                    <  1:2800646/35‑1 (MQ=255)
          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                    <  1:2631984/35‑1 (MQ=38)
          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                    <  1:4543444/35‑1 (MQ=38)
           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                   <  1:6701840/35‑1 (MQ=255)
           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                   <  1:7870632/35‑1 (MQ=255)
            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                  <  1:1930951/35‑1 (MQ=255)
             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                 >  1:2996614/1‑35 (MQ=38)
             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                 >  1:2792825/1‑35 (MQ=255)
              ttCTCATCAATAAATTTATGGCTCAGTTGACTGaa                >  1:515313/1‑35 (MQ=38)
              ttCTCATCAATAAATTTATGGCTCAGTTGACTGaa                >  1:6094863/1‑35 (MQ=38)
              ttCTCATCAATAAATTTATGGCTCAGTTGACTGaa                >  1:7718303/1‑35 (MQ=255)
                ctcATCAATAAATTTATGGCTCAGTTGACTGAaca              >  1:4490723/1‑35 (MQ=255)
                ctcATCAATAAATTTATGGCTCAGTTGACTGAaca              >  1:4995896/1‑35 (MQ=255)
                ctcATCAATAAATTTATGGCTCAGTTGACTGAaca              >  1:847034/1‑35 (MQ=38)
                ctcATCAATAAATTTATGGCTCAGTTGACTGAaca              >  1:2041290/1‑35 (MQ=38)
                 tcatcaATAAATTTATGGCTCAGTTGACTGAacac             <  1:5871668/35‑1 (MQ=255)
                 tcatcaATAAATTTATGGCTCAGTTGACTGAacac             >  1:202145/1‑35 (MQ=255)
                 tcatcaATAAATTTATGGCTCAGTTGACTGAacac             >  1:4230995/1‑35 (MQ=255)
                  catcaATAAATTTATGGCTCAGTTGACTGAAcacc            >  1:140409/1‑35 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           >  1:5424021/1‑35 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           <  1:1926174/35‑1 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           <  1:4498098/35‑1 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           <  1:205935/35‑1 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           <  1:6245103/35‑1 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           >  1:2303664/1‑35 (MQ=255)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           >  1:7107712/1‑35 (MQ=39)
                   atcaATAAATTTATGGCTCAGTTGACTGAAcacca           >  1:7022949/1‑35 (MQ=39)
                    tcaATAAATTTATGGCTCAGTTGACTGAAcaccac          <  1:5626600/35‑1 (MQ=39)
                    tcaATAAATTTATGGCTCAGTTGACTGAAcaccac          >  1:6788436/1‑35 (MQ=39)
                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCt        <  1:2023237/35‑1 (MQ=255)
                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCt        <  1:7685939/35‑1 (MQ=255)
                       aTAAATTTATGGCTCAGTTGGCTGAACACCACCtt       <  1:5880886/35‑1 (MQ=25)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       <  1:2655554/35‑1 (MQ=255)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       <  1:1278243/35‑1 (MQ=255)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       <  1:1650976/35‑1 (MQ=255)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       >  1:4417026/1‑35 (MQ=255)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       >  1:5672242/1‑35 (MQ=38)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       >  1:483093/1‑35 (MQ=38)
                       aTAAATTTATGGCTCAGTTGACTGAACACCACCtt       <  1:4862323/35‑1 (MQ=255)
                        tAAATTTATGGCTCAGTTGACTGAACACCACCttt      <  1:6950261/35‑1 (MQ=255)
                        tAAATTTATGGCTCAGTTGACTGAACACCACCttt      >  1:4041301/1‑35 (MQ=38)
                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTc    >  1:6634014/1‑35 (MQ=38)
                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTc    <  1:4728041/35‑1 (MQ=255)
                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCa   >  1:2004260/1‑35 (MQ=38)
                            tttATGGCTCAGTTGACTGAACACCACCTTTTCAc  >  1:4920467/1‑35 (MQ=38)
                            tttATGGCTCAGTTGACTGAACACCACCTTTTCAc  >  1:6643719/1‑35 (MQ=38)
                            tttATGGCTCAGTTGACTGAACACCACCTTTTCAc  <  1:3838469/35‑1 (MQ=255)
                            tttATGGCTCAGTTGACTGAACACCACCTTTTCAc  <  1:7653063/35‑1 (MQ=255)
                            tttATGGCTCAGTTGACTGAACACCACCTTTTCAc  <  1:806433/35‑1 (MQ=38)
                                |                              
CGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCAC  >  NC_000913/2843380‑2843442

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