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.1% 323.3 / 1.1 68intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (54/12);  total (55/13)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATA  >  REL606/2131383‑2131531
                                                                           |                                                                         
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                           >  2:1627661/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                           >  2:1809963/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                           <  2:2443575/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                           >  2:3582257/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                           <  2:1265098/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          <  1:2914009/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          <  1:4147555/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          >  1:3953927/1‑76 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          >  2:2334256/1‑76 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          <  1:4713476/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                          <  1:1027066/76‑1 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                         >  1:8888/1‑76 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                        <  1:3388452/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                       >  1:2191201/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAACAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                      >  1:4863977/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                      >  2:4265386/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                      >  2:490355/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTTATAGAGGAATGAAAGGCaaaaaaa                                                                     >  2:372972/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                     <  1:1205616/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                     <  2:623802/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                     >  1:3321190/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                    >  1:1496173/1‑76 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                   >  2:531527/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                   >  1:887820/1‑75 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                   <  2:458178/76‑2 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                  <  1:2175816/76‑3 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                 >  1:3899085/1‑73 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                              >  2:4883057/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                              >  2:1467670/1‑76 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                              >  1:249554/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                           >  1:1328659/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                           >  1:828588/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                           >  1:4017385/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAAGCAAAAAAAACGGCCTcc                                                           <  2:3284951/76‑1 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                         >  2:1591635/1‑76 (MQ=255)
                 aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                         <  2:3940249/76‑1 (MQ=255)
                  tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                        >  2:1800447/1‑76 (MQ=255)
                  tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                        >  2:4131386/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                       <  2:1237466/76‑1 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                     <  2:1445799/76‑1 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                     >  1:1473820/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                    <  1:3290980/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                   >  2:1323605/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                   >  2:3712411/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                   <  2:65226/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                   >  1:1517576/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                   >  1:4097306/1‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                  <  1:3712609/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                 <  2:2202510/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                 <  2:1615367/76‑1 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                >  2:2525368/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                >  1:3661061/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                >  1:2688584/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                               <  2:1884616/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                               >  2:4100899/1‑76 (MQ=255)
                             cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                             >  2:4264962/1‑76 (MQ=255)
                               aTCAAAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCCCCCGGAGGGAAGCCGcaa                                           <  2:3517488/76‑1 (MQ=255)
                                tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                          >  2:3944440/1‑76 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                        >  1:1241216/1‑76 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                        >  1:2041374/1‑76 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                       >  1:4366483/1‑76 (MQ=255)
                                      aaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGt                                     >  2:2628966/1‑76 (MQ=255)
                                      aaaCAAAAATATAAAAATTAAAAGAGAAAAAAAAGGAAAAAAAACCGCCCCCCGGGGGGAAAGCCGCAACAAAGTg                                    <  2:2673377/76‑1 (MQ=255)
                                      aaaCAAAAATATAAAAATAAAAAGGGGAATAAAAGGAAAAAAAACCGCCCCCCCGGGGGGAAGCCGCAACAAAGTg                                    <  2:440406/76‑1 (MQ=255)
                                        aCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcg                                  >  2:3372458/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGATTGACTGGCAAAAAAACACGGCCTCCCGGAGGGCAGCCGCAAAAACGTgcg                                  >  2:4810168/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                               >  2:176351/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                               <  1:4891443/76‑1 (MQ=255)
                                             aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                             >  2:1348982/1‑76 (MQ=255)
                                             aaTATAAAAATTAATAGAGGAATAAAAGCAAAAAAAACCGGCCCCCCGGGGGAAAGCCGCAACAAAGTGCGCGTgg                             <  2:3395446/76‑1 (MQ=255)
                                             aaTATAAAAATAAATAGAGAAATAAAAGGAAAAAAAACCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                             <  1:4076377/76‑1 (MQ=255)
                                              atatGAAAATAAATAGAGAAATGAAATAAAAAAAAAACCGCCCTCCCGGAGGAAAGCCGCAACAAAGTGCGCGTgg                             <  1:516569/76‑1 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                           >  1:747598/1‑76 (MQ=255)
                                                ataaaaaTTAATAGAGAAATAAAAGAAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGttt                          <  2:3044559/73‑1 (MQ=255)
                                                aaaaaaaTAAAGAGAGAAAGAAAAGAAAAAAAAACCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGttt                          <  2:1371705/73‑1 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                         >  2:2227232/1‑76 (MQ=255)
                                                 aaaaaaTTAAAGGAGAAATAAAAGAAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGtttt                         <  2:4772010/74‑1 (MQ=255)
                                                  gAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                        >  1:3833763/1‑76 (MQ=255)
                                                  gAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                        >  1:4452969/1‑76 (MQ=255)
                                                  gAAAATAAATAGAGGAATAAAAGGCAAAAAAACCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                        <  1:3281980/76‑1 (MQ=255)
                                                   aaaaTAAATAGGGAAGTAAAGGAAAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTAt                       <  2:2160868/76‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                     >  2:2032575/1‑76 (MQ=255)
                                                     aaTTAAAAGAGGAATGAAAGAAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                     <  2:2094843/76‑1 (MQ=255)
                                                      ataaaTAGAGAAATGAAAGGAAAAAAACCCGGCCTCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                    <  1:2518215/73‑1 (MQ=255)
                                                      aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                    >  1:3099766/1‑76 (MQ=255)
                                                      aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                    >  1:2840996/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAACCGCCCCCCCGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  <  1:4265386/76‑1 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  >  2:2874196/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  >  2:1570763/1‑76 (MQ=255)
                                                        tAAAAGAGGAATAAAAGGAAAAAAAACCGCCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  <  1:2791616/76‑1 (MQ=255)
                                                         aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                 >  2:1105382/1‑76 (MQ=255)
                                                            agagGAATGAAAGAAAAAAAAACCGCCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc              <  1:3283949/76‑1 (MQ=255)
                                                                gAATGAAAGAAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg          <  2:145391/76‑1 (MQ=37)
                                                                  aTGAAAGACAAAAAAACCGGCCTCCCGGGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa        <  2:1058235/76‑1 (MQ=37)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa       >  1:1160307/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc     >  2:4095960/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc     >  1:3308218/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc     >  2:2030404/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGGGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc     >  2:1036054/1‑76 (MQ=38)
                                                                         ggaaaaaaaaaCGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACata  <  1:2239134/74‑1 (MQ=37)
                                                                           |                                                                         
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATA  >  REL606/2131383‑2131531

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