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% 172.0 / NA 43R12H (CGT→CAT) hycHformate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/21);  total (22/21)

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843365‑2843457
                                               |                                             
ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                             <  1:2216618/50‑1 (MQ=255)
 ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                            <  1:856849/50‑1 (MQ=255)
   gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           >  1:3357359/1‑49 (MQ=255)
    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                           <  1:3731122/48‑1 (MQ=255)
    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                         <  1:1099529/50‑1 (MQ=255)
        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                      >  1:2347541/1‑49 (MQ=255)
        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                     <  1:283443/50‑1 (MQ=255)
          cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                   <  1:3268626/50‑1 (MQ=255)
           tCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                  >  1:3816481/1‑50 (MQ=255)
             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                  <  1:2956240/48‑1 (MQ=255)
             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                 >  1:1567608/1‑49 (MQ=255)
              gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               <  1:3821215/50‑1 (MQ=255)
              gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                               >  1:3748898/1‑50 (MQ=255)
                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                             <  1:109801/50‑1 (MQ=255)
                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                             >  1:2912038/1‑50 (MQ=255)
                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            >  1:595088/1‑50 (MQ=255)
                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            >  1:2352681/1‑49 (MQ=255)
                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                            >  1:1295499/1‑49 (MQ=255)
                   cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                          >  1:1942018/1‑50 (MQ=255)
                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                      <  1:1591400/50‑1 (MQ=255)
                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                       <  1:614341/49‑1 (MQ=255)
                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                    >  1:3152194/1‑50 (MQ=255)
                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                   <  1:1893484/50‑1 (MQ=255)
                           ctttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  >  1:1705517/3‑50 (MQ=255)
                           cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  >  1:1251896/1‑50 (MQ=255)
                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                  <  1:1797022/49‑1 (MQ=255)
                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 <  1:108080/50‑1 (MQ=255)
                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:2684115/1‑50 (MQ=255)
                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                 >  1:3512076/1‑50 (MQ=255)
                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:2354869/1‑50 (MQ=255)
                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:2888300/1‑50 (MQ=255)
                              tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:3416980/1‑49 (MQ=255)
                              tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                >  1:2027365/1‑49 (MQ=255)
                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa              <  1:2392022/47‑1 (MQ=255)
                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc            <  1:723806/49‑1 (MQ=255)
                                   tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg          <  1:482559/50‑1 (MQ=255)
                                    caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         >  1:3577953/1‑50 (MQ=255)
                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa         <  1:2693593/49‑1 (MQ=255)
                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt        <  1:509931/50‑1 (MQ=255)
                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc      >  1:3518947/1‑50 (MQ=255)
                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc      <  1:2063152/50‑1 (MQ=255)
                                          aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg   <  1:770593/50‑1 (MQ=255)
                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc  >  1:3682154/1‑50 (MQ=255)
                                               |                                             
CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGC  >  NC_000913/2843365‑2843457

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