Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,733,358 T→C 8.7% Y127H (TAC→CAC)  bglB → cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,733,3580TC 8.7% 120.8 / 4.3 46Y127H (TAC→CAC) bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCGTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCATCG  >  AP009048/3733289‑3733418
                                                                     |                                                            
gCGTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTa                                                             <  1:2024847/71‑1 (MQ=255)
 cGTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTAc                                                            >  1:4060149/1‑71 (MQ=255)
  gTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTAc                                                            >  1:389862/1‑70 (MQ=255)
  gTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTAc                                                            >  1:2130883/1‑70 (MQ=255)
     ttACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCAtt                                                              >  1:3071559/1‑65 (MQ=255)
     ttACGATCGGCTGTTTGATAAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGa                                                          >  1:2769759/1‑69 (MQ=255)
        cGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGc                                                     >  1:1306601/1‑71 (MQ=255)
         gATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGcc                                                    <  1:799538/71‑1 (MQ=255)
               cTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATgg                                              <  1:3192661/71‑1 (MQ=255)
                tGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCa                                                   >  1:3799015/1‑65 (MQ=255)
                  tttGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGc                                            <  1:3470931/70‑1 (MQ=255)
                    tgatgaAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGc                                            >  1:521134/1‑68 (MQ=255)
                    tgatgaAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTgg                                         >  1:3905559/1‑71 (MQ=255)
                         aaaTGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGaaa                                    <  1:4127628/71‑1 (MQ=255)
                         aaaTGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGaaa                                    <  1:3778175/71‑1 (MQ=255)
                            tGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATg                                               >  1:4085106/1‑57 (MQ=255)
                            tGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAAc                                 <  1:3858332/71‑1 (MQ=255)
                              gcgcAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTa                               >  1:539749/1‑71 (MQ=255)
                              gcgcAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTa                               >  1:3205830/1‑71 (MQ=255)
                              gcgcAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTa                               >  1:349735/1‑71 (MQ=255)
                                gcAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTAcg                             <  1:2284390/71‑1 (MQ=255)
                                 cAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTAcgg                            <  1:431095/71‑1 (MQ=255)
                                   ggCGGGGATCAAGCCGCTGGTAACGTTATCCCAtt                                                              <  1:4046944/35‑1 (MQ=255)
                                   ggCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCa                                                   <  1:2967784/46‑1 (MQ=255)
                                    gCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGaaa                                                        >  1:1413077/1‑40 (MQ=255)
                                    gCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCt                                           >  1:1170603/1‑53 (MQ=255)
                                    gCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTAcggcgg                         <  1:3931567/71‑1 (MQ=255)
                                    gCGGGGATCAAGCCGCTGGTAACGTTATCCCATCACGa                                                          <  1:4503777/38‑1 (MQ=255)
                                     cGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGt                        <  1:3317076/71‑1 (MQ=255)
                                      ggggATCAAGCCGCTGGTAACGTTATCCCATTACGAAATg                                                      <  1:2241669/40‑1 (MQ=255)
                                      ggggATCAAGCCGCTGGTAACGTTATCCCATCACGAAATGCCATATg                                               <  1:391357/47‑1 (MQ=255)
                                      ggggATCAAGCCGCTGGTAACGTTATAACATTACGAAATg                                                      >  1:52443/1‑40 (MQ=255)
                                         gATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGc                                            <  1:1799309/47‑1 (MQ=255)
                                         gATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTAc                              >  1:2902811/1‑61 (MQ=255)
                                          aTCAAGCCGCTGGTAACGTTATCCCATCACGAAATGCCATATGGGCTGGTGATAAACTACGGCGGTTGGGc                   >  1:720790/1‑71 (MQ=255)
                                           tCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCt                  <  1:2292555/71‑1 (MQ=255)
                                           tCAAGCCGCTGGTAACGTTATCCCATCACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGc                   >  1:4379588/1‑70 (MQ=255)
                                              aGCCGCTGGTAAGGTTATCCCATTACGAAATGCCATATGGGCTGGTGaaaaa                                  >  1:1155024/1‑52 (MQ=255)
                                               gCCGCTGGTAACGTTATCCCATTACGAAATGCCata                                                 >  1:3706222/1‑36 (MQ=255)
                                               gCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATc              >  1:134352/1‑71 (MQ=255)
                                                ccGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCTGTTGGGCTAATCg             <  1:4057726/71‑1 (MQ=255)
                                                   cTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAg           <  1:3133239/70‑1 (MQ=255)
                                                       tAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGaaa                                    <  1:3675733/41‑1 (MQ=255)
                                                        aaCGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCa     >  1:4104870/1‑71 (MQ=255)
                                                           gTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCATc   >  1:4369642/1‑70 (MQ=255)
                                                           gTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCATCg  <  1:4150782/71‑1 (MQ=255)
                                                               tCCCATTACGAAATGCCATATGGGCTTGTGAAAAAc                                 <  1:98715/36‑1 (MQ=255)
                                                                 ccATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCATCg  <  1:1180820/65‑1 (MQ=255)
                                                                     |                                                            
GCGTTTTACGATCGGCTGTTTGATGAAATGGCGCAGGCGGGGATCAAGCCGCTGGTAACGTTATCCCATTACGAAATGCCATATGGGCTGGTGAAAAACTACGGCGGTTGGGCTAATCGAGCGGTCATCG  >  AP009048/3733289‑3733418

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