Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,478,544 G→T 19.1% Q43K (CAA→AAA)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,478,5440G→T19.1% 2.9 ‑3.1 0.036127malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  new base (7/1):  ref base (12/16):  total (19/17)
Fisher's exact test strand distribution p-value = 4.38e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGG  >  REL606/3478515‑3478579
                             |                                   
cGTTACCGCCACTTTCGTGGCGGTACGTTGCTGcct                               >  1:7499048/1‑34 (MQ=38)
    aCCGCCACTTTCGTGCCGGTACGTTGACGCATcgcg                           >  1:844336/1‑36 (MQ=38)
      cGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTc                         <  1:2998008/36‑1 (MQ=255)
      cGCCACTTTCGTGGCGGGACCTTTATGTATTGCGTc                         >  1:5222183/1‑36 (MQ=16)
        ccACTTTCGTGGCGGTACGTTTATTCATTGCGTCaa                       >  1:1676182/1‑36 (MQ=25)
        ccACTTTCGTGGCGGTACGTTTATTCATCGCGTCaa                       >  1:5147307/1‑36 (MQ=37)
         cACTTTCGTGGCGGTACGTTTATGCATCGCGTCAAg                      >  1:1819185/1‑36 (MQ=255)
         cACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAg                      <  1:5715699/36‑1 (MQ=255)
          aCTTTCGTGGCGGTACTTTTATGCATTGCGTCAAGc                     >  1:6998998/1‑36 (MQ=25)
          aCTTGCGTGGCGGTACGTTGATGCATCGCGTCAAGc                     <  1:4312104/36‑1 (MQ=255)
             ttCGTGGCGGTACGTTTATGCATCTCGTCAAGCaaa                  >  1:246898/1‑36 (MQ=37)
             ttCGTGGCGGTACGTTTATGCATCGCTTCAAGCaaa                  <  1:1812457/36‑1 (MQ=37)
                 tGGCGGTATGTTGATGCATCGCGTCAAGCAAACGcc              <  1:7167133/36‑1 (MQ=255)
                 tGGCGGTACGTTTATTCATCGCGTCAAGCCAACGcc              >  1:4269612/1‑36 (MQ=25)
                  ggCGGTACGTTGATGAATCGCGTCAAGCAAACGCCg             >  1:7065008/1‑36 (MQ=255)
                    cGGTACGTTGATGCATCGCGTCAAGCAAACGCCGtt           >  1:4648379/1‑36 (MQ=255)
                     ggTACGTTGATGCATCGCGTCAAGCAAACGCCGttt          >  1:1288636/1‑36 (MQ=255)
                     ggTACGTTGATGCATCGCGTCAAGCAAACGCCGttt          <  1:6605010/36‑1 (MQ=255)
                     ggTACGTTGATGCATCGCGTCAAGCAAACGCCGttt          <  1:3119456/36‑1 (MQ=255)
                     ggTACGTTGATGCATCGCGTCAAGCAAACGCCGttt          <  1:6142513/36‑1 (MQ=255)
                      gTACGTTGATGCATCTCATCACGCAAACGCCTTTTg         <  1:4756849/36‑1 (MQ=16)
                      gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         <  1:4579587/36‑1 (MQ=255)
                      gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         <  1:186235/36‑1 (MQ=255)
                      gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         <  1:1436792/36‑1 (MQ=255)
                      gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         >  1:2010202/1‑36 (MQ=255)
                       tACGTTGATGCATCGCGTCAAGCAAACGCCGTTTgg        <  1:463603/36‑1 (MQ=255)
                       tACGTTGATGCATCGCGTCAAGCAAACGCCGTTTgg        >  1:1660232/1‑36 (MQ=255)
                        aCGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGt       >  1:6703426/1‑36 (MQ=255)
                         cGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGtt      <  1:1190002/36‑1 (MQ=255)
                         cGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGtt      >  1:1021671/1‑36 (MQ=255)
                          gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     <  1:3879919/36‑1 (MQ=255)
                          gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     <  1:6128312/36‑1 (MQ=255)
                          gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     >  1:3629278/1‑36 (MQ=255)
                          gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     >  1:2878065/1‑36 (MQ=255)
                           ttGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTc    <  1:2371965/36‑1 (MQ=255)
                            tGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCg   >  1:5140085/1‑36 (MQ=255)
                             gATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCgg  >  1:5214178/1‑36 (MQ=255)
                             |                                   
CGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGG  >  REL606/3478515‑3478579

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG