Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,048,224 G→T 14.8% Q109K (CAA→AAA)  yihR ← predicted aldose‑1‑epimerase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,048,2240G→T14.8% 2.7 ‑3.1 0.031325yihRpredicted aldose‑1‑epimerase
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (11/15):  total (16/15)
Fisher's exact test strand distribution p-value = 4.34e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGGGTAACGCTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGG  >  REL606/4048190‑4048259
                                  |                                   
aGGGTAACGCTGGTGGCGGTTAATTCACTAATTTGc                                    <  1:2387611/36‑1 (MQ=255)
aGGGGAACGCTGGTGGCGGTTAATTTACTTAATTgt                                    >  1:101071/1‑35 (MQ=16)
   gTAACGCTGGTGGCGGTTAATTCACTAATTTGCCAg                                 <  1:309973/36‑1 (MQ=255)
      aCGCTGGTGGCGGTTAATTCACTAATTTTCCAGgct                              >  1:2774717/1‑36 (MQ=37)
       cGCTGGTGGCGGTTTATTTTCTGATTTTCCAGtctc                             >  1:4272694/1‑36 (MQ=16)
       cGCTGGTGGCGGCTAATTCCCTAATTTTTCCGtctc                             >  1:2351340/1‑36 (MQ=16)
        gCTGGTGGCGGTTAATTTACTCATTTTCCCGTCTCg                            >  1:1312890/1‑36 (MQ=19)
         cTGGTGGCTGTTAATTCTCTAATTTGCCAGTCTCGc                           <  1:5394980/36‑1 (MQ=37)
         cTGGTGGCGGTTAATTCACTAATTTTCCAGtctctc                           >  1:5146541/1‑34 (MQ=37)
          tggtggCGGTTAATTCACTAATTTGCCAGTCTCGcc                          <  1:1583506/36‑1 (MQ=255)
           ggtggCGGTTAATTCACTAATTTGCCAGTCTCGCCa                         <  1:2870150/36‑1 (MQ=255)
               gCGGTTAATTCACTAATTGGACAGTCTCGACAGGCa                     >  1:1929570/1‑36 (MQ=25)
                cGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCaa                    >  1:1008153/1‑36 (MQ=255)
                 ggTTAATTCACTAATTTGCCAGTCTCGCCAGGCAag                   <  1:816734/36‑1 (MQ=255)
                 ggTTAATTCACTAATTTGCCAGTCTCGCCAGGCAag                   >  1:4774025/1‑36 (MQ=255)
                 ggTTAATTCACTAATTTGCCAGTCTCGCCAGGCAag                   >  1:2055885/1‑36 (MQ=255)
                  gTTAATTCACTAATTTGCCAGTCTCGCCAGGCAagc                  <  1:5097339/36‑1 (MQ=255)
                   ttAATTCACTAATTTGCCAGTCTCTCCAGGCAagca                 >  1:1099019/1‑36 (MQ=255)
                   ttAATTCACTAATTTGCCAGTCTCGCCAGGCAagca                 >  1:148680/1‑36 (MQ=255)
                     aaTTCACTAATTTGCCAGTCTCGCCAGGCAagcagc               <  1:1034801/36‑1 (MQ=255)
                       ttCACTAATTTGCCAGTCTCGCCAGGCAAGCAGccc             >  1:2271780/1‑36 (MQ=255)
                       ttCACTAATTTGCCAGTCTCGCCAGGCAAGCAGccc             <  1:1445753/36‑1 (MQ=255)
                          aCTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTg          <  1:2807174/36‑1 (MQ=255)
                          aCTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTg          <  1:318656/36‑1 (MQ=255)
                           cTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGa         <  1:3005662/36‑1 (MQ=255)
                           cTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGa         <  1:256169/36‑1 (MQ=255)
                            tAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGaa        <  1:4738954/36‑1 (MQ=255)
                             aaTTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAAt       >  1:3810321/1‑36 (MQ=255)
                             aaTTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAAt       <  1:313246/36‑1 (MQ=255)
                              aTTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATa      >  1:4402738/1‑36 (MQ=255)
                              aTTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATa      >  1:4522593/1‑36 (MQ=255)
                                ttGCCAGTCTCGCCAGGCAAGCAGCCCGGGAATAGc    >  1:4011731/1‑36 (MQ=255)
                                  gCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCgg  >  1:1731344/1‑36 (MQ=255)
                                  |                                   
AGGGTAACGCTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGG  >  REL606/4048190‑4048259

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG