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.91.3% 87.8 / 2.9 23intergenic (‑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 . (15/6);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACA  >  REL606/2131383‑2131529
                                                                           |                                                                       
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                         <  2:252085/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                     >  1:345373/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                    >  2:910146/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                    >  1:452014/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                   <  2:452014/76‑1 (MQ=255)
       ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa                                                                  >  2:283135/1‑76 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                 >  2:110201/1‑75 (MQ=255)
          tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                               >  1:1126454/1‑73 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGcc                                                            >  1:85126/1‑76 (MQ=255)
              tAAATGTCGAGAGAACGATCACAAAAACGACAATATGAAAATTAATAGAGGNATGAAAGGCAAAAAAAACGGCCTc                                                          <  2:207061/76‑1 (MQ=255)
                aaTGTCGAGAAAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                        >  1:162301/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                   >  1:735973/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGAAAATATGAAAATTAATAGAGGAATGAAAGGCATAAAAAACGGCCTCCCGGAgg                                                   >  1:865193/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                 >  1:822131/1‑76 (MQ=255)
                          tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGc                                              >  1:485472/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCGAAAAAAACGGCCTCCCGGAGGGAAGcc                                             >  1:682013/1‑76 (MQ=255)
                             cGATCACAAAANCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGc                                           >  1:175789/1‑76 (MQ=255)
                                 cacaAAAACAACAAAATAAAAATAAATAGAGAAATAAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGcaaca                                       <  2:827317/76‑1 (MQ=255)
                                   caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa                                     <  2:735973/76‑1 (MQ=255)
                                            cAATATGAAAATTAAGTGAGGAATGAAAGGCAAGTAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                            <  2:772095/76‑1 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                       <  2:442361/76‑1 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                       >  1:245131/1‑76 (MQ=255)
                                                 tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                       >  2:690572/1‑76 (MQ=255)
                                                  aaaaaTAAAAAGAGAAAGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                     <  2:916121/75‑1 (MQ=255)
                                                   aaaaTTAATAGAGGAATGAAGTGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                     <  2:969732/76‑1 (MQ=255)
                                                       ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGCTTTTATGCAg                 <  1:162773/76‑1 (MQ=255)
                                                       ttAAAAGAGAAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                 <  2:162301/76‑1 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                >  2:231604/1‑76 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGACGCAACAAAGTGCGCGTGGTTTTATGCAGGATg             >  1:921341/1‑76 (MQ=255)
                                                             gacgAATGAAAGGCAAAAAAAACGGCNTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcg           >  2:998004/4‑76 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAGGTGCGCGTGGTTTTATGCAGGATgcgg          >  1:483893/1‑76 (MQ=255)
                                                              agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg          <  2:715015/76‑1 (MQ=255)
                                                                 aaTGAAAGGAAAAACAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt       <  2:721583/76‑1 (MQ=39)
                                                                    gAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa    <  2:17599/76‑1 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc   >  1:10692/1‑76 (MQ=255)
                                                                      aaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACa  >  2:164128/1‑76 (MQ=255)
                                                                           |                                                                       
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACA  >  REL606/2131383‑2131529

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: