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.83.3% 93.2 / 19.8 30intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base . (15/10);  total (18/12)
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.87e-01

AATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAT  >  REL606/2131384‑2131532
                                                                          |                                                                          
aatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaa                                                                           <  1:736563/76‑1 (MQ=255)
   gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa                                                                        <  2:1451751/76‑1 (MQ=255)
       tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                    >  1:1442410/1‑75 (MQ=255)
       tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                    >  2:1561409/1‑75 (MQ=255)
       tGTTTCTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac                                                                    <  2:1739457/76‑2 (MQ=255)
        gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                   >  1:1557653/1‑74 (MQ=255)
        gTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacg                                                                   >  2:661633/1‑74 (MQ=255)
           tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGc                                                                >  2:794810/1‑76 (MQ=255)
           tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAACAACGGcc                                                               >  2:763206/1‑76 (MQ=255)
           tttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGcc                                                               >  2:214867/1‑76 (MQ=255)
            ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                              <  2:619086/76‑1 (MQ=255)
            ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCt                                                              <  1:1004585/76‑1 (MQ=255)
             tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                             <  1:186600/76‑1 (MQ=255)
             tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                             >  1:413569/1‑76 (MQ=255)
              aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTc                                                             <  1:341229/76‑1 (MQ=255)
              aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                            >  2:977273/1‑76 (MQ=255)
                   tCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAg                                                       >  1:757193/1‑76 (MQ=255)
                   tCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAg                                                       >  1:59171/1‑76 (MQ=255)
                     gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                     <  1:701818/76‑1 (MQ=255)
                     gagaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAggg                                                     >  1:1356667/1‑76 (MQ=255)
                      agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:680545/76‑1 (MQ=255)
                      aaaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGa                                                    <  2:252882/74‑1 (MQ=255)
                       gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGAAAAAAAAACGGCCTCCCGGTGGGaa                                                   <  2:222540/76‑1 (MQ=255)
                            cGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                               >  2:1361628/1‑76 (MQ=255)
                               tCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCCCCCGGAGGGAAGCCGcaac                                           <  2:244564/76‑1 (MQ=255)
                                 acaAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaacaa                                         >  1:412640/1‑76 (MQ=255)
                                            aaTATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTgg                              <  2:1442410/76‑1 (MQ=255)
                                               atGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                           >  1:376079/1‑76 (MQ=255)
                                               atGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                           >  1:363975/1‑76 (MQ=255)
                                               atGAAAATAAAAAGAGAAAGAAAAGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGttt                           <  2:293024/76‑1 (MQ=255)
                                                taaaaaTAAAAAGAGAAAAAAAAGGAAAAAAAACCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                          <  2:1096264/74‑1 (MQ=255)
                                                tGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt                          <  1:62984/76‑1 (MQ=255)
                                                  aaaaTAAAAAGAGAAATAAAAGAAAAAAAAAACGGCCCCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt                        <  2:1574849/76‑1 (MQ=255)
                                                     aTTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCa                     >  2:324236/1‑76 (MQ=255)
                                                      atAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg                    <  2:1665700/75‑1 (MQ=255)
                                                       tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                   >  1:1342166/1‑76 (MQ=255)
                                                       tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg                   <  2:155675/76‑1 (MQ=255)
                                                        aaTAGAGAAATAAAAGGCAAAAAAACCGGCCCCCGGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGa                  <  2:1344209/76‑1 (MQ=255)
                                                            gagGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcg              >  1:23937/1‑76 (MQ=255)
                                                             agGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg             >  2:1373150/1‑76 (MQ=255)
                                                                aaTGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggcgg          >  1:766311/1‑75 (MQ=255)
                                                                    aaaGGCACAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc      >  1:1102545/1‑76 (MQ=38)
                                                                         aaaaaaaaacCGGCCGCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACatat  <  1:1741756/75‑1 (MQ=255)
                                                                          |                                                                          
AATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACATAT  >  REL606/2131384‑2131532

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