Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,608,001 T→C 10.6% I180T (ATC→ACC)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,608,0010T→C10.6% 6.4 ‑2.3 0.060539deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (28/25):  total (29/31)
Fisher's exact test strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAA  >  REL606/4607966‑4608035
                                   |                                  
cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                    <  1:3498770/36‑1 (MQ=255)
cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                    <  1:7331592/36‑1 (MQ=255)
cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                    >  1:7056824/1‑36 (MQ=255)
  cgACCCGTGATATTACCGCAACCGTGGACTCCATcc                                  <  1:6224453/36‑1 (MQ=255)
   gACCCGTGATATTACCGCAACCGTGGACTCCATccc                                 <  1:1414368/36‑1 (MQ=255)
   gACCCGTGATATTACCGCAACCGTGGACTCCATccc                                 >  1:6793794/1‑36 (MQ=255)
   gACCCGTGATATTACCGCAACCGTGGACTCCATccc                                 <  1:2838171/36‑1 (MQ=255)
    aCCCGTGATATTACCGCAACCGTGGACTCTATcccc                                >  1:912771/1‑35 (MQ=255)
    aCCCGTGATATTACCGCAACCGTGGACTCCATCCCg                                <  1:1525821/36‑1 (MQ=255)
     cccGTGATATTACCGCAACCGTGGACTCCATCCCGc                               <  1:436420/36‑1 (MQ=255)
     cccGTGATATTACCGCAACCGTGGACTCCATCCCGc                               <  1:3641046/36‑1 (MQ=255)
     cccGTGATATTACCGCAACCGTGGACTCCATCCCGc                               <  1:6558046/36‑1 (MQ=255)
      ccGTGATCTTACCGCAACCTTGGGCTCCATCCCGCt                              <  1:5819521/36‑1 (MQ=25)
      ccGTGATATTACCGCAACCGTGGACTCCATCCCGCt                              <  1:4941886/36‑1 (MQ=255)
      ccGTGATATTACCGCAACCGTGGACTCCATCCCGCt                              <  1:6620969/36‑1 (MQ=255)
       cGTGATATTACCGCAACCTTGGACTCCATCCCGCTg                             <  1:7435556/36‑1 (MQ=255)
       cGTGATATTACCGCAACCGTGGACTCCATCCCGCTg                             <  1:2762285/36‑1 (MQ=255)
       cGTGATATTACCGCAACCGTGGACTCCATCCCGCTg                             <  1:3041995/36‑1 (MQ=255)
       cGTGATATTACCGCAACCGTGGACTCCATCCCGCTg                             >  1:1893099/1‑36 (MQ=255)
        gTGATATTACCTCAACCGTGGACTCCATCCCGCTGa                            <  1:3929400/36‑1 (MQ=255)
        gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                            <  1:4825543/36‑1 (MQ=255)
        gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                            >  1:787554/1‑36 (MQ=255)
        gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                            <  1:6186566/36‑1 (MQ=255)
        gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                            <  1:322781/36‑1 (MQ=255)
        gTGATATTACCGCAACCGGGGACTCCACCCCGCTGa                            >  1:7024806/1‑36 (MQ=37)
         tGATATTACCGCAACCGTGGACTCCATCCCGCTGAt                           <  1:7368578/36‑1 (MQ=255)
         tGATATTACCGCAACCGTGGACTCCATCCCGCTGAt                           <  1:6225945/36‑1 (MQ=255)
          gATATTACCGCAACCGTGGACTCCATCCCGCTGATc                          >  1:3346909/1‑36 (MQ=255)
           atatTACCGCAACCGTGGACTCCATCCCGCTGATCa                         >  1:4845752/1‑36 (MQ=255)
            tatTACCGCAACCGTGGACTCCATCCCGCTGATCAc                        <  1:2002395/36‑1 (MQ=255)
             atTACCGCAACCGTGGACTCCATCCCGCTGATCAcc                       >  1:1991734/1‑36 (MQ=255)
             atTACCGCAACCGTGGACTCCATCCCGCTGATCAcc                       <  1:2797662/36‑1 (MQ=255)
              ttACCGCAACCGTGGACTCCATCCCGCTGATCACCg                      >  1:4748423/1‑36 (MQ=255)
               tACCGCAGCCGTGGACTCCATCCCGCTGATCACCGc                     <  1:6003973/36‑1 (MQ=255)
               tACCGCAACCGTGGACTCCATCCCTCTGATCACCGc                     >  1:2096110/1‑36 (MQ=255)
               tACCGCAACCGTGGACTCCATCCCGCTGATCACCGc                     <  1:5281556/36‑1 (MQ=255)
               tACCGCAACCGTGGACTCCATCCCGCTGATCACCGc                     <  1:5291246/36‑1 (MQ=255)
                aCCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                    >  1:2985322/1‑36 (MQ=255)
                aCCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                    >  1:7449198/1‑36 (MQ=255)
                 ccGCAACCGTGGACTCCATCCCGCTGATCACCGCct                   <  1:4804978/36‑1 (MQ=255)
                   gCAACCGTGGACTCCATCCCGCTGATCACCGCctct                 >  1:202047/1‑36 (MQ=255)
                   gCAACCGTGGACTCCATCCCGCTGATCACCGCctct                 >  1:6399451/1‑36 (MQ=255)
                    cAACCGTGCACCCCACCCCCCTGATCACCGCCTCTa                <  1:3243988/36‑1 (MQ=19)
                     aaCCGTGGACTCCATCCCGCTGATCACCGCCTCTAt               >  1:2500680/1‑36 (MQ=255)
                      cccGTGGCCCCCCCCCCCCTAACCACCGCCTCTAtt              <  1:4531259/35‑1 (MQ=2)
                      aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt              >  1:77326/1‑36 (MQ=255)
                      aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt              >  1:1116704/1‑36 (MQ=255)
                      aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt              >  1:6225799/1‑36 (MQ=255)
                      aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt              >  1:7366047/1‑36 (MQ=255)
                       ccGTGGCCCCCACCCCGCTAACCACCGCCTCTATTc             <  1:7326527/36‑1 (MQ=16)
                       ccGTGCACCCCACCCCGCTGACCACCGCCTCTATTc             <  1:1060493/36‑1 (MQ=16)
                         gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg           >  1:2582853/1‑36 (MQ=255)
                         gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg           >  1:2940173/1‑36 (MQ=255)
                         gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg           >  1:300845/1‑36 (MQ=255)
                         gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg           >  1:5997132/1‑36 (MQ=255)
                            accTCCACCCCGCTGATCACCGCCTCTATTCTGGCg        <  1:2119171/34‑1 (MQ=37)
                             aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga       >  1:6608897/1‑36 (MQ=255)
                             aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga       >  1:3132422/1‑36 (MQ=255)
                             aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga       >  1:3585415/1‑36 (MQ=255)
                             aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga       >  1:5017101/1‑36 (MQ=255)
                              ccccACCCCGCTGATCACCGCCTCTATTCTGGCgaa      <  1:3464810/34‑1 (MQ=37)
                               tCCATCCCGCTGATCACCGCCTCTATTCTGGCgaag     >  1:7117824/1‑36 (MQ=255)
                                  aTCCCGCTGATCACCGCCTCTATTCTGGCGAAGaaa  >  1:4609876/1‑36 (MQ=255)
                                   |                                  
CGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAA  >  REL606/4607966‑4608035

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG