Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,597,450 C→A 20.9% E411* (GAA→TAA)  ydfI ← predicted mannonate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,597,4500C→A20.9% 5.5 ‑3.9 0.0561231ydfIpredicted mannonate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (1/13):  ref base (18/24):  total (19/37)
Fisher's exact test strand distribution p-value = 2.08e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGTGGGTCACTGATTTCTATCGGATTTTCCTGTTCATCAACACCACCGACATAACGCATCCAACCCGCGA  >  REL606/1597415‑1597485
                                   |                                   
cAGTGGGTCACTGATTTCTATCGGATTTTCCTGTtc                                     >  1:2556594/1‑36 (MQ=255)
cAGTGGGTCACTGATTTCTATCGGATTTTCCTGTtc                                     <  1:5825694/36‑1 (MQ=255)
cAGTGGGTCACTGATTTCTATCGGATTTTCCTGTtc                                     <  1:88923/36‑1 (MQ=255)
  gTGGGTCACTGATTTCTATCGGATTTTCCTGTtcat                                   <  1:4764707/36‑1 (MQ=255)
    gggTCACTGATTTCTATCGGATTTTCCTGTtcatca                                 <  1:1297308/36‑1 (MQ=255)
    gggTCACTGATTTCTATCGGATTTTCCTGTtcatca                                 <  1:4884722/36‑1 (MQ=255)
      gTCACTGATTTCTATCGGATTTTCCTGTTCATCAac                               <  1:151073/36‑1 (MQ=255)
      gTCACTGATTTCTATCGGATTTTCCTGTTCATCAac                               >  1:6306383/1‑36 (MQ=255)
      gTCACTGATTTCTATCGGATTTTCCTGTTCATCAac                               <  1:5115321/36‑1 (MQ=255)
      gTCACTGATTTCTATCGGATTTTCCTGTTCATCAac                               <  1:3062619/36‑1 (MQ=255)
       tCACTGATTTCTATCGGATTTTCCTGTTCATCAaca                              <  1:6175082/36‑1 (MQ=255)
       tCACTGATTTCTATCGGATTTTCCTGTTCATCAaca                              >  1:249778/1‑36 (MQ=255)
        cccTGATTTCTATCGGATTTTCCTGTTCATCAacac                             <  1:5589107/34‑1 (MQ=255)
          cTGATTTCTATCGGATTTTCCTGTTCATCAAcacca                           >  1:5946937/1‑36 (MQ=255)
           tGATTTCTATCGGATTTTCCTGTTCATCAAcaccac                          >  1:1297337/1‑36 (MQ=255)
           tGATTTCTATCGGATTTTCCTGTTCATCAAcaccac                          <  1:703693/36‑1 (MQ=255)
            gATTTCTATCGGATTTTCCTGTTCATCAACaccacc                         >  1:4944883/1‑36 (MQ=255)
            gATTTCTATCGGATTTTCCTGTTCATCAACaccacc                         >  1:4137355/1‑36 (MQ=255)
              tttCTATCGGATTTTCCTGTTCATCAACACCACCGa                       <  1:2044293/36‑1 (MQ=255)
              tttCTATCGGATTTTCCTGTTCATCAACACCACCGa                       <  1:1996375/36‑1 (MQ=255)
               ttCTATCGGATTTTCCTGTTCATCAACACCACCGAc                      <  1:3625833/36‑1 (MQ=255)
                tCTATCGGATTTTCCTGTTCATCAACACCACCGACa                     >  1:7344354/1‑36 (MQ=255)
                tCTATCGGATTTTCCTGTTCATCAACACCACCGACa                     >  1:3165799/1‑36 (MQ=255)
                  tATCGGATTTTCCTGTTCATAAACACCACCGACATa                   <  1:5380710/36‑1 (MQ=255)
                  tATCGGATTTCCCTGTTAATCAACACCACCGACATa                   <  1:3915613/36‑1 (MQ=37)
                   aTCGGATTTTCCTGTTCATCAACACCACCGACATaa                  >  1:3577060/1‑36 (MQ=255)
                   aTCGGATTTTCCTGTTCATAAACACCACCGACATaa                  <  1:518107/36‑1 (MQ=255)
                   aTCGGATTTTCCTGTTAATCAACACCACCGACATaa                  >  1:12968/1‑36 (MQ=255)
                   aTCGAATTTTCCTGTTAATCAACACCACCGACATaa                  <  1:3349980/36‑1 (MQ=37)
                    tCGGATTTTCCTGTTAATAAACACCACCGACATAAc                 <  1:6619789/36‑1 (MQ=37)
                    tCGAATTTTCCTGTTAAAAAACACCACCGACATAAc                 <  1:4631379/36‑1 (MQ=19)
                     cGGATTTTCCTGTTCATCAACACCACCGACATAACg                <  1:2252191/36‑1 (MQ=255)
                     cGGATTTTCCTGTTAATAAACACCACCGACATAACg                <  1:4656679/36‑1 (MQ=37)
                      ggcttttCCTGTTAATCAACACCACCTCCATAACGc               <  1:6554524/33‑1 (MQ=16)
                       gATTTTCCTGTTCATCAACACCACCGACATAACGCa              <  1:1094113/36‑1 (MQ=255)
                        aTTTTCCTGTTCATAAACACCACCGACATAACGCAt             <  1:7223482/36‑1 (MQ=255)
                        aTTTTCCTGTTAATAAACACCACCGACATAACGCAt             <  1:4763629/36‑1 (MQ=38)
                         ttttCCTGTTCATCAACACCACCGACATAACGCATc            >  1:1507253/1‑36 (MQ=255)
                         ttttCCTGTTCATAAACACCACCGACATAACGCATc            <  1:7085297/36‑1 (MQ=255)
                         ttttCCTATTCATCAACACCACCGACATAACGCATc            <  1:5104738/36‑1 (MQ=255)
                          tttCCTGTTCATCAACACCACCGACATAACGCATcc           >  1:6649080/1‑36 (MQ=255)
                          tttCCTGTTCATCAACACCACCGACATAACGCATcc           >  1:5861350/1‑36 (MQ=255)
                          tttCCTGTTCAAAAACACCACCGACATAACGCATcc           <  1:1994703/36‑1 (MQ=37)
                          ttcccTGTTCATAAACACCCCCGACATAACGCATcc           <  1:2581768/33‑1 (MQ=25)
                           ttCCTGTTCATAAACACCACCGACATAACGCATCCa          <  1:810660/36‑1 (MQ=255)
                           ttCCTGTTCATAAACACCACCGACATAACGCATCCa          <  1:1798974/36‑1 (MQ=255)
                           ttCCTGTTCATAAACACCACCGACATAACGCATCCa          <  1:3519011/36‑1 (MQ=255)
                           ttCCTGTTAATCAACACCACCGACATAACGCATCCa          <  1:7337575/36‑1 (MQ=255)
                           ttCCTGTTAAAAAACACCACCGACATAACGCATCCa          <  1:1462462/36‑1 (MQ=25)
                            tCCTGTTAATCAACACCACCGACATAACGCATCCaa         <  1:3684344/36‑1 (MQ=255)
                            tCCTGTTAATAAAAACAACCGACATAACGCATCCaa         <  1:7254885/36‑1 (MQ=19)
                            tCCTATTAAACAACACCACCGACATAACGCATCCaa         <  1:3764287/36‑1 (MQ=37)
                             ccTGTACAAAAACACCACCGACATAACGCATCCAAc        <  1:4862233/36‑1 (MQ=25)
                               ttttCAACAAAACCACCGACATAACGCATCCAAccc      <  1:427820/34‑1 (MQ=37)
                               tGTTCATCAACACCACCGACATAACGCATCCAAccc      >  1:2972941/1‑36 (MQ=255)
                               tGTTCATCAACACCACCGACATAACGCATCCAAccc      >  1:4169729/1‑36 (MQ=255)
                               tGTTCATCAACAACACCGACATAACGCATCCAAccc      >  1:1954272/1‑36 (MQ=255)
                               tGTTCACCAACACCACCGACATAACGCATCCAAccc      <  1:642463/36‑1 (MQ=255)
                               tGTTAAACAACACCACCGACATAACGCATCCAAccc      <  1:6382181/36‑1 (MQ=37)
                                gTTCATCAACACCACCGACATAACGCATCCAACCcg     >  1:1005040/1‑36 (MQ=255)
                                 ttCATCAACACCACCGACATAACGCATCCAACCcgc    >  1:7428291/1‑36 (MQ=255)
                                 ttCATCAACACCACCGACATAACGCATCCAACCcgc    >  1:4087233/1‑36 (MQ=255)
                                  tcatcaACACCACCGACATAACGCATCCAACCcgcg   >  1:755200/1‑36 (MQ=255)
                                   catcaACACCACCGACATAACGCCTCCAACCCGCGa  >  1:5318103/1‑36 (MQ=255)
                                   catcaACACCACCGACATAACGCATCCAACCCGCGa  >  1:3940621/1‑36 (MQ=255)
                                   |                                   
CAGTGGGTCACTGATTTCTATCGGATTTTCCTGTTCATCAACACCACCGACATAACGCATCCAACCCGCGA  >  REL606/1597415‑1597485

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG