Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,131,571 T→A 5.7% intergenic (‑230/+70) 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,5710T→A 5.7% 2.6 ‑3.2 ‑0.960intergenic (‑230/+70)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (24/32):  total (28/32)
Fisher's exact test strand distribution p-value = 4.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.35e-01

CCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAA  >  REL606/2131536‑2131605
                                   |                                  
ccAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGt                                    <  1:7377388/36‑1 (MQ=255)
ccAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGt                                    <  1:6139036/36‑1 (MQ=255)
ccAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGt                                    >  1:384296/1‑36 (MQ=255)
 cAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTc                                   >  1:7431136/1‑36 (MQ=255)
 cAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTc                                   <  1:6873381/36‑1 (MQ=255)
 cAGCGCTGCGTAGGCATGATAAGACGCTTCAGAGTc                                   >  1:7498595/1‑36 (MQ=255)
 cAGCGCTGCGTAGGCANGATAAGACGCTTCAGCGTc                                   <  1:2140142/36‑1 (MQ=255)
    cgcTGCGTAGGCATGATAAGACGCTTCAGCGTCGCa                                <  1:3615096/36‑1 (MQ=255)
      cTGCGTAGGCATGATAAGACGCTTCAGCGTCGCATc                              >  1:1103257/1‑36 (MQ=255)
       tGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCa                             <  1:821657/36‑1 (MQ=255)
       tGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCa                             >  1:6804786/1‑36 (MQ=255)
        gCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                            >  1:6638387/1‑36 (MQ=255)
        gCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                            >  1:1205812/1‑36 (MQ=255)
        gCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                            >  1:2372831/1‑36 (MQ=255)
        gCGTAGGCATGATAAGACGCTTCAGACACGCATccg                            >  1:2448681/1‑34 (MQ=17)
          gTAGGCATGATAAGACCATTCAGCGTCGCATAAGGc                          >  1:5749824/1‑36 (MQ=14)
           tAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                         <  1:1224943/36‑1 (MQ=38)
           tAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                         >  1:4335882/1‑36 (MQ=38)
            aGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                        <  1:692875/36‑1 (MQ=38)
              gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                      <  1:3764386/36‑1 (MQ=255)
              gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                      <  1:6671208/36‑1 (MQ=255)
               cATGATAAGATGCNTCAGCGTCGCATCAGGCATAgg                     >  1:4428406/1‑36 (MQ=255)
                 tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                   <  1:6846999/36‑1 (MQ=255)
                 tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                   <  1:6604558/36‑1 (MQ=255)
                 tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                   <  1:1843797/36‑1 (MQ=255)
                 tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                   <  1:6123701/36‑1 (MQ=255)
                 tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                   <  1:4369716/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:6670155/36‑1 (MQ=255)
                      aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg              <  1:4356025/36‑1 (MQ=255)
                      aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg              >  1:4516250/1‑36 (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:797325/1‑36 (MQ=255)
                            ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg        >  1:5149809/1‑36 (MQ=255)
                            ttCAGCGTCGCATCAGGCAGAGGTTGCCGGATgcgg        >  1:384614/1‑36 (MQ=255)
                            tgcAGCGTCGCATCAGGCATGGGGTGCCGGATgcgg        <  1:3180988/34‑1 (MQ=11)
                             tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc       <  1:5766060/36‑1 (MQ=255)
                             tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc       >  1:59192/1‑36 (MQ=255)
                             tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc       <  1:2509702/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:1871282/36‑1 (MQ=40)
                               aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt     <  1:4773041/36‑1 (MQ=40)
                               aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt     >  1:3998741/1‑36 (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:4946421/1‑36 (MQ=255)
                                  gTCGCATCAGGCATAGGTTGCCGGATGCGGCGTaaa  <  1:54143/36‑1 (MQ=255)
                                   |                                  
CCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAA  >  REL606/2131536‑2131605

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