Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,070,984 T→C 5.9% I630T (ATT→ACT)  glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,070,9840TC 5.9% 289.7 / 7.7 102I630T (ATT→ACT) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base T (47/49);  new base C (4/2);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 111 total aligned reads displayed.

TTTTGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGATTGTG  >  AP009048/4070915‑4071053
                                                                     |                                                                     
ttttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAtt                                                                      >  1:7278759/1‑71 (MQ=255)
ttttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAtt                                                                      >  1:6529780/1‑71 (MQ=255)
ttttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAtt                                                                      >  1:6676174/1‑71 (MQ=255)
ttttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAtt                                                                      >  1:2845256/1‑71 (MQ=255)
 tttGATCGCAATTAACGCCATGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTc                                                                     >  1:2582300/1‑71 (MQ=255)
 tttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAt                                                                       <  1:7930430/69‑1 (MQ=255)
 tttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTc                                                                     >  1:2274432/1‑71 (MQ=255)
 tttGATCGCAATTAACG‑‑ACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTccc                                                                   >  1:6128497/1‑71 (MQ=255)
  ttGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTcc                                                                    >  1:7018638/1‑71 (MQ=255)
     aTCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTcccc                                                                  >  1:3896118/1‑70 (MQ=255)
       cGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCAt                                                               >  1:5460491/1‑71 (MQ=255)
        gCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCAtt                                                              <  1:6210915/71‑1 (MQ=255)
        gCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCAtt                                                              >  1:3363929/1‑71 (MQ=255)
        gCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCAtt                                                              <  1:2154538/71‑1 (MQ=255)
        gCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCAtt                                                              >  1:2155683/1‑71 (MQ=255)
          aaTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCg                                                            >  1:934418/1‑71 (MQ=255)
          aaTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCg                                                            >  1:1521085/1‑71 (MQ=255)
          aaTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCg                                                            >  1:2782101/1‑71 (MQ=255)
          aaTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCg                                                            >  1:6269138/1‑71 (MQ=255)
          aaTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCg                                                            >  1:2317453/1‑71 (MQ=255)
           aTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTc                                                                     <  1:5345761/61‑1 (MQ=255)
           aTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGc                                                           >  1:1384798/1‑71 (MQ=255)
           aTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGc                                                           >  1:6793701/1‑71 (MQ=255)
           aTTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGc                                                           >  1:1586168/1‑71 (MQ=255)
            ttAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCt                                                          >  1:2841118/1‑71 (MQ=255)
            ttAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCt                                                          <  1:4528444/71‑1 (MQ=255)
            ttAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCt                                                          <  1:7770552/71‑1 (MQ=255)
               aCGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTg                                                         >  1:5705385/1‑69 (MQ=255)
                 gCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCTCTGgaga                                                     <  1:349367/71‑1 (MQ=255)
                 gCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCACCCATTCGCTGgaga                                                     <  1:6280408/71‑1 (MQ=255)
                    aCGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCAt                                                                       <  1:6265523/50‑1 (MQ=255)
                     cGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAgg                                                   >  1:3871670/1‑69 (MQ=255)
                       cTTGATGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCAGGCCATTCCCCCATTCg                                                            <  1:5313353/58‑1 (MQ=39)
                             gTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGt                                         <  1:4979237/71‑1 (MQ=255)
                             gTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGt                                         <  1:2396357/71‑1 (MQ=255)
                                aCAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCa                                           >  1:6151480/1‑66 (MQ=255)
                                aCAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCa                                           >  1:2783351/1‑66 (MQ=255)
                                aCAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGAt                                      <  1:3237279/71‑1 (MQ=255)
                                aCAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGAt                                      <  1:6700468/71‑1 (MQ=255)
                                 cAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTc                                                             <  1:1224310/47‑1 (MQ=255)
                                  agagATTGTTTTACCTGCTGGGGAGTGGCAcgccattccgccattc                                                             <  1:1231140/46‑1 (MQ=255)
                                   gagaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGga                                                       >  1:378656/1‑51 (MQ=255)
                                   gagaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTAc                                   <  1:3620276/71‑1 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTc                                                                     >  1:1317952/1‑36 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGc                                                           <  1:318645/46‑1 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGc                                                           <  1:7086343/46‑1 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGa                                                  >  1:1238355/1‑55 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACg                                  >  1:5853160/1‑71 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACg                                  >  1:1065879/1‑71 (MQ=255)
                                    agaTTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACg                                  >  1:4412035/1‑71 (MQ=255)
                                       ttGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCt                               >  1:2257584/1‑71 (MQ=255)
                                        tGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTg                              >  1:7672307/1‑71 (MQ=255)
                                          ttttACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAAcc                                             <  1:1675617/54‑1 (MQ=255)
                                          ttttACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTc                            >  1:883312/1‑71 (MQ=255)
                                          ttttACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTc                            >  1:2206448/1‑71 (MQ=255)
                                          ttttACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTc                            <  1:6207334/71‑1 (MQ=255)
                                            ttACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          >  1:7323207/1‑71 (MQ=255)
                                             tACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          <  1:1161142/70‑1 (MQ=255)
                                             tACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          <  1:756771/70‑1 (MQ=255)
                                             tACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          <  1:5724549/70‑1 (MQ=255)
                                              aCCTGCTGGGGAGTGGCACGCTATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTgg                         <  1:3325704/70‑1 (MQ=255)
                                              aCCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          <  1:5094912/69‑1 (MQ=255)
                                               ccTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCa                       <  1:6528276/71‑1 (MQ=255)
                                               ccTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCa                       <  1:362831/71‑1 (MQ=255)
                                                ctgctgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGt                             <  1:6613039/64‑1 (MQ=255)
                                                ctgctgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAg                      >  1:5038890/1‑71 (MQ=255)
                                                 tgctgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTc                            <  1:1725612/64‑1 (MQ=255)
                                                   ctgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGa                   <  1:6149224/71‑1 (MQ=255)
                                                    tgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGa                                                  <  1:7749815/39‑1 (MQ=255)
                                                    tgGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGAc                  >  1:5461095/1‑71 (MQ=255)
                                                     ggggAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCa                       <  1:7150982/65‑1 (MQ=255)
                                                      gggAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGAc                  <  1:7855066/69‑1 (MQ=255)
                                                      gggAATGGCGTGCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGAc                  <  1:110766/69‑1 (MQ=38)
                                                      gggAATGGCGTGCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGAc                  <  1:3169231/69‑1 (MQ=38)
                                                       ggAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTg               <  1:1477304/71‑1 (MQ=255)
                                                        gAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAAcc                                             <  1:789456/40‑1 (MQ=255)
                                                         aGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGca             >  1:6288287/1‑71 (MQ=255)
                                                          gTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGca             >  1:6627056/1‑70 (MQ=255)
                                                           tGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGAt                                      <  1:6836202/44‑1 (MQ=255)
                                                           tGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGcaca           <  1:1719632/71‑1 (MQ=255)
                                                           tGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGcaca           <  1:3601144/71‑1 (MQ=255)
                                                           tGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGcaca           <  1:2169638/71‑1 (MQ=255)
                                                             gCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTAc                                   >  1:4615087/1‑45 (MQ=255)
                                                             gCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAggg                    >  1:7182959/1‑60 (MQ=255)
                                                             gCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAggg                    >  1:1924364/1‑60 (MQ=255)
                                                             gCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAggg                    >  1:7173559/1‑60 (MQ=255)
                                                              cACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAg                                          >  1:7877001/1‑37 (MQ=255)
                                                              cACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGt                                         >  1:5530073/1‑38 (MQ=255)
                                                              cACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCa                       >  1:1319621/1‑56 (MQ=255)
                                                              cACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACgg        <  1:2247559/71‑1 (MQ=255)
                                                                 gCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGAtt     <  1:4904224/71‑1 (MQ=255)
                                                                 gCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGAtt     >  1:5898155/1‑71 (MQ=255)
                                                                 gCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          >  1:7486637/1‑50 (MQ=255)
                                                                 gCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          >  1:2010340/1‑50 (MQ=255)
                                                                 gCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          >  1:1826812/1‑50 (MQ=255)
                                                                 gCCACTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTg                          >  1:1724506/1‑50 (MQ=255)
                                                                   cATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGa       <  1:292549/67‑1 (MQ=255)
                                                                   cATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGa       <  1:3144533/67‑1 (MQ=255)
                                                                   cATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGATtgt   <  1:6461994/71‑1 (MQ=255)
                                                                    aTTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGATtgtg  >  1:1815192/1‑71 (MQ=255)
                                                                     |                                                                     
TTTTGATCGCAATTAACGCCACGCTTGAGGTAACAGAGATTGTTTTACCTGCTGGGGAGTGGCACGCCATTCCCCCATTCGCTGGAGAGGATAACCCAGTGATTACGGCTGTCTGGCAGGGACCTGCACACGGATTGTG  >  AP009048/4070915‑4071053

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