Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,448,543 +C intergenic (+48/‑90) ytfQ → / → ytfR putative LACI‑type transcriptional regulator/YtfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,5421.C97.9% 229.8 / ‑0.5 96intergenic (+47/‑91)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (62/32);  total (63/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGC  >  NC_000913/4448478‑4448610
                                                                  |                                                                     
aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                             >  1:2187051/1‑76 (MQ=25)
 aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                            <  2:1255121/76‑1 (MQ=25)
  aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                           >  2:987746/1‑76 (MQ=25)
  aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                           >  1:1706126/1‑76 (MQ=25)
  aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                           >  2:2503734/1‑76 (MQ=25)
  aaaGATGGGGTATTGGGGTGTGCTATGCCTGATGCCGATTCGTAGGCCGGTTAAGGCGCTCGCGCCGCATCCGGCg                                                           <  2:1612175/76‑1 (MQ=17)
    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcatt                                                         >  2:3415414/1‑73 (MQ=21)
    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAt                                                         >  1:1373035/1‑76 (MQ=25)
    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAt                                                         >  1:3215376/1‑76 (MQ=25)
     tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                        <  2:2890995/76‑1 (MQ=25)
      atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                       >  2:3567460/1‑76 (MQ=25)
      atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                       >  2:3224235/1‑76 (MQ=25)
      atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGACgg                                                       >  1:3516982/1‑76 (MQ=21)
       tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                      <  1:1179815/76‑1 (MQ=25)
        ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                     <  2:1465640/76‑1 (MQ=25)
         gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCCTCCGGCGATGGGGc                                                    >  2:238068/1‑76 (MQ=21)
         gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                    >  1:3740216/1‑76 (MQ=25)
         gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                    <  1:3695471/76‑1 (MQ=25)
          ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                   >  2:3921908/1‑76 (MQ=25)
           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                  >  1:3845456/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGCTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGTGGCACt                                                 >  2:1456540/1‑76 (MQ=18)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                 >  1:2109300/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                 >  2:1437884/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                 >  2:2098347/1‑76 (MQ=25)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCCTCCGGCGATGGTGCACTg                                                >  2:1716693/1‑76 (MQ=21)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                >  1:3413047/1‑76 (MQ=25)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                >  1:701596/1‑76 (MQ=25)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                >  1:3051415/1‑76 (MQ=25)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                >  1:3025234/1‑76 (MQ=25)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                              <  2:4124228/76‑1 (MQ=21)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                              >  2:4122953/1‑76 (MQ=25)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                              >  2:3018356/1‑76 (MQ=25)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                              >  2:2721019/1‑76 (MQ=25)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCTGCATCCGGCGATGGTGCACTGaa                                              >  2:2385072/1‑76 (MQ=25)
                gAGGGTTGCTATGCCTGATGCTGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                             >  2:1955528/1‑76 (MQ=21)
                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                             >  1:1472479/1‑76 (MQ=25)
                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                             <  2:2406757/76‑1 (MQ=21)
                   ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                          >  1:145960/1‑76 (MQ=25)
                   ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                          <  2:3762852/76‑1 (MQ=25)
                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                         >  1:2226689/1‑76 (MQ=25)
                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                         >  1:52655/1‑76 (MQ=25)
                     ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                        >  2:1575129/1‑76 (MQ=21)
                     ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                        >  2:3338077/1‑76 (MQ=25)
                     ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                        <  2:3516662/76‑1 (MQ=21)
                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                       >  1:222234/1‑76 (MQ=25)
                       gggaTGCCTGATGCCGGTTCGTGGGCCGGTTAAGGCGCTCGCGCCGCGTCCGGCGATGGTGCACTGAAGCCTGATg                                      <  2:4349330/73‑1 (MQ=18)
                       gCTGTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                      <  2:4567735/76‑1 (MQ=21)
                        cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGCGGGGGCACTGCAGCCCGCTGc                                     >  2:919307/1‑76 (MQ=9)
                        cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                     >  1:4668642/1‑76 (MQ=21)
                        cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                     <  1:1906713/76‑1 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                    >  1:4854678/1‑76 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                    >  1:371384/1‑76 (MQ=21)
                          aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGTTGCGa                                   <  1:1392294/76‑1 (MQ=21)
                              cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                               <  1:4004728/76‑1 (MQ=21)
                              cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                               >  2:1604754/1‑76 (MQ=21)
                               tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                             >  2:4416917/1‑61 (MQ=18)
                               tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                             <  1:4416917/61‑1 (MQ=18)
                                gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                             <  1:416030/76‑1 (MQ=21)
                                gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                             <  2:2688923/76‑1 (MQ=21)
                                gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                             <  1:4464802/76‑1 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                           >  2:4218502/1‑76 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                           <  2:408277/76‑1 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                           >  1:3279991/1‑76 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                           >  1:1102401/1‑76 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                           <  2:4224790/76‑1 (MQ=21)
                                   gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                          <  1:3930208/76‑1 (MQ=21)
                                   gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                          >  1:2577619/1‑76 (MQ=21)
                                    ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGTGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                         >  1:3243150/1‑76 (MQ=21)
                                     cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                        >  1:3981010/1‑76 (MQ=21)
                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcgc                       >  2:2526542/1‑75 (MQ=21)
                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                       >  1:4090780/1‑76 (MQ=21)
                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                       <  2:810724/76‑1 (MQ=21)
                                       aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                      <  1:2769159/76‑1 (MQ=21)
                                        ttCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                     >  2:1364117/1‑76 (MQ=21)
                                        ttCGTAGGCCGGATAAGGAGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTCCCGCGTCt                     >  2:1838914/1‑76 (MQ=18)
                                           gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                  >  2:3770963/1‑76 (MQ=21)
                                            tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                 >  1:2642467/1‑76 (MQ=21)
                                            tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                 >  1:2942965/1‑76 (MQ=21)
                                              ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt               >  1:2980077/1‑76 (MQ=18)
                                              ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt               >  1:4049408/1‑76 (MQ=18)
                                               gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg              >  1:3410944/1‑76 (MQ=18)
                                                ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc             <  1:99975/76‑1 (MQ=21)
                                                ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc             <  1:1370359/76‑1 (MQ=18)
                                                ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc             <  2:3290030/76‑1 (MQ=21)
                                                 cGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAAGCTTACCGCGTCTGATCCTGcc            >  2:3954980/1‑76 (MQ=18)
                                                  ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt           >  2:2688720/1‑76 (MQ=21)
                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGACGCCTGATGCGACGCTTACCGCGTCTTATCATGTCtc          >  2:974616/1‑75 (MQ=18)
                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCtc          >  2:938096/1‑75 (MQ=21)
                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa          <  2:4161577/76‑1 (MQ=18)
                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa          <  2:1627324/76‑1 (MQ=21)
                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa          <  2:2566370/76‑1 (MQ=21)
                                                   gATAAGGCGCGCGCGCGGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa          <  2:775602/76‑1 (MQ=21)
                                                     tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt        <  1:1994800/76‑1 (MQ=21)
                                                      aaGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg       <  2:1417574/76‑1 (MQ=21)
                                                      aaGGCGCTCGAGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg       >  2:871565/1‑76 (MQ=21)
                                                       aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg      >  2:984341/1‑76 (MQ=21)
                                                         gcgcgcgcCGCATCCGGCGATGGGGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc  <  2:2142809/76‑1 (MQ=25)
                                                                  |                                                                     
AAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGC  >  NC_000913/4448478‑4448610

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