Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,131,584 A→C 4.9% intergenic (‑243/+57) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,131,5840A→C 4.9% 2.4 ‑2.3 ‑0.862intergenic (‑243/+57)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (26/32):  total (29/33)
Fisher's exact test strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.76e-01

GCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCC  >  REL606/2131550‑2131619
                                  |                                   
gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                                    <  1:6671208/36‑1 (MQ=255)
gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                                    <  1:3764386/36‑1 (MQ=255)
 cATGATAAGATGCNTCAGCGTCGCATCAGGCATAgg                                   >  1:4428406/1‑36 (MQ=255)
   tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                 <  1:4369716/36‑1 (MQ=255)
   tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                 <  1:6846999/36‑1 (MQ=255)
   tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                 <  1:1843797/36‑1 (MQ=255)
   tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                 <  1:6604558/36‑1 (MQ=255)
   tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                 <  1:6123701/36‑1 (MQ=255)
     aTAAGACGCTTCTGCGTCGCATCAGGCATAGGTTGc                               <  1:1498379/36‑1 (MQ=255)
     aTAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGc                               >  1:5811566/1‑36 (MQ=255)
       aaGACGCTTCAGCGTCGCATCAGGCATAGGTTGACg                             >  1:2283650/1‑36 (MQ=255)
        aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                            <  1:4356025/36‑1 (MQ=255)
        aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                            >  1:4516250/1‑36 (MQ=255)
        aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                            <  1:6670155/36‑1 (MQ=255)
         gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGa                           <  1:2109609/36‑1 (MQ=255)
          tcGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAt                          >  1:1664109/2‑36 (MQ=17)
          aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                          >  1:2984003/1‑36 (MQ=255)
          aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                          <  1:3955704/36‑1 (MQ=255)
           cGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATg                         <  1:6665006/36‑1 (MQ=255)
           cGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAta                         >  1:6432677/1‑35 (MQ=17)
           cGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAta                         >  1:4501669/1‑35 (MQ=17)
             cctCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg                       <  1:5001276/34‑1 (MQ=255)
             cTTCAGCTTCGCATCAGGCATAGGTTGCCGGATgcg                       <  1:4814064/36‑1 (MQ=255)
             cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg                       <  1:5941730/36‑1 (MQ=255)
             cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg                       >  1:6896467/1‑36 (MQ=255)
              ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg                      >  1:5149809/1‑36 (MQ=255)
              ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg                      >  1:797325/1‑36 (MQ=255)
              ttCAGCGTCGCATCAGGCAGAGGTTGCCGGATgcgg                      >  1:384614/1‑36 (MQ=255)
              tgcAGCGTCGCATCAGGCATGGGGTGCCGGATgcgg                      <  1:3180988/34‑1 (MQ=11)
               tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc                     <  1:2509702/36‑1 (MQ=255)
               tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc                     >  1:59192/1‑36 (MQ=255)
               tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc                     <  1:5766060/36‑1 (MQ=255)
               tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc                     >  1:1734859/1‑36 (MQ=255)
               tCAGCGTCGCATCAGGCATAGGTTGCCGGATgaggc                     >  1:5710470/1‑36 (MQ=255)
                cAGCGTCGCATCAGGCATAGGTTGCCGGATgcggcg                    >  1:6654778/1‑36 (MQ=255)
                 aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt                   >  1:3998741/1‑36 (MQ=40)
                 aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt                   <  1:4773041/36‑1 (MQ=40)
                 aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt                   <  1:1871282/36‑1 (MQ=40)
                 aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt                   <  1:1609486/36‑1 (MQ=40)
                  gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa                  <  1:4648616/36‑1 (MQ=255)
                  gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa                  <  1:4823808/36‑1 (MQ=255)
                  gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa                  <  1:5713506/36‑1 (MQ=255)
                    gTCGCATCAGGCATAGGTTGCCGGATGCGGCGTaaa                <  1:54143/36‑1 (MQ=255)
                    gTCGCATCAGGCATAGGTTGCCGGATGCGGCGTaaa                >  1:4946421/1‑36 (MQ=255)
                      cGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACg              >  1:6434255/1‑36 (MQ=255)
                      cGCATCAGGCATACGTTGCCGGATGCGGCGTAAACg              >  1:3250586/1‑36 (MQ=255)
                        cATCAGGCATAGGTTGCCGGATGCGGCGTAAACGcc            <  1:712012/36‑1 (MQ=255)
                        cATCAGGCATAGGTTGCCGGATGCGGCGTAAACGcc            >  1:4191562/1‑36 (MQ=255)
                         aTCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCt           <  1:2811470/36‑1 (MQ=255)
                         aTCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCt           <  1:6633207/36‑1 (MQ=255)
                          tCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCtt          <  1:4205536/36‑1 (MQ=255)
                           cAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTa         >  1:5707181/1‑36 (MQ=255)
                            aGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTAt        >  1:4776511/1‑36 (MQ=255)
                            aGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTAt        >  1:4603756/1‑36 (MQ=255)
                              gCATAGGTTGCCGGATGCGGCGTAAACGCCTTATcc      <  1:3357031/36‑1 (MQ=40)
                                aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt    <  1:2048824/36‑1 (MQ=38)
                                aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt    >  1:5736483/1‑36 (MQ=38)
                                aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt    >  1:3424653/1‑36 (MQ=38)
                                aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt    >  1:3620755/1‑36 (MQ=38)
                                aTAGGTTGCCGGATGCGGCGTAAACCCCTTATCCGt    >  1:1481627/1‑36 (MQ=255)
                                aTAGGTTGCCGGAGGCGGTGTAAACGCCTTTTCCga    >  1:2848833/1‑35 (MQ=11)
                                 tAGATTGCCGGATGCGGCGTAAACGCCTTATCCGTc   <  1:32392/36‑1 (MQ=12)
                                  aGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTcc  >  1:1619652/1‑36 (MQ=35)
                                  |                                   
GCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCC  >  REL606/2131550‑2131619

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG