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

TGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATA  >  REL606/2131386‑2131531
                                                                        |                                                                         
tgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                         >  1:1776845/1‑75 (MQ=255)
 gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaa                                                                        >  2:990204/1‑75 (MQ=255)
  aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                       <  1:406019/75‑1 (MQ=255)
   tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                      <  1:1366746/75‑1 (MQ=255)
      gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                   >  1:1850649/1‑74 (MQ=255)
        ttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                                 <  1:701542/75‑4 (MQ=255)
          ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                              >  1:722715/1‑75 (MQ=255)
           tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                             >  2:1468196/1‑75 (MQ=255)
            aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                            >  1:431324/1‑75 (MQ=255)
             aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                           >  2:1446181/1‑75 (MQ=255)
              aTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                          >  1:1035660/1‑75 (MQ=255)
               tGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCg                                                         >  1:1086892/1‑75 (MQ=255)
                gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCg                                                         >  2:1503494/1‑75 (MQ=255)
                gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCgg                                                        >  1:229536/1‑75 (MQ=255)
                   gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                     >  1:492842/1‑75 (MQ=255)
                    agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                    >  2:122406/1‑75 (MQ=255)
                      aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                  >  1:563685/1‑75 (MQ=255)
                        aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                                <  2:613459/75‑1 (MQ=255)
                          cGATCAAAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                              <  2:452143/75‑1 (MQ=255)
                            aTCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                             >  1:422344/1‑75 (MQ=255)
                              cacaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                          <  2:253744/75‑1 (MQ=255)
                              cacaAAAACGACAATATGAAAATTAATAGAGGAATAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaac                                          <  1:2025937/75‑1 (MQ=255)
                               acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaaca                                         >  2:1193237/1‑75 (MQ=255)
                                  aaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                      <  2:461851/75‑1 (MQ=255)
                                   aaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                      >  2:1649575/1‑75 (MQ=255)
                                     aCGACAAAATGAAAATAAATAGAGAAATAAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTgc                                   <  2:1192022/75‑1 (MQ=255)
                                       gACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                 >  1:1409924/1‑75 (MQ=255)
                                        aCAAAATAAAAATAAATAGAGAAATGAAAGCAAAAAAAACCGCCCCCCCGGAGGGAAGCCGCAACAAAGTGcgcg                                <  2:476574/75‑1 (MQ=255)
                                          aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                              >  1:526138/1‑75 (MQ=255)
                                           atatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGCgg                             >  1:673128/1‑75 (MQ=255)
                                             ataaaaaTTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                           <  2:1657828/72‑1 (MQ=255)
                                              tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                          >  1:348140/1‑75 (MQ=255)
                                                aaaaTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                        >  1:900357/1‑75 (MQ=255)
                                                   aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                     >  2:1277294/1‑75 (MQ=255)
                                                     tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                   <  2:526138/75‑1 (MQ=255)
                                                     tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                   >  1:517183/1‑75 (MQ=255)
                                                      aaTAGAGGAATGAAAGGCAAAAAAAAAGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  >  1:1306114/1‑75 (MQ=255)
                                                      aaTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                  >  1:1019916/1‑75 (MQ=255)
                                                       aTAGAGGAATGAAAGGCAAAAAAACCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                 <  2:1789666/75‑1 (MQ=255)
                                                        tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGag                >  2:689013/1‑74 (MQ=255)
                                                        tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt                >  2:2140518/1‑75 (MQ=255)
                                                            ggAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcg             >  2:1306497/1‑75 (MQ=255)
                                                            ggAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg            >  1:1909490/1‑75 (MQ=255)
                                                               aTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt         >  2:702489/1‑75 (MQ=255)
                                                                agAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTa        <  1:1961678/74‑1 (MQ=255)
                                                                   aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc     >  2:1727439/1‑75 (MQ=255)
                                                                       gaaaaaaaacCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACata  <  1:736025/74‑1 (MQ=255)
                                                                        |                                                                         
TGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATA  >  REL606/2131386‑2131531

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