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,4120CT91.6% 439.1 / 22.4 131R12H (CGT→CAT) hycHformate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (3/8);  new base T (50/70);  total (53/78)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 131 total aligned reads displayed.

ACCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843364‑2843458
                                                |                                              
 ccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                              >  2:3344834/1‑50 (MQ=255)
  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             >  2:4218646/1‑50 (MQ=255)
  ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGc                                             <  1:4993899/50‑1 (MQ=255)
   tgctgcGCCTCGGCGGGCGTGGCCTACTCCTTATCCATAAATTTATGGCt                                            >  2:3010121/1‑50 (MQ=255)
     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            >  1:119454/1‑48 (MQ=255)
     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                            <  2:931202/48‑1 (MQ=255)
     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                           <  2:1898149/49‑1 (MQ=255)
     ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  2:4635762/50‑1 (MQ=255)
      tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         >  2:1636354/1‑50 (MQ=255)
       gcgcCTCGGCGGGGGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  2:2074453/1‑50 (MQ=255)
       gcgcCTCGGCGGGCGTGGCCTCGTTCTCATCAATAAATTTATGGCTCAGt                                        >  1:1983455/1‑50 (MQ=255)
        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAGTAAATTTATGGCTCAGt                                        <  1:1267482/49‑1 (MQ=255)
        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:4311323/50‑1 (MQ=255)
        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAg                                         <  2:4390062/48‑1 (MQ=255)
        cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGt                                        <  1:769471/49‑1 (MQ=255)
         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:5264904/50‑1 (MQ=255)
         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  2:3918043/50‑1 (MQ=255)
         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:5354476/50‑1 (MQ=255)
         gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:2960384/50‑1 (MQ=255)
          ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                                     <  1:5065657/50‑1 (MQ=255)
          ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                                     <  2:2083081/50‑1 (MQ=255)
           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  1:2065134/50‑1 (MQ=255)
           cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    <  2:3471691/50‑1 (MQ=255)
            tCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  2:477706/1‑50 (MQ=255)
              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:3825638/48‑1 (MQ=255)
              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  <  1:4616555/49‑1 (MQ=255)
              ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:826740/50‑1 (MQ=255)
               tcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  2:4908287/49‑1 (MQ=255)
               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  2:304552/50‑1 (MQ=255)
               gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  1:3178554/50‑1 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGatc                               <  2:4728913/50‑3 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  2:3599339/1‑49 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  2:5358378/1‑49 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:1858261/1‑49 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  1:2420533/1‑49 (MQ=255)
                cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAac                               <  1:1940346/50‑1 (MQ=255)
                 gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              >  2:4862990/1‑50 (MQ=255)
                 gggcgTGGCAGCGTTCTCATCAATAAATTGACGGCTCAGTTGACTGAaca                              <  1:1245916/50‑1 (MQ=255)
                  gtcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  2:4097073/3‑50 (MQ=255)
                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGCacac                             >  1:2823301/1‑50 (MQ=255)
                  ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  2:3704817/50‑1 (MQ=255)
                   gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacc                            <  2:2525806/50‑1 (MQ=255)
                    cgTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAAcacca                           <  2:3664455/50‑1 (MQ=255)
                     gTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAAcaccac                          >  1:168799/1‑50 (MQ=255)
                      tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         >  1:3877950/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:1177750/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:4389909/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:2216983/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:5371394/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  2:835627/50‑1 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:1638916/1‑50 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:3652568/49‑1 (MQ=255)
                        gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  2:2267632/49‑1 (MQ=255)
                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:3370820/1‑50 (MQ=255)
                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  2:4172663/50‑1 (MQ=255)
                         cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:426649/50‑1 (MQ=255)
                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  2:1533175/50‑1 (MQ=255)
                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  2:1895235/50‑1 (MQ=255)
                          aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  2:4349402/1‑50 (MQ=255)
                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  2:4583456/1‑49 (MQ=255)
                           tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                    <  1:4131546/50‑1 (MQ=255)
                           tCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTc                    <  1:4075199/50‑1 (MQ=255)
                            cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  2:2709806/1‑50 (MQ=255)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCGCCTTTTCa                   <  1:2082702/49‑1 (MQ=39)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   <  1:4745893/49‑1 (MQ=255)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  2:3707057/1‑50 (MQ=255)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  2:1075349/1‑50 (MQ=255)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  <  1:1506941/50‑1 (MQ=255)
                             gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTGTCa                   <  2:2907069/49‑1 (MQ=39)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCTCCTTTTCACt                 <  2:667743/50‑1 (MQ=255)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:4354348/1‑50 (MQ=255)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 <  1:3907117/50‑1 (MQ=255)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  2:1150123/1‑50 (MQ=255)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  2:3767211/1‑50 (MQ=255)
                              ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:2648648/1‑50 (MQ=255)
                                 tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  2:2939351/1‑50 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa               <  1:5182051/47‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  2:3893378/48‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:3481822/49‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:5154901/49‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  1:2428061/49‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:4358965/49‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc             <  2:2678472/49‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:3229189/1‑50 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  2:4358284/50‑1 (MQ=255)
                                   atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  1:2948847/1‑50 (MQ=255)
                                   atcaATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            <  2:2605348/50‑1 (MQ=255)
                                    tcaATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg            >  2:2358160/1‑49 (MQ=255)
                                      aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:3946571/1‑49 (MQ=255)
                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:325880/50‑1 (MQ=255)
                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  2:1912216/50‑1 (MQ=255)
                                        tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  2:2605324/50‑1 (MQ=255)
                                         aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  2:2733086/1‑50 (MQ=255)
                                         aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  2:2551723/1‑50 (MQ=255)
                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:5122010/1‑50 (MQ=255)
                                          aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      <  1:4910005/49‑1 (MQ=255)
                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    <  1:5346489/50‑1 (MQ=255)
                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    >  1:1978717/1‑50 (MQ=255)
                                           aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg    >  2:2088910/1‑50 (MQ=255)
                                             ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  2:1248443/1‑50 (MQ=255)
                                                |                                              
ACCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843364‑2843458

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