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.98.0% 490.7 / 1.7 100intergenic (‑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 . (80/18);  total (81/19)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533
                                                                           |                                                                           
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             >  1:5546948/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             >  2:6398698/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                             <  2:1797400/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  2:7445370/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  2:3734598/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  1:7295228/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  2:2417869/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                            <  1:7746895/76‑1 (MQ=255)
  atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                           >  1:6583259/1‑76 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                          >  1:2960012/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                         >  1:5493378/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        <  1:1244467/76‑1 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        <  1:1345228/76‑1 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                        >  1:1954669/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  2:7983818/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  1:6230495/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      <  2:2374742/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      <  2:3141724/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                      >  2:5828629/1‑76 (MQ=255)
       gtGTTTTTAAATGGCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGTAATGAAAGCCaaaaaaaa                                                                      <  2:2374741/75‑1 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     <  2:2039798/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                     >  1:4558293/1‑75 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                    <  2:1546543/76‑3 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                   <  1:5856132/76‑4 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                   >  2:6889636/1‑73 (MQ=255)
            tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                                >  1:6857563/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  1:2807369/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  2:1668394/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  1:172008/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                               >  1:4908882/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                              >  2:1449018/1‑76 (MQ=255)
               aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                             <  2:4994443/76‑1 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGcccccc                                                             >  1:7072963/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  2:3025582/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  2:6945092/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                            >  2:5772370/1‑76 (MQ=255)
                  tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                          >  2:7433925/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                         >  1:4792415/1‑76 (MQ=255)
                    tCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAg                                                        >  1:745289/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:2506813/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       >  1:950930/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       <  2:3759118/76‑1 (MQ=255)
                     agagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                       <  1:4642857/75‑1 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      >  1:6975521/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      <  2:6075012/76‑1 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      >  2:3542779/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                      >  1:6272373/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                   <  1:5102032/76‑1 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                  >  1:1931989/1‑76 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                  <  2:4733610/74‑1 (MQ=255)
                            aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                  >  1:4733610/1‑74 (MQ=255)
                               aTCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaa                                             >  1:6966912/1‑76 (MQ=255)
                                 cacaAAAACAAAAATAAAAAAATAAATAGAGAAATGAAAGGCAAAAAAACCGGCCCCCCGGAGGGAAGCCGcaaca                                           <  1:1126697/76‑1 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                         >  1:5620350/1‑76 (MQ=255)
                                   aaaaaaCAAAAATATAAAAATAAAAAGAGGAAAAAAAGGAAAAAAAACCGCCCCCCGGGGGGAAAGCCGCAACaaa                                         <  2:2164016/75‑1 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                        >  1:7651717/1‑76 (MQ=255)
                                      aaaCGAAAATATGAAAATTAATAGAGGAAAGAAAGGAAAAAAAAACGGCCCCCCGGGGGAAAGCCGCAACAAAGTg                                      <  2:5764187/76‑1 (MQ=255)
                                      aaaCAAAAATATAAAAATAAAAAGAGAAATAAAAGAAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTg                                      <  2:62628/76‑1 (MQ=255)
                                        aCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcg                                    >  2:2339482/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                   >  2:4109728/1‑76 (MQ=255)
                                          aaaaaTGTGAAAATAAATAGGGGAATAAAAGAAAAAAAAAACGGCCCCCGGGAGGAAAGCCGCAACAAAGTGcgcg                                  <  2:2684856/73‑1 (MQ=255)
                                           aaaaTATAAAAATTAATAGAGGAATAAAAGGAAAAAAAAACGCCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGt                                 <  2:7383493/74‑1 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:7773512/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  1:2384172/1‑76 (MQ=255)
                                           aCAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGt                                 >  2:5503203/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  2:5255480/1‑76 (MQ=255)
                                            cAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                                >  2:2463913/1‑76 (MQ=255)
                                             aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                               >  1:6979660/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                             >  2:3866495/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATAAAAGGAAAAAAACACGGCCCCCCGGGGGGAAGCCGCAACAAAGTGCGCGTGGtt                             <  2:6169714/76‑1 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAAACGCAACAAAGAGCGCGTGGgttt                           >  2:6360565/1‑76 (MQ=255)
                                                  gAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                          >  1:1213529/1‑76 (MQ=255)
                                                  gAAAATTAATAGAGAAAGAAAAGGCAAAAAAACCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                          <  1:3444704/76‑1 (MQ=255)
                                                   aaaaTAAATAGAGAAATGAAAGGAAAAAAAACCGCCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                         <  1:2071374/76‑1 (MQ=255)
                                                    aaaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATg                        >  1:6740969/1‑76 (MQ=255)
                                                    aaaTTAAGAGAGGAATAAAAGGAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATg                        <  2:6575499/76‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                       >  2:3398938/1‑76 (MQ=255)
                                                      ataaaTAGGGGAATGAAAGGCAAAAAAACCGGCCCCCGGGGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      <  2:3903795/73‑1 (MQ=255)
                                                      aTTAAAAGAGAAATGAAAGAAAAAAAAACCGGCCCCCGGGAGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                      <  2:2870989/76‑1 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:7156159/1‑76 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:3196320/1‑76 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                     >  2:7698170/1‑76 (MQ=255)
                                                          aTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                  >  1:4654424/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                 >  2:2849325/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  1:180323/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  1:438266/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc                >  2:3576006/1‑76 (MQ=255)
                                                              agaaaTGAAAGGAAAAAAAACCGGCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg              <  1:6396547/73‑1 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg              >  1:1669277/1‑76 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg              >  2:7214168/1‑76 (MQ=255)
                                                                gAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            >  1:5317823/1‑76 (MQ=255)
                                                                aaaTAAAAGGAAAAAAAACCGCCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcg            <  1:2819835/75‑1 (MQ=25)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt           >  2:6628453/1‑76 (MQ=255)
                                                                  aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGAAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa          >  1:3964169/1‑76 (MQ=39)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa         >  1:136441/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  2:3058660/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc       >  1:6106256/1‑76 (MQ=255)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAACCa      >  1:1964247/1‑76 (MQ=38)
                                                                        agaaaaaaaaagCGGCCCCCGGGGGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACat     <  1:1127931/73‑1 (MQ=37)
                                                                           aaaaaaaacCGGCCTCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAtg  <  1:1566617/76‑1 (MQ=255)
                                                                           |                                                                           
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATATG  >  REL606/2131383‑2131533

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