Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,666,286 T→C 18.8% D284D (GAT→GAC xylR → DNA‑binding transcriptional activator, xylose‑binding

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,666,2860T→C18.8% 3.4 ‑4.4 ‑0.035852xylRDNA‑binding transcriptional activator, xylose‑binding
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (19/8):  total (21/14)
Fisher's exact test strand distribution p-value = 1.18e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.85e-01

GAACGGCGCTCAACAGATTATCGCTCGCTGACCGATCCCGCCGTTATTCAGGCCATGCATTACATTCGT  >  REL606/3666251‑3666319
                                   |                                 
gAACGGCGCTCAACAGATTATCGCTCGCTGGCCGAt                                   >  1:3758424/1‑36 (MQ=255)
gAACGGCGCTCAACAGATTATCGCTCGCTGAACGAt                                   >  1:1505694/1‑36 (MQ=255)
 aaCGGCGCTCAACAGATTATCGCTCGCTGACCGATc                                  >  1:348873/1‑36 (MQ=255)
  aCGGCGCTCAACAGATTATCGCTCGCTGACCGATcc                                 <  1:1386488/36‑1 (MQ=255)
      cgcTCAACAGATTATCGCTCGCTGACCGATCccgcc                             <  1:779934/36‑1 (MQ=255)
      cgcTCAACAGATTATCGCTCGCTGACCGATCccgcc                             >  1:4975724/1‑36 (MQ=255)
      cgcTCAACAGATTATCGCTCGCTGACCGATCccgcc                             <  1:2185780/36‑1 (MQ=255)
       gcTCAACAGATTATCGCTCGCTTATCGATCccgccg                            >  1:2976073/1‑36 (MQ=37)
       gcTCAACAGATTATCGCTCGCTGACCGATCccgccg                            <  1:2106368/36‑1 (MQ=255)
       gcTCAACAGATTATCGCTCGCTCAACGACTccgccg                            >  1:4274957/1‑36 (MQ=16)
        caaaaCAAATTACCGCTCGCTGCCCGATCCCGCCGt                           <  1:2130010/33‑1 (MQ=16)
         tCAACAGATTATCGCTCGCTGACCGATCCCGCCGtt                          >  1:2858217/1‑36 (MQ=255)
         tCAACAGATTATCGCTCGCTGACCGATCCCGCCGtt                          >  1:2867781/1‑36 (MQ=255)
         tCAACAGATTATCGCTCACTGACCGATCccgccgat                          >  1:5071427/1‑34 (MQ=38)
           aaCAGATTATCGCTCGCTGACCGATCCCGCCGTTAt                        >  1:2403074/1‑36 (MQ=255)
           aaCAGATTATCGCTCGCTGACCGATCCCGCCGTTAt                        >  1:3447255/1‑36 (MQ=255)
           aaCAAATTATCCCTCGTTGCCCCTTCCCGCCGTTAt                        <  1:928879/36‑1 (MQ=2)
            aaaGATTCCCGCTCGCTGACCACTCCCGCCGTTAtt                       <  1:2921684/34‑1 (MQ=16)
            aCAGATTATCGCTCGCGGACCGAGCTCTCCGCTAtg                       >  1:2982188/1‑35 (MQ=12)
             cAGATTATCGCTCGCTGACCGATCCCGCCGTTATTc                      >  1:1066100/1‑36 (MQ=255)
              aGATTATCGCTCGCTGACCGATCCCGCCGTTATTCa                     >  1:5501304/1‑36 (MQ=255)
               gATTATCGCTCGCTGACCGATCCCGCCGTTATTCAg                    >  1:2430352/1‑36 (MQ=255)
                  taccgctccctGACCGACCCCGCCGTTATTCAGGcc                 <  1:257247/33‑1 (MQ=25)
                    ccgctcgctAACCCACCCCGCCGTTATTCAGGCCAt               <  1:1332713/35‑1 (MQ=21)
                      gctcgctGACCGATCCCGCCGTTATTCAGGCCATGc             >  1:3296975/1‑36 (MQ=255)
                        ccgctGACAGACCCCGCCGTTATTCAGGCCATGCAt           <  1:5838486/35‑1 (MQ=37)
                          gctGCCCGACCCCGCCGTTATTCAGGCCATGCATTa         <  1:2732214/36‑1 (MQ=37)
                          gctGACCGATCCCGCCGTTATTCAGGCCATGCCTTa         >  1:416489/1‑36 (MQ=255)
                          gctGACCGATCCCGCCGTTATTCAGGCCATGCATTa         >  1:3280467/1‑36 (MQ=255)
                          gctGACCGATCCCGCCGTTATTCAGGCCATGCATTa         >  1:4951628/1‑36 (MQ=255)
                          gctGACCGATCCCGCCGTTATTCAGGCCATGCATTa         <  1:531196/36‑1 (MQ=255)
                          gctGACCGATCCCGCCCTTATTCAGGCCATGCATTa         >  1:949689/1‑36 (MQ=255)
                          gctGACCGATACAGCCGTTGTTTAGGCCATGCAGTa         >  1:2053510/1‑36 (MQ=16)
                            tGACCGATCCCGCCGTTATTCAGGCCATGCATTACa       >  1:5530380/1‑36 (MQ=255)
                            cccccGACCCCGCCGTTATTCAGGCCATGCATTACa       <  1:1903796/33‑1 (MQ=37)
                             gACCGACCCCGCCGTTATTCAGGCCATGCATTACAt      <  1:2903746/36‑1 (MQ=255)
                                 gATCCCGCCGTTATTCAGGCCATGCATTACATTCGt  >  1:5652170/1‑36 (MQ=255)
                                   |                                 
GAACGGCGCTCAACAGATTATCGCTCGCTGACCGATCCCGCCGTTATTCAGGCCATGCATTACATTCGT  >  REL606/3666251‑3666319

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