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

AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAT  >  REL606/2131383‑2131532
                                                                           |                                                                          
aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa                                                                            <  1:334083/76‑1 (MQ=255)
 aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                           >  2:1466197/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATTAAAGGCaaaaa                                                                        <  1:56373/76‑1 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                        >  1:817382/1‑76 (MQ=255)
    gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                        >  1:936083/1‑76 (MQ=255)
     aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa                                                                       >  1:335851/1‑76 (MQ=255)
      tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa                                                                      >  1:399245/1‑76 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                    <  1:754941/76‑2 (MQ=255)
        tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                    >  1:1193738/1‑75 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                   >  1:1025385/1‑74 (MQ=255)
         gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                   <  2:400837/76‑3 (MQ=255)
           ttttAAATGTCGAGATANCGATCACAAAAACGACAATATGAAAATTAATAGATGAATGAAAGGAAACAAAAACGGc                                                                >  1:214132/1‑76 (MQ=255)
             ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                              >  1:1099278/1‑76 (MQ=255)
              tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                             >  1:119153/1‑76 (MQ=255)
                aaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTcc                                                            >  2:1254343/1‑76 (MQ=255)
                   gTCGAGATAACGATCACTAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                                        >  1:1697986/1‑76 (MQ=255)
                   gTCGAGATAACGATCACAAAAACGACAATATGAAGATTAATNGAGGAATGAAAGGCAAAANAANCGGCCTCCCGGa                                                        >  1:1396242/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                      >  2:45820/1‑76 (MQ=255)
                     cGAGATAACGATCACAAAAACGACAATATGAAAAATAATAGNGGAATGAAAGGCAAAAAANACNGCCTCCCGGAgg                                                      >  1:1406798/1‑76 (MQ=255)
                     cGAGATAACGACCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAgg                                                      >  1:98362/1‑76 (MQ=255)
                      gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                     <  2:1063147/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATCGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    >  1:1440738/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    >  2:1131878/1‑76 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:1046643/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:817381/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:971132/76‑1 (MQ=255)
                       agaTAACGATCACAAAAACGACAATATAAAAATTCATAGAGGAAAGAANCGATAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:9468/76‑1 (MQ=255)
                       acctAACGCTCACAAAAACTACACTATCAAAATTAATAGATGAATAAATGTCTAAAAAAACGGCCTCCCTTAGGGa                                                    >  2:17349/4‑76 (MQ=255)
                        gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                                   <  2:1547185/76‑1 (MQ=255)
                         aTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAg                                                  >  1:1280288/1‑76 (MQ=255)
                          tAACGTTCACAAAAACGACAATATGAAAATTAATAGAGGAACGAAAGGCAAAAANAACGGCCTGCNGGATGGAAGc                                                 <  2:1682733/76‑1 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                >  1:662720/1‑76 (MQ=255)
                           aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGcc                                                >  2:549362/1‑76 (MQ=255)
                                    aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                                       >  2:1352547/1‑76 (MQ=255)
                                       aaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgc                                    >  2:675760/1‑76 (MQ=255)
                                         cGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTgcgc                                  >  2:582512/1‑76 (MQ=255)
                                               tataAAAATTAAGAGAGGAAGAAAAGCAAAAAAAACCGCCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                            <  2:1469331/76‑1 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                            >  1:1055381/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                            >  1:1344928/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtt                            >  1:378999/1‑76 (MQ=255)
                                               tatGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGCAAGCCGCAACAAAGTGCGCGTGGtg                            >  1:263519/1‑75 (MQ=255)
                                                atNAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGAntt                           <  2:1394016/76‑1 (MQ=255)
                                                  gAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa                         >  2:1025657/1‑76 (MQ=255)
                                                   aaaaTTAATAGAGAAATGAAAGACANATAAACCGGCATCCCGGAGGGNAGCAGCAACAAAGTGCGCGTGGTTTTAt                        <  1:239872/76‑1 (MQ=255)
                                                     aaTTACTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGc                      >  1:111642/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCACAAAAAAAGGCCTCCCGGAGGGAAGCCGCAACAGAGTGCGCGTGGTTTTATGCAgg                   >  2:132540/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                   >  2:421872/1‑76 (MQ=255)
                                                        tAATAGAGGAATGAACGATAAAAAAAACGGCCTCCCGGAGGGATGCCGCAACAAGGGGCGCGTGATTTTATGCAga                   >  1:20140/1‑75 (MQ=255)
                                                           tAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATg                >  1:1301690/1‑76 (MQ=255)
                                                           aagagGAATGAAAGNCAAAAAAACCGGCCTCCCGGAGGGAAGCCGCTACAAAGTGCGCGTGGTTTTATCCNGTATg                <  2:1718657/75‑1 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc               >  2:308348/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc               >  2:1028276/1‑76 (MQ=255)
                                                            agagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgc               >  2:1530312/1‑76 (MQ=255)
                                                               ggAATGAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc            <  2:535424/76‑1 (MQ=255)
                                                                  aaaaaaGGAAAAAAAAACCGGCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt          <  1:1033644/73‑1 (MQ=38)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa        >  1:397415/1‑76 (MQ=255)
                                                                   tGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaa        >  1:528743/1‑76 (MQ=255)
                                                                     aaaGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc      >  1:611040/1‑76 (MQ=255)
                                                                          gaaaaaaaaaCGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACatat  <  2:430026/75‑1 (MQ=255)
                                                                           |                                                                          
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAT  >  REL606/2131383‑2131532

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: