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.96.8% 141.2 / ‑1.0 31intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (27/3);  total (28/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACAT  >  REL606/2131383‑2131530
                                                                           |                                                                        
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                          <  2:976684/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                          <  2:1977957/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                          <  2:35456/76‑1 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                          >  1:1829938/1‑76 (MQ=255)
aaatgaTTTGTTTTTAAATGTCGAGATAACGAGCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                          <  2:35478/76‑1 (MQ=255)
   tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                       >  1:1766086/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                      >  2:1864134/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                    <  1:1814005/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                    <  1:635268/76‑1 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                    <  1:1517848/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                   >  1:1560134/1‑76 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                   >  2:1633776/1‑76 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                  <  1:1363365/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                  <  1:543265/76‑2 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                 >  1:1570010/1‑74 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                >  1:1519739/1‑73 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                            >  1:965980/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                            >  1:423002/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                           >  1:1411481/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                         >  2:29379/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                      >  1:533516/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                    >  1:1748196/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                   >  1:376096/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                >  2:89286/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                >  2:737642/1‑76 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                >  2:1266988/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                               >  1:1550227/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                              <  2:325780/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                              >  1:1231481/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                              >  1:393786/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAAGCAAAAAACACGGCCTCCCGGAGGGAAGcc                                              <  2:110065/76‑1 (MQ=255)
                             cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                            >  2:1673241/1‑76 (MQ=255)
                              gATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                            >  2:123187/1‑76 (MQ=255)
                               aTCAAAAAAACGAAAATATGAAAATTAATAAAGGAATAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaa                                          <  2:1829938/76‑1 (MQ=255)
                                  acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                       >  2:1118288/1‑76 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                      >  2:1290532/1‑76 (MQ=255)
                                   aaaaaaCGACAATATGAAAATTAATAGAGGAATAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                      <  2:1460840/75‑1 (MQ=255)
                                   aaaaaaCGAAAATATAAAAATAAATAGAGGAAAAAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                      <  2:522096/75‑1 (MQ=255)
                                      aaaCGACAATAGAAAAATTAAAAAAGGAAGAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTg                                   <  2:220929/76‑1 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                >  1:336378/1‑76 (MQ=255)
                                          gACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGcgcg                               >  2:932550/1‑76 (MQ=255)
                                          gACAATATAAAAATTAATAGAGAAATAAAAGGCAAAAAAACCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGcgcg                               <  1:280091/76‑1 (MQ=255)
                                             aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAGCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                            <  2:1837173/76‑1 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                        >  1:543080/1‑76 (MQ=255)
                                                  gAAAATTAATAGAGGAATAAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                       <  2:1307930/76‑1 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                    >  1:249399/1‑76 (MQ=255)
                                                     aaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                    >  1:1479106/1‑76 (MQ=255)
                                                          aTAGAGGAATGAAAGAAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt               <  2:1677946/76‑1 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg              >  1:721136/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTCATGCAGGATg              >  1:261886/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc             >  1:1453285/1‑76 (MQ=255)
                                                            agagAAATAAAAGACAAAAAAAGCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc             <  2:752917/76‑1 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt        >  1:1086394/1‑76 (MQ=255)
                                                                 aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt        >  2:277230/1‑76 (MQ=255)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa      >  2:189174/1‑76 (MQ=255)
                                                                     aaaGGGCAAAAAGACCGGCCTCCGGGGGGAAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa     <  1:638300/76‑1 (MQ=25)
                                                                      aaagCAAAAAAAAACGGCCGCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc    <  1:1961973/73‑1 (MQ=255)
                                                                        aagaaaaaaacaCGGGCTCCCGGGGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACat  <  1:144440/73‑1 (MQ=37)
                                                                           |                                                                        
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACAT  >  REL606/2131383‑2131530

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