Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,029,577 C→A 18.2% L13F (TTG→TTT ugd ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,029,5770C→A18.2% 2.4 ‑3.2 0.03739ugdUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  new base (1/7):  ref base (17/12):  total (18/19)
Fisher's exact test strand distribution p-value = 4.22e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACCTCATGATTTTGTGCGATTAGAAGCCCGTTTGACAAGCCGACATAGCCAGTACCGGAAATGGTGATT  >  REL606/2029542‑2029610
                                   |                                 
aCCTCATGATTTTGTGCGATTAGAAGCCCGTTTGAc                                   <  1:1497253/36‑1 (MQ=255)
 ccTCAGGATTTTGTGCGATTAGAAGCCCGTTTGACa                                  <  1:2904400/36‑1 (MQ=255)
  cTCATGATTTTGTGCGATTAGAAGCCCGTTTGACaa                                 <  1:249545/36‑1 (MQ=255)
    cATGATTTTGTGCGATTAGAAGCCCGTTTGACAAGc                               <  1:4957413/36‑1 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              >  1:5724439/1‑36 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              >  1:5965981/1‑36 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              >  1:2276099/1‑36 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              <  1:726865/36‑1 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              >  1:6444873/1‑36 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              <  1:6697501/36‑1 (MQ=255)
     aTGATTTTGTGCGATTAGAAGCCCGTTTGACAAGcc                              <  1:4539244/36‑1 (MQ=255)
      tGATTTTGTGCGATTAGAAGCCCGTTTGACAAGCCg                             >  1:4483380/1‑36 (MQ=255)
       gATTTTGTGCGATTAGAAGCCCGTTTGACAAGCCGa                            >  1:4032192/1‑36 (MQ=255)
        aTTTTGTGCGATTAGAAGCCCGTTTGACAAGCCGAc                           <  1:576194/36‑1 (MQ=255)
         ttttGTGAGATTAGAAGCCCGTTTGACAAGCCGACa                          <  1:795770/36‑1 (MQ=255)
          tttGTGCGATTAGAAGCCCGTTGGACAAGCCGACAt                         <  1:4041213/36‑1 (MQ=255)
           ttGTGCGATTAGAAGCCCGTTTGAAAAGCCGACATa                        >  1:3137862/1‑36 (MQ=255)
            tgtgCGATTAGAAGCCCGTTTGACAAGCCGACATAg                       >  1:1433358/1‑36 (MQ=255)
             gtgCGATTAGAAGCCCGTTTGACAAGCCGACATAGc                      <  1:6471339/36‑1 (MQ=255)
             gtgCGATTAGAAGCCCGTTTGACAAGCCGACATAGc                      <  1:639665/36‑1 (MQ=255)
             gtgCGATTAGAAGCCCGTTTGACAAGCCGACATAGc                      <  1:5593071/36‑1 (MQ=255)
               gCGATTAGAAGCCCGTTTGACAAGCCGACATAGCCa                    >  1:1059807/1‑36 (MQ=255)
                   ttAGAAGCCCTTTTGAAAAGCCGACATAGCCAGTAc                <  1:511726/36‑1 (MQ=37)
                    tAGAAGCCCGTTTGACAAGTCGACATAGCCAGTAcc               >  1:2690501/1‑36 (MQ=255)
                    tAGAAGCCCGTTTGACAAACCGACATAGCCAGTAcc               >  1:4172471/1‑36 (MQ=255)
                    tAGAAGCCCGTTTGAAAAGCCGACATAGCCAGTAcc               <  1:7399829/36‑1 (MQ=255)
                     aGAAGCCCGTTTGACAAGCCGACATAGCCAGTACCg              >  1:5116327/1‑36 (MQ=255)
                     aGAAGCCCGTTTGACAAGCCGACATAGCCAGTACCg              >  1:1393433/1‑36 (MQ=255)
                       aaccccTTTAGAAAAGCCGACATAGCCAGTACCGGa            <  1:6819663/33‑1 (MQ=16)
                       aaGCCCTTTTGAAAAGCCGAAATAGCCAGTACCGGa            <  1:879354/36‑1 (MQ=25)
                        aGCCCGTTTGAAAAGCCGACATAGCCAGTACCGGaa           <  1:3608563/36‑1 (MQ=255)
                         ccccTTTTAAAAAGCAGACATAGCCAGTACCGGaaa          <  1:1803383/35‑1 (MQ=12)
                             gTTTGACAAGCCGACATAGCCAGTACCGGAAATGGt      >  1:5581791/1‑36 (MQ=255)
                             gTTTGACAAGCCGACATAGCCAGTACCGGAAATGGt      >  1:67638/1‑36 (MQ=255)
                              tttGACAAGCCGATATAGCCAGTACCGGAAATGGTg     >  1:4381303/1‑36 (MQ=255)
                                tGACAAGCCGACATAGCCAGTACCGGAAATGGTGAt   >  1:525500/1‑36 (MQ=255)
                                tGAAAAGCCGAAATAGCCAGTACCGGAAATGGTGAt   <  1:2277532/36‑1 (MQ=25)
                                 gACAAGCCGACATAGCCAGTACCGGAAATGGTGAtt  >  1:5502268/1‑36 (MQ=255)
                                   |                                 
ACCTCATGATTTTGTGCGATTAGAAGCCCGTTTGACAAGCCGACATAGCCAGTACCGGAAATGGTGATT  >  REL606/2029542‑2029610

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG