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.95.0% 465.5 / 16.5 100intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base . (86/9);  total (89/11)
Fisher's exact test for biased strand distribution p-value = 9.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533
                                                                           |                                                                           
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            >  1:573294/1‑76 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                           >  1:2819773/1‑76 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                           <  1:3745799/76‑1 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                           >  2:4684900/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         >  2:6692901/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        >  1:770693/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                       >  1:4245810/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                       >  1:5873948/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      <  1:7016359/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  1:6265578/1‑76 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     >  2:7725814/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     >  2:3530031/1‑75 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                    <  1:2830493/76‑3 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                    <  1:5458670/76‑3 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                    >  1:1186284/1‑74 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:479536/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:6746454/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:7469987/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                <  1:1089149/76‑1 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               <  1:7256751/76‑1 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  2:225507/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  2:4848351/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  1:4877236/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  2:600351/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  1:2103776/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  1:2384923/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  1:5136031/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑‑AAAAAAACGGCCTcc                                                             >  1:4883214/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  1:1275653/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  1:594160/1‑76 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                           <  1:1640660/76‑1 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCgg                                                          >  1:3048574/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAg                                                        <  2:3048574/76‑1 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  2:7507653/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGCAAGGCAAAAAAACACGGCCTCCCGGAgg                                                       >  2:3905219/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      >  2:6381749/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      >  1:2279266/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  1:4860109/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:7536178/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  1:2538035/1‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    <  2:1310426/76‑1 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    >  2:1403624/1‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    <  2:7004920/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 <  2:657675/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 >  2:3770311/1‑76 (MQ=255)
                                tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                            >  2:6425227/1‑76 (MQ=255)
                                tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                            >  1:5248830/1‑76 (MQ=255)
                                aaacaAAAGCGACAATATGAAAATGAAGAGAGGAAAGAAAGCCAAAAAACACGCCCTCGGGGGGGAAAGCCGcaac                                            <  2:5344877/74‑1 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  1:4202817/1‑76 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  2:2225133/1‑76 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                         >  1:2464005/1‑76 (MQ=255)
                                   aaaaaaCAAAATTAAGAAATTTAATAGAGAAAAAAAGGCAAAAAAAACCGGCCCCCCGGGGGAAAGCCGCAACaaa                                         <  2:4149319/75‑1 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                        >  1:32412/1‑76 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                        >  2:2675343/1‑76 (MQ=255)
                                      aaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTg                                      >  1:4529936/1‑76 (MQ=255)
                                      aaaCGAAAATATAAAAATAAAAAGAGAAATAAAAGGAAAAAAAACCGGCCCCCCGGAGGAAAGCCGCAACAAAGTg                                      <  2:876909/76‑1 (MQ=255)
                                      aaaCAAAAAAATAAAAATAAAAAGAGGAAGAAAAGAAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTg                                      <  2:191241/76‑1 (MQ=255)
                                       aaagACAATAGGAAAAATAAAAAAGAAATAAAAGGAAAAAAAAAACGGCCTCCCGGGGGGAAGCCGCAACAAAGTg                                      <  1:5540030/73‑1 (MQ=255)
                                       aaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgc                                     >  1:3320774/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  1:450292/1‑76 (MQ=255)
                                          gACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGcgcg                                  >  1:2620442/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:5615322/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:4935513/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  1:429748/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:4215382/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  1:2574018/1‑76 (MQ=255)
                                                ataaaaaTAAAAAGAGAAATAAAAGAAAAAAAAAACGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGttt                            <  2:1364781/73‑1 (MQ=255)
                                                atGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                            >  2:3553340/1‑76 (MQ=255)
                                                   aaaaTAAATAGAGGAATAAAAGGAAAAAAAACCGCCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  1:4799888/76‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       >  1:1099234/1‑76 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       >  2:5250135/1‑76 (MQ=255)
                                                      ataaaTAAAGGAATGAAAGTCAAAAAACCCGCCCTCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      <  1:7468061/73‑1 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:3659600/1‑76 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:4091059/1‑76 (MQ=255)
                                                       ttAAAAGGGGAACAAAAGGCAAAAAAACCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     <  1:3716957/76‑1 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:2513781/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:68096/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:6857745/1‑76 (MQ=255)
                                                         aaaagagAAATGAAAGGCAAAAAAAACGGCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                   <  2:3236047/73‑1 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                   >  1:4132477/1‑76 (MQ=255)
                                                          aTAGAGTAATGAAAGGAAAAAAACCCGCCCGCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                  <  1:654231/76‑1 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAACCGGCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg              <  2:3041835/76‑1 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg              >  1:19793/1‑76 (MQ=255)
                                                               tgAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc             >  2:5855957/2‑76 (MQ=255)
                                                               ggAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc             >  2:4860834/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  1:1142302/1‑76 (MQ=255)
                                                                  aggAAAGGAAAAAAAACCGCCCTCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          <  1:954738/74‑1 (MQ=25)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          >  1:7306298/1‑76 (MQ=255)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          >  2:2472344/1‑76 (MQ=255)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        >  1:7652582/1‑76 (MQ=255)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        >  1:1351490/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  1:6853642/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  2:6013958/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  1:1611435/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  2:7531261/1‑76 (MQ=255)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACa      >  2:6335320/1‑76 (MQ=255)
                                                                        ggCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGAGGCGTAAACat     >  1:7459279/1‑76 (MQ=255)
                                                                          aaaaaaaaagCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACatat   <  1:4441039/75‑1 (MQ=37)
                                                                           aaaaaaaacCGGCCTCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAtg  <  2:4926832/76‑1 (MQ=255)
                                                                           aaaaaaaacCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAtg  <  1:6435810/76‑1 (MQ=37)
                                                                           |                                                                           
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533

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