Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,594,079 G→T 15.3% L356L (CTC→CTA dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,594,0790G→T15.3% 2.3 ‑4.1 ‑0.0251068dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (6/15):  total (10/15)
Fisher's exact test strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.41e-01

AACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCC  >  REL606/1594050‑1594113
                             |                                  
aaCACCGTGTTTAATTCAAAATATGGCTGGAGCttg                              >  1:932512/1‑34 (MQ=25)
    ccGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCt                          <  1:2292036/36‑1 (MQ=255)
     cGTGTATAATTCAAAATATGGCTTGAGCTGCGCCtc                         <  1:5009500/36‑1 (MQ=255)
      gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCtca                        <  1:5821215/36‑1 (MQ=255)
         tttAATTCAAAATATGGCTTTAGCTTCTCCtcatca                     >  1:364143/1‑36 (MQ=25)
         tttAATTCAAAATATGGCTTTAGCTTCGCCtcatcc                     >  1:4438857/1‑35 (MQ=37)
         tttAATTCAAAATATGGCTTGAGCTGCGCCtcatca                     <  1:632118/36‑1 (MQ=255)
          ttAATTCAAAATATGGCTTGAGCTGCGCCTCATCaa                    <  1:614941/36‑1 (MQ=255)
             aTTCAAAATATGGCTTTTGCTGCGCCTCTTCAAGag                 >  1:5799384/1‑35 (MQ=21)
             aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGAt                 <  1:3314462/36‑1 (MQ=255)
             aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGAt                 <  1:3366699/36‑1 (MQ=255)
             aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGAt                 <  1:4545565/36‑1 (MQ=255)
              ttCAAAATATGGCTTTATCTGCGCCTCCTCAAGATc                >  1:2534345/1‑36 (MQ=25)
              ttCAAAATATGGCTTGAGCTGCGCCTCATCAAGATc                <  1:4343368/36‑1 (MQ=255)
              ttCAAAATATGGCTTGAGCTGCGCCTCATCAAGATc                <  1:852363/36‑1 (MQ=255)
                cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCaa              <  1:4668633/36‑1 (MQ=255)
                    atatGGCTTGAGCTGCGCCTCATCAAGATCAAAttt          <  1:3116137/36‑1 (MQ=255)
                    atatGGCTTGAGCTGCGCCTCATCAAGATCAAAttt          <  1:3915158/36‑1 (MQ=255)
                    atatGGCTTGAGCTGAGCCTCATCAAGATCAAAttt          <  1:5496953/36‑1 (MQ=255)
                     tatGGCTTGAGCTGGGGCTCAGTAAGATCAAAATtc         >  1:5455115/1‑36 (MQ=16)
                        ggCTTGAGCTGCGCCTCATCAAGATCAAATTtctac      >  1:5851983/1‑34 (MQ=255)
                         gCTTGAGCTGCGCCTCATCAAGATCAAATTTCTccc     >  1:108346/1‑36 (MQ=255)
                          cTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccg    >  1:1784403/1‑36 (MQ=255)
                           ttGAGCTGCGCCTCATCAAGCTCAAATTTCTCccgc   <  1:438200/36‑1 (MQ=255)
                            tGAGCTGCGCCTCATCAAGATCAAATTTCTCccgcc  >  1:1006115/1‑36 (MQ=255)
                             |                                  
AACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCC  >  REL606/1594050‑1594113

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG