Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,607,9970T→C 6.0% 0.8 ‑2.2 ‑0.064535deoAthymidine phosphorylase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (31/29):  total (31/33)
Fisher's exact test strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAG  >  REL606/4607963‑4608032
                                  |                                   
cTACGCGACACGTGATATTACTGCAACCGGGGACTc                                    >  1:2323492/1‑36 (MQ=25)
  aCGCGACCCGTGATATTACCGCAACCGTGGACTCCa                                  >  1:3054995/1‑36 (MQ=255)
  aCGCGACCCGTGATATTACCGCAACCGTGGACTCCa                                  <  1:3020384/36‑1 (MQ=255)
  aCGCGACCCGTGATATTACCGCAACCGTGGACTCCa                                  >  1:2495494/1‑36 (MQ=255)
   cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                 <  1:3498770/36‑1 (MQ=255)
   cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                 >  1:7056824/1‑36 (MQ=255)
   cgcgACCCGTGATATTACCGCAACCGTGGACTCCAt                                 <  1:7331592/36‑1 (MQ=255)
     cgACCCGTGATATTACCGCAACCGTGGACTCCATcc                               <  1:6224453/36‑1 (MQ=255)
      gACCCGTGATATTACCGCAACCGTGGACTCCATccc                              <  1:1414368/36‑1 (MQ=255)
      gACCCGTGATATTACCGCAACCGTGGACTCCATccc                              <  1:2838171/36‑1 (MQ=255)
      gACCCGTGATATTACCGCAACCGTGGACTCCATccc                              >  1:6793794/1‑36 (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:6558046/36‑1 (MQ=255)
        cccGTGATATTACCGCAACCGTGGACTCCATCCCGc                            <  1:3641046/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:6186566/36‑1 (MQ=255)
           gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                         >  1:787554/1‑36 (MQ=255)
           gTGATATTACCGCAACCGTGGACTCCATCCCGCTGa                         <  1:322781/36‑1 (MQ=255)
           gTGATATTACCGCAACCGGGGACTCCACCCCGCTGa                         >  1:7024806/1‑36 (MQ=37)
            tGATATTACCGCAACCGTGGACTCCATCCCGCTGAt                        <  1:6225945/36‑1 (MQ=255)
            tGATATTACCGCAACCGTGGACTCCATCCCGCTGAt                        <  1:7368578/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:2797662/36‑1 (MQ=255)
                atTACCGCAACCGTGGACTCCATCCCGCTGATCAcc                    >  1:1991734/1‑36 (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:5291246/36‑1 (MQ=255)
                  tACCGCAACCGTGGACTCCATCCCGCTGATCACCGc                  <  1:5281556/36‑1 (MQ=255)
                   aCCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                 >  1:7449198/1‑36 (MQ=255)
                   aCCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                 >  1:2985322/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:6225799/1‑36 (MQ=255)
                         aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt           >  1:1116704/1‑36 (MQ=255)
                         aCCGTGGACTCCATCCCGCTGATCACCGCCTCTAtt           >  1:77326/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:300845/1‑36 (MQ=255)
                            gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg        >  1:2940173/1‑36 (MQ=255)
                            gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg        >  1:2582853/1‑36 (MQ=255)
                            gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg        >  1:5997132/1‑36 (MQ=255)
                               accTCCACCCCGCTGATCACCGCCTCTATTCTGGCg     <  1:2119171/34‑1 (MQ=37)
                                aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga    >  1:3132422/1‑36 (MQ=255)
                                aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga    >  1:3585415/1‑36 (MQ=255)
                                aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga    >  1:5017101/1‑36 (MQ=255)
                                aCTCCATCCCGCTGATCACCGCCTCTATTCTGGCga    >  1:6608897/1‑36 (MQ=255)
                                  tCCATCCCGCTGATCACCGCCTCTATTCTGGCgaag  >  1:7117824/1‑36 (MQ=255)
                                  |                                   
CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAG  >  REL606/4607963‑4608032

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