Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,863,739 C→A 30.6% D111Y (GAT→TAT)  bglB ← cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,863,7390C→A30.6% 2.2 ‑2.1 0.019331bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (5/8):  total (5/14)
Fisher's exact test strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TACCAGCGGCTTGATCCCCGCCTGCGCCATTTCATCAAACAGCCGATCGTAAAACGCTAACCCCGCTTC  >  REL606/3863704‑3863772
                                   |                                 
tACCAGCGGCTTGATCCCCGCCTGCGCCATTtcatc                                   <  1:3362568/36‑1 (MQ=255)
 aCCAGCGGCTTGATCCCCGCCTGCGCCATTtcatca                                  >  1:5150895/1‑36 (MQ=255)
   cAGCGGCTTGATCCCCGCCTGCGCCATTTCATCaaa                                <  1:582576/36‑1 (MQ=255)
   cAGCGGCTTGATCCCCGCCTGCGCCATTTCATCaaa                                <  1:4597501/36‑1 (MQ=255)
     gCGGCTTGATCCCCGCCTGCGCCATTTCATCAAACa                              <  1:73470/36‑1 (MQ=255)
       ggCTTGATCCCCGCCTGCGCCATTTCATCAAACAGc                            <  1:1509154/36‑1 (MQ=255)
       ggCTTGATCCCCGCCTGCGCCATCTCATCAAACAGc                            >  1:1881783/1‑36 (MQ=255)
        gCTTGATCCCCGCCTGCGCCATTTCATCAAACAGcc                           <  1:5554058/36‑1 (MQ=255)
        gCTTGATCCCCGCCTGCGCCATTTCATCAAACAGcc                           <  1:3859758/36‑1 (MQ=255)
          ttGATCCCCGCCTGCGCCATTTCATCAAACAGCCGa                         <  1:3472094/36‑1 (MQ=255)
          ttGATCCCCGCCTGCGCCATTTCATCAAACAGCCGa                         >  1:2945044/1‑36 (MQ=255)
            gtcccccGCCCGCGCAATTTCATCAAACAGCCGATc                       <  1:5819036/33‑1 (MQ=37)
            gATCCCCGCCTGCGCCATTTCATCAAACAGCCGATc                       >  1:1511002/1‑36 (MQ=255)
                  cGCCTGCGCCATTTCATAAAACAGCCGATCGTaaaa                 <  1:5930849/36‑1 (MQ=255)
                   gCCTGCGCCATTTCATCAAACAGCCGATCGTAAAAc                >  1:3870457/1‑36 (MQ=255)
                      tGCGCAATTACATAAAACAGCCGATCGTAAAACGCt             <  1:4707563/36‑1 (MQ=25)
                        caccATTACAAAAAACAGCCGATCGTAAAACGCTaa           <  1:1774418/34‑1 (MQ=21)
                          cgaTTTAATAAAACACGGGATCGTAAAACGCTAAcc         <  1:4525412/34‑1 (MQ=12)
                          ccATTACATAAAACAGCCGATCGTAAAACGCTAAcc         <  1:188897/36‑1 (MQ=37)
                             tttCATCAAACAGACGATCGTAAAACGCTAACCCCg      >  1:3073113/1‑36 (MQ=255)
                             tttAATAAAACAGCCGATCGTAAAACGCTAACCCCg      <  1:5896600/36‑1 (MQ=37)
                               tcatcaAACAGCCGATCGTAAAACGCTAACCCCGCt    >  1:4078856/1‑36 (MQ=255)
                                aagcaACCACCCGATCGTAAAACGCTACCCCCGCtt   <  1:5672824/33‑1 (MQ=19)
                                 atcaAACAGCCGATCGTAAAACGCTAACCCCGCTTc  >  1:1315929/1‑36 (MQ=255)
                                 accaACCAGCCGATCGTAAAACGCTAACCCCGCTTc  <  1:5565122/34‑1 (MQ=37)
                                   |                                 
TACCAGCGGCTTGATCCCCGCCTGCGCCATTTCATCAAACAGCCGATCGTAAAACGCTAACCCCGCTTC  >  REL606/3863704‑3863772

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG