Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,806,030 C→A 17.9% E324* (GAA→TAA)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,806,0300C→A17.9% 3.0 ‑2.4 ‑0.037970astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (20/10):  total (22/15)
Fisher's exact test strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.87e-01

TGCCAGGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTC  >  REL606/1805996‑1806063
                                  |                                 
tgcCAGGCAGTAACCACCTGCTGTGCGGCCTGTTCa                                  >  1:2330623/1‑36 (MQ=255)
   cAGGCAGTAACCACCTGCTGTGCGGCCTGTTCAGaa                               >  1:3237899/1‑36 (MQ=255)
    aGGCAGTAACCACCTGCTGTGCGGCCTGTTCAGaaa                              >  1:3958845/1‑36 (MQ=255)
    aGGCAGTAACCACCTGCTGTGCGGCCTGTTCAGaaa                              <  1:2096694/36‑1 (MQ=255)
     ggCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAAt                             <  1:2500274/36‑1 (MQ=255)
      gCAGTAACGACCTGCTGTGCGGCCTGTTCAGAAATc                            <  1:7466310/36‑1 (MQ=255)
      gCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATc                            <  1:2087551/36‑1 (MQ=255)
      gCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATc                            >  1:6713708/1‑36 (MQ=255)
       cAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCa                           >  1:1164751/1‑36 (MQ=255)
        aGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAg                          >  1:3727028/1‑36 (MQ=255)
        aGTAACCACCTGCTGTGCGGCATGTTCAGAAATCAg                          >  1:1261862/1‑36 (MQ=255)
         gTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGc                         <  1:4485406/36‑1 (MQ=255)
          tAACCACCTGCTATGCGGCCTGTTCAGAAATCAGcc                        <  1:5916301/36‑1 (MQ=255)
           aaCCACCTGCTGTGCGGCCTGTTCAGAAATCAGccc                       >  1:7375884/1‑36 (MQ=255)
            accaccTGCTGGGCGGCCTGTTCAGAAATCAGCCCg                      >  1:1184609/1‑36 (MQ=255)
             ccaccTGCTGTGCGGCCTGTTCAGAAATCAGCCCGc                     >  1:2736736/1‑36 (MQ=255)
              caccTACTGTGCGGCCTTTTCAGAAATCAGCCCGcc                    <  1:173417/36‑1 (MQ=37)
               accTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCa                   >  1:1699876/1‑36 (MQ=255)
                ccTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCaa                  <  1:605102/36‑1 (MQ=255)
                  agctgTGCGGCTTGTTCAGAAAACAGCCCGCCAATa                <  1:6750578/35‑1 (MQ=37)
                   gctgTGCGGCCTGTTCAGAAATCAGCCCGCCAATaa               >  1:2415323/1‑36 (MQ=255)
                   gctgTGCGGCCTGTTAAAAAATCAGCCCGCCAATaa               <  1:5387135/36‑1 (MQ=37)
                    ctgTGCGGCCTGTTCAGAAATCATCCAGCCAATaaa              >  1:6031788/1‑36 (MQ=37)
                    ctgTGCGGCCTGTTCAGAAATCAGCCCGCCAATaaa              >  1:3700810/1‑36 (MQ=255)
                     tgtgCGGCTTGTTAAAAAATCACCCCGCCAATAAAc             <  1:3295211/36‑1 (MQ=19)
                     tgtgCGGCCTGTTCAGAAATCAGCCCGCCAATAAAc             >  1:4578362/1‑36 (MQ=255)
                     tgtgCGGCCTGTTAAGAAATCAGCCCGCCAATAAAc             <  1:91604/36‑1 (MQ=255)
                     tgtaCGGCATGTACAAAAACCAGCCCGCCAATAAAc             <  1:4077083/36‑1 (MQ=12)
                         cGGCCTGTTAAGAAATCAGCCCGCCAATAAACGGCt         <  1:2208523/36‑1 (MQ=255)
                         cGGCCTGTAAAAAAATCAGCCCGCCAATAAACGGCt         <  1:2706157/36‑1 (MQ=21)
                          ggCCTGTTAAGAAATCAGCCCGCCAATAAACGGCTg        >  1:515297/1‑36 (MQ=255)
                            ccTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCg      >  1:3620125/1‑36 (MQ=255)
                             cTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCgg     >  1:474497/1‑36 (MQ=255)
                             cTGTTAAGAAATCAGCCCGCCAATAAACGGCTGCgg     >  1:4192385/1‑36 (MQ=255)
                              tGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGt    >  1:5298339/1‑36 (MQ=255)
                               gTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGtt   >  1:2373794/1‑36 (MQ=255)
                                ttCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTc  >  1:1376873/1‑36 (MQ=255)
                                  |                                 
TGCCAGGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTC  >  REL606/1805996‑1806063

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