Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,255,385 A→T 11.1% V488E (GTA→GAA)  yhbX ← predicted hydrolase, inner membrane

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,255,3850A→T11.1% 3.2 ‑3.8 ‑0.0481463yhbXpredicted hydrolase, inner membrane
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (18/24):  total (24/24)
Fisher's exact test strand distribution p-value = 2.19e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.76e-01

TTCCTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATCCACGCATTAATAAGGTAATTATTGTA  >  REL606/3255351‑3255419
                                  |                                  
ttCCTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTAc                                   <  1:2731059/36‑1 (MQ=255)
 tCCTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTAcc                                  <  1:902393/36‑1 (MQ=255)
   cTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTAcccc                                <  1:3934016/36‑1 (MQ=255)
   cTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTAcccc                                <  1:1917109/36‑1 (MQ=255)
   cTCAAGCGTTTGCGGCTGGTGCGGCTTTTTTTcccc                                >  1:5626848/1‑36 (MQ=19)
    tCAAGCGTTTGCGGCTGTTTCGGCTTTTTTTcaccc                               >  1:1949479/1‑36 (MQ=16)
    tCAAGCGTTTGCGGCTGTTCCGGCTTTGTTAccccc                               <  1:5417347/36‑1 (MQ=255)
    ctaaGCGTTTGCGGCTGTTCCGGCTTTGTTAccccc                               <  1:7811447/34‑1 (MQ=255)
      aaGCGTTTGCGGCTGTTTCGGCTTTTTTTtccccct                             >  1:3752580/1‑34 (MQ=12)
      aaGCGTTTGCGGCTGTTCCGGCTTTTTTAcccccct                             >  1:1453135/1‑34 (MQ=37)
      aaGCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCAt                             <  1:4327481/36‑1 (MQ=255)
       aGCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATc                            <  1:4765381/36‑1 (MQ=255)
       aGCGTTTGCGGCTGGTTTGGGTTTTTTTTCCCCATc                            >  1:4914743/1‑36 (MQ=2)
        gCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATcc                           <  1:6249813/36‑1 (MQ=255)
        gCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATcc                           <  1:1454579/36‑1 (MQ=255)
        gCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATcc                           <  1:3355611/36‑1 (MQ=255)
        gCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATcc                           <  1:3592770/36‑1 (MQ=255)
         cGTTTGCGGCTGTTCCGGGTTTTTTATCCCCCTccc                          >  1:5118883/1‑35 (MQ=16)
         cGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATCCa                          <  1:5719222/36‑1 (MQ=255)
          gTTTGCGGCTGTTTCGGGTTTTTTACCCCCATCCAc                         >  1:374357/1‑36 (MQ=25)
          gTTTGCGGCTGTTTCGGCTTTGTTTCCCCCCCCCAc                         >  1:3368911/1‑36 (MQ=16)
          gTTTGCGGCTGTTCCGGCTTTGTTACCCCCATCCAc                         <  1:3673279/36‑1 (MQ=255)
           tttGCGTCTGTTCCGGCTTTGTTACCCCCATCCACg                        <  1:1727122/36‑1 (MQ=255)
           tttGCGGCTGTTCCGGCTTTTTTTCCCCCCTCCACg                        >  1:6205295/1‑36 (MQ=25)
            ttGCGGCTGTTCCGGCTTTGTTACCCCCATCCACGc                       <  1:2490409/36‑1 (MQ=255)
            ttGCGGCTGTTCCGGCTTTGTTACCCCCATCCACGc                       <  1:7365272/36‑1 (MQ=255)
                ggCTGTTCCGGCTTTGTTACCCCCATCCACGCATta                   >  1:4862713/1‑36 (MQ=255)
                 gCTGTTCCGGCTTTGTTACCCCCATCCACGCATtaa                  <  1:7681434/36‑1 (MQ=255)
                 gCTGTTCCGGCTTTGTTACCCCCATCCACGCATtaa                  <  1:113410/36‑1 (MQ=255)
                 gCTGTTCCGGCTTTGTTACCCCCATCCACGCATtaa                  <  1:2106883/36‑1 (MQ=255)
                  cTGTTCCGGCTTTGTTACCCCCATCCACGCATtaat                 >  1:2993502/1‑36 (MQ=255)
                   tGTTCCGGCTTTGTTACCCCCATCCACGCATtaata                >  1:3461813/1‑36 (MQ=255)
                    gTTCCGGCTTTGTTACCCCCATCCACGCATtaataa               >  1:2245213/1‑36 (MQ=255)
                      tCCGGCTTTGTTACCCCCATCCACGCATTAATAAgg             >  1:6749893/1‑36 (MQ=255)
                        cGGCTTTGTTACCCCCATCCACGCATTAATAAGGTa           >  1:6591454/1‑36 (MQ=255)
                         ggCTTTGTTACCCCCATCCACGCATTAATAAGGTaa          >  1:6731616/1‑36 (MQ=255)
                         ggCTTTGTTACCCCCATCCACGCATTAATAAGGTaa          >  1:1684129/1‑36 (MQ=255)
                           cTTTGTTACCCCCATCCACGCATTAATAAGGTAatt        <  1:6308909/36‑1 (MQ=255)
                           cTTTGTTACCCCCATCCACGCATTAATAAGGTAatt        <  1:77640/36‑1 (MQ=255)
                            tttGTTACCCCCATCCACGCATTAATAAGGTAatta       >  1:4454136/1‑36 (MQ=255)
                              tGTTACCCCCATCCACGCATTAATAAGGTAattatt     >  1:263293/1‑36 (MQ=255)
                              tGTTACCCCCATCCACGCATTAATAAGGTAattatt     >  1:5523260/1‑36 (MQ=255)
                              tGTTACCCCCATCCACGCATTAATAAGGTAattatt     <  1:4073763/36‑1 (MQ=255)
                               gTTACCCCCATCCACGCATTAATAAGGTAATTATTg    <  1:3976425/36‑1 (MQ=255)
                               gTTACCCCCATCCACGCATTAATAAGGGAATTATTg    >  1:7274242/1‑36 (MQ=255)
                                ttACCCCCATCCACGCATTAATAAGGTAATTATTGt   >  1:3064290/1‑36 (MQ=255)
                                ttACCCCCATCCACGCATTAATAAGGTAATTATTGt   >  1:1864461/1‑36 (MQ=255)
                                 tACCCCCATCCACGCATTAATAAGGTAATTATTGTa  >  1:7678705/1‑36 (MQ=255)
                                  |                                  
TTCCTCAAGCGTTTGCGGCTGTTCCGGCTTTGTTACCCCCATCCACGCATTAATAAGGTAATTATTGTA  >  REL606/3255351‑3255419

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG