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.9% 441.4 / 1.4 93intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (75/16);  total (75/18)
Fisher's exact test for biased strand distribution p-value = 3.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533
                                                                           |                                                                             
aaatgaTTTGTTTTTAAGTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                               >  1:1784747/1‑76 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTTATAGAGGAATGAAAGGCaa                                                                              <  1:2163733/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                              <  2:5072740/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                              <  1:4740897/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                              <  1:3278072/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                              >  2:5581461/1‑76 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                             <  2:1201552/76‑1 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                             >  2:5092383/1‑76 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                            <  2:352664/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                            <  2:1901024/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                            >  1:2716175/1‑76 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                            <  1:2083164/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                           >  2:3110692/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                           >  1:5792957/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                           >  2:1905181/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                          >  2:245948/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                         <  1:4796883/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                         <  1:4929489/76‑1 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                       <  2:4229452/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                       <  2:5695078/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                       >  1:3377307/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                       >  2:421246/1‑75 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACGAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                      >  2:2674409/1‑74 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                      >  2:5010337/1‑74 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                     >  2:3679274/1‑73 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                     <  1:4135106/76‑4 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                 >  1:200577/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                 >  1:3457728/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                 >  2:1593784/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                                 <  1:1945989/76‑1 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                               >  1:5669226/1‑76 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                            <  2:3851802/76‑1 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                            <  2:4526010/76‑1 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                          >  2:4795479/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                          >  1:5651433/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                          >  1:4720434/1‑76 (MQ=255)
                    tCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAg                                                         <  2:3288682/76‑1 (MQ=255)
                    tCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAg                                                         >  1:5338808/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAGGGCAAAAAAAACGGCCTCCCGGAgg                                                        <  1:830955/76‑1 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                        >  2:5969378/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                        >  1:749056/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  1:4248880/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:2064978/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     >  1:403819/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:5739448/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  1:746857/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                     <  2:2440233/76‑1 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAGACGGCCTCCCGGAGGGaa                                                    >  2:791988/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   >  1:2918861/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   >  2:2194767/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   <  2:2530541/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   >  1:2942713/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   >  2:3568451/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                  >  1:1959838/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 <  2:4924991/76‑1 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                 <  1:4843025/76‑1 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                                <  1:3784284/76‑1 (MQ=255)
                             cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                               >  1:5820120/1‑76 (MQ=255)
                              gATCAAAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCCCCGGGAGGAAAGCCGca                                              <  2:2034220/76‑1 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  1:2132713/1‑76 (MQ=255)
                                 cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                           >  1:2013720/1‑76 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                          >  1:3775566/1‑76 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                          >  2:2035468/1‑76 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                         >  2:2867700/1‑76 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                        >  1:724126/1‑76 (MQ=255)
                                       aaCGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgc                                     <  1:4647873/76‑1 (MQ=255)
                                       aaCGACAATATGAAAATTAAAAGAGGAATGAAAGGCAAAAAAAACGGCCACCCGGAGGAAAGCCGCAACAAAGTgc                                     <  1:3142404/76‑1 (MQ=255)
                                        aCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcg                                    >  1:2902648/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  1:3900902/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  2:5536409/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  2:4664818/1‑76 (MQ=255)
                                            aaaTATGAAAATTAAGAGAGGAAGAAAAGGAAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTg                                <  2:3333387/75‑1 (MQ=255)
                                            aaaTATGAAAATAAAGAGAGAAATAAAAGGAAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                <  1:3480491/75‑1 (MQ=255)
                                              atatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGCGGGAAGCCGCAACAAAGTGCGCGTGGt                              >  1:3816643/1‑76 (MQ=255)
                                              atatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGt                              >  2:5322279/1‑76 (MQ=255)
                                                atGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                            >  1:4588848/1‑76 (MQ=255)
                                                aaaaaaaTAAAAAGAGAAAAGAAAGGAAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGttt                            <  2:954993/73‑1 (MQ=255)
                                                 tGAAAATAAATGGAGGAAGGAAAGAAAAAAAAACCGGCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGtttt                           <  2:3882094/76‑1 (MQ=255)
                                                  gAAAATTAATAGAGGAATGAAAGGCAAAAAAAACAGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                          >  1:513213/1‑76 (MQ=255)
                                                   aaaaTAAATAGAGAAATGAAACCAAAAAAAACCGCCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  2:4530326/76‑1 (MQ=255)
                                                    aaaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATg                        >  1:5471267/1‑76 (MQ=255)
                                                     aagtAATAGAGAAATGAAAGACAAAAAAACCGCCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       <  1:2610732/73‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       >  1:5232294/1‑76 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:2879181/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                    >  2:2998895/1‑76 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                   >  1:1121557/1‑76 (MQ=255)
                                                          aTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                  >  1:2817602/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGGGCGTGGTTTTATGCAGGATgc                >  2:2939254/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  2:288153/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTCTATGCAGAATgc                >  2:1424735/1‑76 (MQ=255)
                                                              agaaaTGAAAGG‑‑AAAAAAACCGGCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc             <  2:2882962/73‑1 (MQ=37)
                                                                gAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGTAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            >  2:5379586/1‑76 (MQ=39)
                                                                aaaTAAAAGGAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            <  1:4954947/75‑1 (MQ=25)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa         >  1:1004245/1‑76 (MQ=255)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa        >  1:4851261/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  1:1322131/1‑76 (MQ=255)
                                                                        aggaaaaaaaagCGGCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACat     <  1:1864991/73‑1 (MQ=37)
                                                                          aaaaaaaaaaGCGCCCCCCGGAGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAtata    <  1:4757304/75‑1 (MQ=25)
                                                                           aaaaaaaagCGCCCACCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAtg  <  2:338614/76‑1 (MQ=37)
                                                                           aaaaaaaaaCCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACatat   <  2:4847671/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 < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: