Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 683,496 A→C V65G (GTT→GGT)  nagC DNA‑binding transcriptional dual regulator, repressor of N‑acetylglucosamine

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606683,4960A→C100.0% 194.8 60V65G (GTT→GGT) nagCDNA‑binding transcriptional dual regulator, repressor of N‑acetylglucosamine
Reads supporting (aligned to +/- strand):  new base (24/36):  ref base (0/0):  total (24/36)

GCGCGGCGGCCCCCAGTGGAGGCCTGCTGATCAACCTCTTTGATCAGCCCGCGTTCGATAAGCTGA  >  REL606/683463‑683528
                                 |                                
gcgcGGCGGCCCCCAGTGGAGGCCTGCTGATCACcc                                <  1:14596/36‑1 (MQ=255)
 cgcgGCGGCCCCCAGTGGAGGCCTGCTGATCACCct                               <  1:3601440/36‑1 (MQ=255)
 cgcgGCGGCCCCCAGTGGAGGCCTGCTGATCACCct                               <  1:6727403/36‑1 (MQ=255)
 cgcgGCGGCCCCCAGTGGAGGCCTGCTGATCACCct                               <  1:3166465/36‑1 (MQ=255)
 cgcgGCGGCCCCCAGTGGAGGCCTGCTGATCACCct                               <  1:5298743/36‑1 (MQ=255)
  gcggcgGCCCCCAGTGGAGGCCTGCTGATCACCctc                              <  1:5477864/36‑1 (MQ=255)
  gcggccGCCCCCAGTGGAGGCCTGCTGATCACCctc                              <  1:4646999/36‑1 (MQ=25)
   cggcggCCCCCAGTGGAGGCCTGCTGATCACCctct                             <  1:5742518/36‑1 (MQ=255)
   cggcggCCCCCAGTGGAGGCCTGCTGATCACCctct                             >  1:5676174/1‑36 (MQ=255)
    ggcggcCCCCAGTGGAGGCCTGCTGATCACCCTCtt                            <  1:2292446/36‑1 (MQ=255)
      cggcCCCCAGTGGAGGCCTGCTGATCACCCTCTTTg                          >  1:6644795/1‑36 (MQ=255)
      ccgcCCCCAGTGGAGGCCTGCTGATCACCCTCTTTg                          <  1:6655776/34‑1 (MQ=25)
       ggcCCCCAGTGGAGGCCTGCTGATCACCCTCTTTGa                         <  1:6580951/36‑1 (MQ=255)
       ggcCCCCAGTGGAGGCCTGCTGATCACCCTCTTTGa                         <  1:5640709/36‑1 (MQ=255)
       ggcCCCCAGTGGAGGCCTGCTGAGCACCCTCTTTGa                         <  1:5876441/36‑1 (MQ=25)
         cccccAGTGGAGGCCTGCTGATCACCCTCTTTGATc                       <  1:6859015/36‑1 (MQ=255)
         cccccAGTGGAGGCCTGCTGATCACCCTCTTTGATc                       <  1:3163302/36‑1 (MQ=255)
          ccccGGTGGAGGCCTGCTGATCACCCTCTTTGATCa                      <  1:4588575/36‑1 (MQ=25)
          ccccAGTGGAGGCCTGCTGATCACCCTCTTTGATCa                      <  1:6147572/36‑1 (MQ=255)
          ccccAGTGGAGGCCTGCTGATCACCCTCTTTGATCa                      <  1:1690733/36‑1 (MQ=255)
           cccAGTGGAGGCCTGCTGATCACCCTCTTTGATCAg                     >  1:322723/1‑36 (MQ=255)
           cccAGGGGAGGCCTGCTGATCACCCTCTGGGATCAg                     <  1:4795073/36‑1 (MQ=18)
             cggTGGAGGCCTGCTGATCACCCTCTTTGATCAGcc                   <  1:2410003/34‑1 (MQ=25)
             cAGTGGAGGCCTGCTGATCACCCTCTTTGATCAGcc                   >  1:2333921/1‑36 (MQ=255)
             cAGTGGAGGCCTGCTGATCACCCTCTTTGATCAGcc                   >  1:2160074/1‑36 (MQ=255)
             cAGGGGAGGCCTGCTGATCACCCTCTTTGATCAGcc                   <  1:4870768/36‑1 (MQ=25)
              aGTGGAGGCCTGCTGATCACCCTCTTTGATCAGccc                  >  1:1814721/1‑36 (MQ=255)
              aGTGGAGGCCTGCTGATCACCCTCTTTGATCAGccc                  <  1:6248248/36‑1 (MQ=255)
              aGTGGAGGCCTGCTCATCACCCTCTTTTATCAAccc                  >  1:2412454/1‑36 (MQ=18)
               gTGGAGGCCTGCTGATCACCCTCTTTGATCAGCCcg                 <  1:1754453/36‑1 (MQ=255)
               gTGGAGGACTGCTGATCACCCTCTTTGATCAGCCcg                 >  1:208724/1‑36 (MQ=25)
                tGGAGGCCTGCTGATCACCCTCTTTGATCAGCCcgc                >  1:3505569/1‑36 (MQ=255)
                 ggAGGCCTGCTGATCACCCTCTTTGATCAGCCcgcg               >  1:4745169/1‑36 (MQ=255)
                  gAGGCCTGCTGATCACCCTCTTTGAGCAGCCCGCGt              <  1:5597434/36‑1 (MQ=25)
                   agcccTGCTGATCACCCTCTTTGATCAGCCCGCGtt             <  1:7206759/33‑1 (MQ=25)
                   aGGGCTGCTGATCACCCTCTTTGATCAGCCCGCGtt             <  1:5720218/36‑1 (MQ=25)
                     gCCTGCTGATCACCCTCTTTGATCAGCCCGCGTTCg           <  1:6483497/36‑1 (MQ=255)
                     gCCTGCTGATCACCCTCTTTGATCAGCCCGCGTTCg           <  1:3319300/36‑1 (MQ=255)
                       ctgctgGTCACCGTCTGTGATCAGCCCGCGTTCGAt         <  1:95532/36‑1 (MQ=18)
                       ctgctgATCACCCTCTTTGATCAGCCCGCGTTCGAt         >  1:5736530/1‑36 (MQ=255)
                       ctgctgATCACCCTCTTTGATCAGCCCGCGTTCGAt         <  1:2049290/36‑1 (MQ=255)
                       ctgctgATCACCCTCTTTGATCAGCCCGCGTTCGAt         <  1:4983649/36‑1 (MQ=255)
                        tgctgGTCACCCTCTTTGATCAGCCCGCGTTCGATa        <  1:5636239/36‑1 (MQ=25)
                        tgctgATCACCCTCTTTGATCAGCCCGCGTTCGATa        >  1:6155337/1‑36 (MQ=255)
                        tgctgATCACCCTCTTTGATCAGCCCGCGTTCGATa        >  1:989125/1‑36 (MQ=255)
                        tgctgATCACCCTCTTTGATCAGCCCGCGTTCGATa        >  1:5797172/1‑36 (MQ=255)
                        tgctgATCACCCTCTTTGATCAGCCCGCGTTCGATa        >  1:2964142/1‑36 (MQ=255)
                          ctgATCACCCTCTTTAATCAGCCCGCGTTCGATAAg      >  1:1616952/1‑36 (MQ=25)
                           tgATCACCCTCTTTGATCGGCCCGCGTTCGATAAGc     <  1:1096693/36‑1 (MQ=25)
                           tgATCACCCTCTTTGATCAGCCCGCGTTCGATAAGc     >  1:5809441/1‑36 (MQ=255)
                            gATCACGCTCTTTGATCAGCCCGCGTTCGATAAGCt    <  1:622419/36‑1 (MQ=25)
                            gATCACCCTCTTTGATCAGCCCGCGTTTGATAAGCt    >  1:3141840/1‑36 (MQ=25)
                            gATCACCCTCTTTGATCAGCCCGCGTTCGATATGCt    >  1:4083917/1‑36 (MQ=25)
                            gATCACCCTCTTTGATCAGCCCGCGTTCGATAAGCt    >  1:5381037/1‑36 (MQ=255)
                            gATCACCCTCTTTGATCAGCCCGCGTTCGATAAGCt    >  1:476776/1‑36 (MQ=255)
                            gATCACCCTCTTTGATCAGCCCGCGTTCGAAAAGCt    >  1:6515816/1‑36 (MQ=25)
                             aTCACCCTCTTTGATCAGCCCGCGTTCGATAAGCTg   >  1:5890251/1‑36 (MQ=255)
                             aTCACCCTCTTTGATCAGCCCGCGTTCGATAAGCTg   <  1:5789558/36‑1 (MQ=255)
                             aTCACCCTCTTTGATCAGCCCGCGTTCGATAAGCTg   <  1:4059865/36‑1 (MQ=255)
                              tCACCCTCTTTGATCAGCCCGCGTTCGTTAAGCTGa  >  1:6952739/1‑36 (MQ=25)
                                 |                                
GCGCGGCGGCCCCCAGTGGAGGCCTGCTGATCAACCTCTTTGATCAGCCCGCGTTCGATAAGCTGA  >  REL606/683463‑683528

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG