Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,131,466 (A)9→8 intergenic (‑125/+175) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,131,4580A.99.0% 522.6 / ‑3.0 103intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (84/18);  total (84/19)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533
                                                                           |                                                                            
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                              >  1:2242656/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                              >  1:6888735/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                          >  1:2371875/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                          >  1:416419/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                          <  2:5581297/76‑1 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                         <  2:2470524/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                        <  2:5833247/76‑1 (MQ=255)
        tGTTTTTAAATGTCGTGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                      >  1:153041/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                      >  1:5821644/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                      <  2:5386175/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                      >  2:3219438/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                      >  1:5701496/1‑75 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                     >  2:3478112/1‑74 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                     >  1:3505863/1‑74 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                 >  2:3090890/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                 >  1:150697/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                 >  1:5649532/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                >  1:6067660/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                >  2:4000246/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                >  1:3375804/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                >  1:3313853/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                >  2:488618/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                               >  2:5330975/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                               <  2:5503404/76‑1 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                               >  2:355597/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                              <  2:6489362/76‑1 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                              >  1:5855007/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                              >  2:6810369/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                             >  1:5555507/1‑76 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                           <  2:617978/76‑1 (MQ=255)
                  tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                          <  2:2432494/76‑1 (MQ=255)
                  cgTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                          >  2:2850491/2‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:1964909/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:4985687/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:2037111/1‑76 (MQ=255)
                     cGAGATAACGAGCACAAAAACGACAATATGCAAATTAATAGAGGAACGAAAGGCAAAAAAATCGGCCTCCCGGacg                                                       <  2:2713648/76‑3 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      <  2:6684161/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  2:5664133/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:100847/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  1:2712415/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:4045778/76‑1 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    >  2:6144784/1‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    <  1:3077192/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   <  1:1428588/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   <  1:6868676/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   >  2:6338232/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                  >  1:521702/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 >  2:4616716/1‑76 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                                >  1:1526321/1‑76 (MQ=255)
                             cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                               >  1:6552441/1‑76 (MQ=255)
                             cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                               >  1:3194286/1‑76 (MQ=255)
                             cGATCAAAAAAACGACAATATAAAAATAAATAGAGGAATAAAAGGAAAAAAAAACGGCCCCCCGGAGGAAAGCCGc                                               <  2:3035771/76‑1 (MQ=255)
                              gATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGca                                              >  1:6693797/1‑76 (MQ=255)
                                tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                            >  1:6352830/1‑76 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  2:670825/1‑76 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                          >  2:4401480/1‑76 (MQ=255)
                                   aaaaaaCGAAAATATAAAAATAAAAAGAGGAAGGAAAGGAAAAAAAACCGCCCCCCCGGGGGAAAGCCGCAACaaa                                         <  2:167103/75‑1 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCAGCAACAAAg                                        >  2:4647129/1‑76 (MQ=255)
                                     aaaaCGACAATATGAAAATTAATAGAGGAATAAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGt                                       <  1:5427186/76‑1 (MQ=255)
                                     aaaaCGACAATATGAAAATTAATAGAGAAATGAAAGGAAAAAAAACCGGCCCCCCGGAGGGAAGCCGCAACAAAGt                                       <  1:4548055/76‑1 (MQ=255)
                                     aaaaCAAAAAAATAAAAATTAAAAGAGAAATAAAAGAAAAAAAAACCGCCCCCCCGGGGGAAAGCCGCAACAAAGt                                       <  2:4859588/76‑1 (MQ=255)
                                     aaaaCAAAAAAATAAAAAACAATAGAGGATGAAAAAAAAAAAAAAGCGCCCCCCGGGAGGGAAGCCGCAACAAAGt                                       <  2:4693052/76‑1 (MQ=255)
                                      aaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTg                                      <  1:1728545/76‑1 (MQ=255)
                                        aCGACAATATGAAAATTAATAGAGGAATAAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcg                                    <  1:5006492/76‑1 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  2:2539894/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:3210635/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  2:4907203/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  2:2052155/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  1:4227476/1‑76 (MQ=255)
                                             aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                               >  1:3702220/1‑76 (MQ=255)
                                             aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                               >  2:1014086/1‑76 (MQ=255)
                                              atatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGt                              >  2:2584790/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                             >  1:6209377/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                             >  2:4186025/1‑76 (MQ=255)
                                                atGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                            >  1:1015777/1‑76 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                           >  2:1434088/1‑76 (MQ=255)
                                                 aaaaaaTAAATAGAGAAAGAAAAGAAAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGtttt                           <  2:698293/74‑1 (MQ=255)
                                                   aaaaTTAAGAGGGGAATGAAAGAAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  2:897410/76‑1 (MQ=255)
                                                    aaaTTAATAGAGAAATAAAAGGCAAAAAAAACGGCCCCCCGGGGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATg                        <  1:1602962/76‑1 (MQ=255)
                                                    aaaTAAAGAGAGGAATGAAGGAAAAAAAAAACGCCCCCCCGGGGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATg                        <  2:1494751/76‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       >  2:2114372/1‑76 (MQ=255)
                                                      aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      >  1:2237922/1‑76 (MQ=255)
                                                      aTTAAGAGAGAAAGGAAAGGAAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      <  2:5535336/76‑1 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:396707/1‑76 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  1:82480/1‑76 (MQ=255)
                                                          aaagagGAATAAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                  <  1:1313493/74‑1 (MQ=255)
                                                          aTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                  >  2:6855256/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 >  1:1004582/1‑76 (MQ=255)
                                                           aagagAAATAAAAGGAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 <  2:1090839/75‑1 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  2:1276038/1‑76 (MQ=255)
                                                             gagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcg               >  1:731577/1‑76 (MQ=255)
                                                                gAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            >  1:6336025/1‑76 (MQ=255)
                                                                gAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            >  1:1666046/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  1:4771751/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  2:1533804/1‑76 (MQ=255)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          >  2:5897029/1‑76 (MQ=255)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGGGGTTTTATGCAGGATGCGGCGTa          >  1:3991159/1‑76 (MQ=255)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa         >  1:5210230/1‑76 (MQ=255)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACa      >  2:5474878/1‑76 (MQ=255)
                                                                          gaaaaaaaaaCCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACata    <  2:5658111/75‑1 (MQ=25)
                                                                           |                                                                            
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533

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