Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,131,564 T→G 5.5% intergenic (‑223/+77) 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,5640T→G 5.5% 3.6 ‑2.7 ‑1.155intergenic (‑223/+77)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (20/32):  total (23/32)
Fisher's exact test strand distribution p-value = 6.75e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.85e-02

ATATGTGCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGG  >  REL606/2131529‑2131599
                                   |                                   
atatGTGCCAGCGCTGCGTAGGCATGATAAGACGCt                                     <  1:2273578/36‑1 (MQ=255)
atatGTGCCAGCGCTGATTAGGCATGATAAGACGCt                                     <  1:7375820/36‑1 (MQ=25)
  atGTGCCAGCGCTGCGTAGGCATGATAAGACGCTTc                                   <  1:6228169/36‑1 (MQ=255)
  atGTGCCAGCGCTGCGTAGGCATGATAAGACGCTTc                                   <  1:1639654/36‑1 (MQ=255)
  atGTGCCAGCGCTGCGTAGGCATGATAAGACGCTTc                                   <  1:2048155/36‑1 (MQ=255)
     tgCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGc                                <  1:430607/36‑1 (MQ=255)
     tgCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGc                                <  1:3146535/36‑1 (MQ=255)
     tgCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGc                                <  1:2516271/36‑1 (MQ=255)
     tgCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGc                                >  1:2175768/1‑36 (MQ=255)
     tgCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGc                                <  1:1067985/36‑1 (MQ=255)
     tgCCAGCGCTGCGTACGCATGATAAGACGCTTCAGc                                >  1:4990030/1‑36 (MQ=255)
      gCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCg                               >  1:6084697/1‑36 (MQ=255)
      gCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCg                               >  1:3644325/1‑36 (MQ=255)
      gCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCg                               <  1:2950052/36‑1 (MQ=255)
      gCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCg                               <  1:6256898/36‑1 (MQ=255)
       ccAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGt                              <  1:6139036/36‑1 (MQ=255)
       ccAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGt                              <  1:7377388/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:1205812/1‑36 (MQ=255)
               gCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                      >  1:6638387/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:6671208/36‑1 (MQ=255)
                     gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                <  1:3764386/36‑1 (MQ=255)
                      cATGATAAGATGCNTCAGCGTCGCATCAGGCATAgg               >  1:4428406/1‑36 (MQ=255)
                        tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt             <  1:1843797/36‑1 (MQ=255)
                        tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt             <  1:6846999/36‑1 (MQ=255)
                        tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt             <  1:6604558/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)
                                  cTTCAGCTTCGCATCAGGCATAGGTTGCCGGATgcg   <  1:4814064/36‑1 (MQ=255)
                                  cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg   >  1:6896467/1‑36 (MQ=255)
                                  cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg   <  1:5941730/36‑1 (MQ=255)
                                   ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg  >  1:5149809/1‑36 (MQ=255)
                                   ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg  >  1:797325/1‑36 (MQ=255)
                                   ttCAGCGTCGCATCAGGCAGAGGTTGCCGGATgcgg  >  1:384614/1‑36 (MQ=255)
                                   |                                   
ATATGTGCCAGCGCTGCGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGG  >  REL606/2131529‑2131599

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG