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.97.8% 679.6 / 3.9 139intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (116/20);  total (117/22)
Fisher's exact test for biased strand distribution p-value = 6.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533
                                                                           |                                                                           
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             <  1:3124263/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             <  1:4813178/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             <  1:7408188/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  1:3869204/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          <  1:6154569/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          <  1:3823944/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          <  1:6604071/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          <  1:7154279/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          <  1:1506521/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          >  2:3586875/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         <  1:5559122/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         <  1:1459075/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         >  1:4631438/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         >  2:1826887/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         >  2:2037643/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        >  1:220910/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        >  1:2555497/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                       <  2:7990122/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                       <  1:4130769/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  1:3945873/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  2:5995307/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      <  1:298468/76‑1 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     >  2:4796242/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     >  1:2382529/1‑75 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                   >  2:4291654/1‑73 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  2:5442048/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:4906302/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:5349833/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  2:7456996/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  1:797558/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              <  1:6370531/76‑1 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  2:728789/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  2:2344997/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  1:6224359/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  2:2074264/1‑76 (MQ=255)
                  tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCGTCCCgg                                                          <  1:1377512/76‑1 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                         >  2:1347675/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                         >  1:791510/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                         <  1:7077575/76‑1 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:7635481/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  2:3849581/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:39221/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑‑AAAAAAACGGCCTCCCGGAggg                                                      >  1:6823161/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAAGCAAAAAAAACGGCCTCCCGGAgg                                                       <  2:7152279/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  2:6595820/1‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    <  2:3978287/76‑1 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                    <  1:6318194/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 >  1:1923307/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 >  2:1166077/1‑76 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                                <  2:1923366/76‑1 (MQ=255)
                              gATCAAAAAAACGACAATATAAAAATTAATAGAGAAATGAAAGGAAAAAAAAACGGCCCCCCGGAGGAAAGCCGca                                              <  2:5341696/76‑1 (MQ=255)
                               aTCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGcaa                                             <  1:4069168/76‑1 (MQ=255)
                               aTCAAAAAAACGACAATATGAAAATTAATAGAGGAATAAAAGGAAAAAAAACCGGCCCCCGGGAGGGAAGCCGcaa                                             <  2:3305571/76‑1 (MQ=255)
                                tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                            >  1:4412991/1‑76 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  1:1026106/1‑76 (MQ=255)
                                  aaaaaaaCGCCAATATAAAAATAAATGGGCGAAGAAAAGCAAAAAAAACGGGCCCCCGGGAGGAAAGCCGcaacaa                                          <  2:7192810/74‑1 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                         >  2:6203954/1‑76 (MQ=255)
                                   aaaaaaCGAAAATAGAAAAATTAATAGAGGAATGAAGGAAAAAAAAAGCGCCCCCCGGGGGGAAAGCCGCAACaaa                                         <  2:4991245/75‑1 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                        >  1:6030078/1‑76 (MQ=255)
                                         caaCAATATGAAAATAAATAGAGGAATAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   <  1:2909693/74‑1 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  1:6456740/1‑76 (MQ=255)
                                         cGAAAATAAAAAAATTAAAAGAGAAATAAAAGGAAAAAAAACCGGCCCCCCGGGGGAAAGCCGCAACAAAGTgcgc                                   <  2:5741993/76‑1 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  1:3194306/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  1:7608281/1‑76 (MQ=255)
                                            aaaTATAAAAATAAATAGAGAAATAAAAGGAAAAAAAACCGCCCCCCGGGGGGAAAGCCGCAACAAAGTGCGCGTg                                <  2:2224599/75‑1 (MQ=255)
                                             aaTATAAAAATTAATAGAGGAATGAAAGGCAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTgg                               <  2:7807206/76‑1 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                             >  2:6299392/1‑76 (MQ=255)
                                                 taaaaaTTAATAGAGGAATGAAAGGAAAAAAAAGCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGtttt                           <  2:1410432/74‑1 (MQ=255)
                                                   aaaaTTAAGAGAGAAATGAAAGGAAAAAAAACCGCCCCCGGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  2:3473115/76‑1 (MQ=255)
                                                   aaaaTTAAGAGAGAAATAAAAGGAAAAAAAAACGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  2:1538972/76‑1 (MQ=255)
                                                      aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      >  1:327507/1‑76 (MQ=255)
                                                      aTTAAAAGAAGAATGAAAGAAAAAAAAAGTGACCTCCGGGGGAGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      <  2:2568586/76‑1 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGCTTTATGCAg                     >  1:4539627/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:1598462/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:1723927/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  1:5483052/1‑76 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                   >  1:125657/1‑76 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                   >  2:7703822/1‑76 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAACTGCGCGTGGTTTTATGCAGGa                   >  1:5815840/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 >  1:6703566/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 >  2:5509170/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 >  1:4724313/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  1:1376679/1‑76 (MQ=255)
                                                               ggaaggaaAGGAAAAAAAACCGGCCCCCCGGGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc             <  2:7894325/76‑1 (MQ=37)
                                                               ggAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc             >  1:6285581/1‑76 (MQ=255)
                                                                gAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            >  1:2951358/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  2:2713721/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  2:3005397/1‑76 (MQ=255)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          >  2:6740409/1‑76 (MQ=255)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa         >  1:753689/1‑76 (MQ=255)
                                                                    gAAAGGCAAAAAACACGCCCTCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        <  1:1969842/76‑1 (MQ=37)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        >  2:5991447/1‑76 (MQ=255)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        >  1:3005580/1‑76 (MQ=255)
                                                                    gAAAGGAAAAAAAACCGGCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        <  1:5249342/76‑1 (MQ=37)
                                                                     aaaGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       <  1:3744136/76‑1 (MQ=38)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACa      >  2:4171210/1‑76 (MQ=255)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAACCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCAGCGTAAACa      >  2:1672344/1‑76 (MQ=38)
                                                                      aaGGAAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       <  2:5556120/76‑1 (MQ=37)
                                                                         agaaaaaaaaaCGGCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACata    <  1:5541791/74‑1 (MQ=38)
                                                                           aaaaaaaacCGGCCCCCCGGGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAtg  <  2:5954158/76‑1 (MQ=37)
                                                                           |                                                                           
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533

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