Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,500,852 A→G I641T (ATT→ACT)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,500,8520AG99.0% 1231.1 / ‑0.7 292I641T (ATT→ACT) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (170/119);  total (171/121)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 302 total aligned reads displayed.

CCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTAC  >  REL606/3500779‑3500925
                                                                          |                                                                         
ccATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTc                                                                         <  2:4829203/76‑1 (MQ=255)
 cATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCa                                                                        >  2:7105691/1‑76 (MQ=255)
  aTGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCAc                                                                       >  2:2706986/1‑76 (MQ=255)
  aTGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCAc                                                                       >  2:6506724/1‑76 (MQ=255)
  aTGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCAc                                                                       >  1:3054265/1‑76 (MQ=255)
    gACTACCTCCTTTTTTATCATCTCTGAAACCCCCACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTg                                                                     <  2:6649828/76‑1 (MQ=255)
     aCTAACTCCTTTTTTAT‑‑TCTCTGGAACACACACAATCCGTGTGCAGGCCCCTGCCAGACAGCCGTAGTCACTGGGt                                                                  <  2:1337756/76‑1 (MQ=255)
      cTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTggg                                                                   >  2:1841354/1‑76 (MQ=255)
       tAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGt                                                                  >  1:3559429/1‑76 (MQ=255)
       tAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGt                                                                  >  1:8124055/1‑76 (MQ=255)
        aaCTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGtt                                                                 >  2:6052035/1‑76 (MQ=255)
          cTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTAt                                                               >  1:7026622/1‑76 (MQ=255)
           tCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATc                                                              <  1:4829552/76‑1 (MQ=255)
           tCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATc                                                              >  1:6229513/1‑76 (MQ=255)
           tCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATc                                                              >  1:5998352/1‑76 (MQ=255)
              ttttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCctc                                                           <  1:1944548/76‑1 (MQ=255)
               tttttATCATCTCTGGAACACACACAACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCctct                                                          <  2:4423328/76‑1 (MQ=255)
                ttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCtctc                                                         >  2:4646426/1‑76 (MQ=255)
                ttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCtctc                                                         >  1:5688757/1‑76 (MQ=255)
                ttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCtctc                                                         >  2:6120725/1‑76 (MQ=255)
                 tttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTcc                                                        >  2:6827367/1‑76 (MQ=255)
                 tttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTcc                                                        <  2:7991599/76‑1 (MQ=255)
                 tttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTcc                                                        >  2:95024/1‑76 (MQ=255)
                 tttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTGATGCTCTcc                                                        >  2:4957190/1‑76 (MQ=255)
                  ttATCCTCTCTGGAACACACACAATCCGGGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCa                                                       <  1:243105/76‑1 (MQ=255)
                  ttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCa                                                       >  2:3490817/1‑76 (MQ=255)
                  ttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCa                                                       >  1:369322/1‑76 (MQ=255)
                  ttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCa                                                       >  1:4172890/1‑76 (MQ=255)
                  ttATCATCTCTGGAACACACACAACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCa                                                       <  1:1188221/76‑1 (MQ=255)
                   tATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAg                                                      >  2:7269132/1‑76 (MQ=255)
                    atcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGc                                                     >  1:803878/1‑76 (MQ=255)
                    atcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGc                                                     >  1:3471821/1‑76 (MQ=255)
                    atcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGc                                                     >  2:7443722/1‑76 (MQ=255)
                    atcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGc                                                     >  1:2624600/1‑76 (MQ=255)
                    atcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGc                                                     >  2:2950762/1‑76 (MQ=255)
                      catcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGa                                                   >  1:7605403/1‑76 (MQ=255)
                      catcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGa                                                   >  1:5057308/1‑76 (MQ=255)
                      catcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGa                                                   >  1:4673348/1‑76 (MQ=255)
                       atcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGaa                                                  >  2:2834045/1‑76 (MQ=255)
                        tctcTGAAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAAt                                                 <  1:2162768/76‑1 (MQ=255)
                         ctctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATg                                                >  2:2235104/1‑76 (MQ=255)
                          tctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATgg                                               >  1:58158/1‑76 (MQ=255)
                          tctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATgg                                               >  1:1363658/1‑76 (MQ=255)
                          tctGAAACACACACAACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATgg                                               <  1:6397583/76‑1 (MQ=255)
                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATggg                                              >  2:3253418/1‑76 (MQ=255)
                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATggg                                              >  1:1002088/1‑76 (MQ=255)
                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATggg                                              >  2:751618/1‑76 (MQ=255)
                            tGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATgggg                                             >  1:1557337/1‑76 (MQ=255)
                             ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATggggg                                            >  2:2161638/1‑76 (MQ=255)
                             ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATggggg                                            >  2:1623333/1‑76 (MQ=255)
                              gaccacacacaACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGa                                           <  1:6577410/73‑1 (MQ=255)
                              gAACACACACAACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGa                                           <  1:8093607/76‑1 (MQ=255)
                               aaCACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGaa                                          >  2:596979/1‑76 (MQ=255)
                                acccacacaACCCGTGTGCAGTTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAAt                                         <  1:5223793/73‑1 (MQ=255)
                                acacacacAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAAt                                         >  1:4910001/1‑76 (MQ=255)
                                acacacacAACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAAt                                         <  1:2310112/76‑1 (MQ=255)
                                 cacacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATg                                        <  2:3941388/76‑1 (MQ=255)
                                 cacacacaATCCGTGTACAGGTCCCTGCCAGCCAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATg                                        <  2:7675606/76‑1 (MQ=255)
                                  acacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATgg                                       >  1:7222510/1‑76 (MQ=255)
                                    acacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCg                                     >  2:5777564/1‑76 (MQ=255)
                                      acaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTg                                   >  1:2718340/1‑76 (MQ=255)
                                      acaACCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTg                                   <  1:3797044/76‑1 (MQ=255)
                                        aaTCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGcc                                 >  2:5179434/1‑76 (MQ=255)
                                          tCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCAc                               <  1:5326080/76‑1 (MQ=255)
                                          tCCGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCAc                               <  1:6806535/76‑1 (MQ=255)
                                            cGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTc                             >  2:8170326/1‑76 (MQ=255)
                                            cGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTc                             >  1:5818842/1‑76 (MQ=255)
                                            cGTGTGCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTc                             >  2:4043067/1‑76 (MQ=255)
                                               gtgCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTcccc                          <  2:4132238/76‑1 (MQ=255)
                                                tgCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCCATCCCca                         >  2:5510120/1‑76 (MQ=255)
                                                tgCAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCca                         >  1:7365421/1‑76 (MQ=255)
                                                 gCAGGCCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcag                        <  2:5867456/76‑1 (MQ=255)
                                                  cAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagc                       >  2:3058582/1‑76 (MQ=255)
                                                  cAGGTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagc                       <  2:6286997/76‑1 (MQ=255)
                                                     gccccTGCCAGCCAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCTCCAGCAgg                    <  1:3323433/74‑1 (MQ=255)
                                                     gTCCCTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAgg                    >  2:77808/1‑76 (MQ=255)
                                                      ccccTCCCGGCCAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGt                   <  2:203958/75‑1 (MQ=255)
                                                      ccccGGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGt                   <  1:3958541/75‑1 (MQ=255)
                                                       cccTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGATTGGCGTGCCACTCCCCAGCAggga                  >  2:3597641/1‑74 (MQ=255)
                                                       cccTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTa                  >  2:1572046/1‑76 (MQ=255)
                                                        ccTGCCAAGCCAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTa                  <  2:7482437/76‑1 (MQ=255)
                                                        ccTGCCAGCCAGCCGTAGTCACTGGGTTATCCTCACCAGCGAATGGGGGAATGGCGTGCCACCCCCCAGCAGGTaa                 <  2:6948933/76‑1 (MQ=255)
                                                        ccTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaa                 <  1:3265876/76‑1 (MQ=255)
                                                        ccTGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaa                 <  2:1883384/76‑1 (MQ=255)
                                                        ccTCCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaa                 <  1:5771972/76‑1 (MQ=255)
                                                         cggccagccaGCCGTAGTCACTGGTTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaa                <  1:431540/74‑1 (MQ=255)
                                                         cTGCCAGCCAGCCTTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaa                <  1:7711773/76‑1 (MQ=255)
                                                          tGCCAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaaa               >  2:8254073/1‑76 (MQ=255)
                                                          tGCAAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaaa               <  2:2536169/76‑1 (MQ=255)
                                                           gccagccaGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAAc              <  1:4594728/76‑1 (MQ=255)
                                                            ccAGACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTACCACTCCCCAGCAGGTAAAACa             >  2:1053563/1‑76 (MQ=255)
                                                                gccAGCCTTATTCCCTGGTTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtc          <  2:1235772/74‑1 (MQ=255)
                                                                gACAGCCGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtc          <  1:3897517/76‑1 (MQ=255)
                                                                      cGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt    <  1:6054190/76‑1 (MQ=255)
                                                                      cGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt    >  2:2354876/1‑76 (MQ=255)
                                                                      cGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt    >  1:1690912/1‑76 (MQ=255)
                                                                      cGTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt    <  1:864884/76‑1 (MQ=255)
                                                                      cGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt    >  2:3420794/1‑76 (MQ=255)
                                                                       gTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTa   <  1:1820296/76‑1 (MQ=255)
                                                                       gTAGTCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTa   <  2:640008/76‑1 (MQ=255)
                                                                          |                                                                         
CCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTAC  >  REL606/3500779‑3500925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: