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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,478,5440G→T24.3% 5.6 ‑3.9 ‑0.046127malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  new base (12/0):  ref base (22/12):  total (34/12)
Fisher's exact test strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.40e-01

TGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGG  >  REL606/3478512‑3478579
                                |                                   
tGGCGTTACCGCCACTTTCGTGGCGGTACGTTGTTg                                  >  1:6785959/1‑36 (MQ=255)
       aCCGCCACTTTCGTGGCGGTACGTTTATTCATtgcg                           >  1:7906104/1‑36 (MQ=25)
       aCCGCCACTTTCGTGGCGGTACGTTGTTTCCTcgcg                           >  1:3630429/1‑36 (MQ=25)
        ccGCCACTTTCGTGGCGGTACGTTGATGCATCGCGt                          <  1:2007048/36‑1 (MQ=255)
         cGCCACTTTCGTGGCGGTACTGTGATTTATCGCCTc                         >  1:6198641/1‑36 (MQ=16)
         cGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTc                         <  1:1867633/36‑1 (MQ=255)
         cGCCACTTTCGTGGCGGTACCTTTATGCATCGCGTc                         >  1:4187145/1‑36 (MQ=37)
          gCCACTTTCGTGGCGGTTCGGTTGTGCATTGCCTcc                        >  1:3232837/1‑35 (MQ=9)
          gCCACTTTCGTGGCGGTACGTTTATGCATCGCGTcc                        >  1:2622404/1‑35 (MQ=255)
          gCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCa                        <  1:797486/36‑1 (MQ=255)
           ccACTTTCGTGGCGGGACGTTTATGTTTCGCGGCac                       >  1:4736337/1‑35 (MQ=12)
            cACTTTCGTGGCGGTACGTTTATGGATTGCGTCAAg                      >  1:6106819/1‑36 (MQ=25)
             aCTTTCGTGGCGGTACGTTTTTGCTTCGAGTCAGGc                     >  1:1070792/1‑36 (MQ=16)
             aCTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGc                     <  1:4999374/36‑1 (MQ=255)
              cTTTCGTGGCGGTACGTTTATTCATCGCGTCAAGCa                    >  1:7807465/1‑36 (MQ=37)
              cTTTCGTGGCGGTACGTTTATGCATCGCGTCAAGCa                    >  1:2507183/1‑36 (MQ=255)
              cTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCa                    <  1:2795255/36‑1 (MQ=255)
               tttCGTGGCGGTACGTTTATTCATCTCGTCAAGCaa                   >  1:1285547/1‑36 (MQ=25)
               tttCGTGGCGGTACGTTTATTCATCGCGTCAAGcca                   >  1:296983/1‑34 (MQ=25)
               tttCGTGGCGGTACGTTGATGCATCGCGTCCAGcca                   >  1:5926933/1‑34 (MQ=37)
               tttCGTGGCGGTACGTTGATGCATCGCGTCAAGCaa                   >  1:1043590/1‑36 (MQ=255)
                ttCGTGGCGGTACGTTGATGCATCGCGTCAAGCaaa                  <  1:5131087/36‑1 (MQ=255)
                  cGTGGCGGTACGTTTATGCATCGCGTCAAGCAAACg                >  1:4426644/1‑36 (MQ=255)
                  cGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACg                <  1:2166658/36‑1 (MQ=255)
                    tGGCGGTACGTTGATGCATCGCGTCAAGCAAACGcc              >  1:3937916/1‑36 (MQ=255)
                    tGGCGGTACGTTGATGCATCGCGTCAAGCAAACGcc              <  1:7189643/36‑1 (MQ=255)
                     ggCGGTACGTTGATGCATCGCGTCAAGCAAACGCCg             >  1:1952584/1‑36 (MQ=255)
                     ggCGGTACGTTGATGCATCCCGTCAAGCAAACGCCg             >  1:788171/1‑36 (MQ=255)
                      gCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGt            >  1:6767359/1‑36 (MQ=255)
                      gCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGt            >  1:5645175/1‑36 (MQ=255)
                       cGGTACGTTGATGCATCGCGTCAAGCAAACGCCGtt           >  1:958887/1‑36 (MQ=255)
                       cGGTACGTTGATGCATCGCGTCAAGCAAACGCCGtt           <  1:6650668/36‑1 (MQ=255)
                       cGGTACGTTGATGCATCGCGTCAAGCAAACGCCGtt           >  1:4169091/1‑36 (MQ=255)
                        ggTACGTTGATGCATCGCGTCAAGCAAACGCCGttt          <  1:4328053/36‑1 (MQ=255)
                         gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         >  1:6576094/1‑36 (MQ=255)
                         gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTg         >  1:4879062/1‑36 (MQ=255)
                          tACGTTGATGCATCGCGTCAAGCAAACGCCGTTTgg        >  1:5748754/1‑36 (MQ=255)
                          tACGTTGATGCATCGCGTCAAGCAAACGCCGTTTgg        <  1:8308848/36‑1 (MQ=255)
                           aCGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGt       >  1:4804460/1‑36 (MQ=255)
                             gTTGATGCATCGCGTCAAGCAAACGCCNNTTGGttt     >  1:1359125/1‑36 (MQ=255)
                             gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     >  1:5122302/1‑36 (MQ=255)
                             gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     >  1:3619337/1‑36 (MQ=255)
                             gTTGATGCATCGCGTCAAGCAAACGCCGTTTGGttt     >  1:8016454/1‑36 (MQ=255)
                              ttGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTc    >  1:1474004/1‑36 (MQ=255)
                              ttGATGCATCGCGNCAAGCANACGCCGTTTGGTTTc    >  1:1554638/1‑36 (MQ=40)
                               tGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCg   <  1:7223360/36‑1 (MQ=255)
                                gATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCgg  >  1:1955740/1‑36 (MQ=255)
                                |                                   
TGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGG  >  REL606/3478512‑3478579

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG