Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,174,628 C→G 14.0% G402G (GGC→GGG dld → D‑lactate dehydrogenase, FAD‑binding, NADH independent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,174,6280C→G14.0% 2.6 ‑3.9 0.0371206dldD‑lactate dehydrogenase, FAD‑binding, NADH independent
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (13/18):  total (19/18)
Fisher's exact test strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAAC  >  REL606/2174593‑2174662
                                   |                                  
catcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGc                                    <  1:3656101/36‑1 (MQ=255)
    tgctgTTAAAAATGGCGGGCGATGGAGTGGGCGAAg                                >  1:938437/1‑36 (MQ=255)
       tgTTAAAAATGGCGGGCGATGGCGTGGGGGGAGCCa                             >  1:4120257/1‑36 (MQ=25)
       tgTTAAAAATGGCGGGCGATGGCGTGGGGGGAGCCa                             >  1:7297327/1‑36 (MQ=25)
       tgTTAAAAATGGCGGGCGATGGCGTGGGGGAAGCCa                             >  1:3040718/1‑36 (MQ=255)
        gTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCaa                            >  1:5633308/1‑36 (MQ=255)
          tAAAAATGGCGGGCGATGGGGTGGGGGAAGGCAAAt                          >  1:3498300/1‑36 (MQ=21)
           aaaaaTGGCGGGCGATGGCGTGGGCGAAGCCAAATc                         <  1:1691009/36‑1 (MQ=255)
           aaaaaTGGCGGGCGATGGCGTGGGCGAAGCCAAATc                         <  1:395499/36‑1 (MQ=255)
           aaaaaTGGCGGGCGACAGCGTGGGGGAAGGCAAATc                         >  1:4547690/1‑36 (MQ=16)
           aaaaaTGGCGGGCCGTGGCGGGGGCGAAGCCATATc                         >  1:3634006/1‑36 (MQ=19)
            aaaaTGGCGGGCGATGGCGTGGGGGAAGCCAAATCg                        >  1:5611536/1‑36 (MQ=255)
             aaaTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGt                       >  1:3556047/1‑36 (MQ=255)
             aaaTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGt                       <  1:7505305/36‑1 (MQ=255)
              aaTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTg                      <  1:4091068/36‑1 (MQ=255)
               aTGGGGGGCGATGGCGTGGGCGAAGCCAAATCGTgg                     <  1:3692400/36‑1 (MQ=255)
                  gcgggcgATGGCGTGGGCGAAGCCAAATCGTGGCtg                  <  1:3136642/36‑1 (MQ=255)
                   cggtcgATGGCGTGGGCGAAGCCAAATCGTGGCtgg                 <  1:6258509/36‑1 (MQ=255)
                   cgggcgATGGCGTGGGCGAAGCCAAATCGTGGCtgg                 <  1:7326093/36‑1 (MQ=255)
                    gggcgATGGCGTGGGCGAAGCCAAATCGTGGCtggt                <  1:7647084/36‑1 (MQ=255)
                    gggcgATGGCGTGGGCGAAGCCAAATCGTGGCtggt                <  1:1280351/36‑1 (MQ=255)
                     ggcgATGGCGTGGGCGAAGCCAAATCGTGGCtggtg               <  1:4180907/36‑1 (MQ=255)
                      gcgCTGGCGTGGGCGAAGCCAAATCGTGGCtggtgg              <  1:4518231/36‑1 (MQ=255)
                       cgATGGCGTGGGCGAAGCCAAATCGTGGCtggtggg             >  1:1960390/1‑35 (MQ=255)
                        gATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGat            >  1:2715558/1‑36 (MQ=255)
                        gATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGat            <  1:3720756/36‑1 (MQ=255)
                         aTGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGatt           <  1:4536745/36‑1 (MQ=255)
                         aTGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGttt           >  1:970785/1‑36 (MQ=37)
                            gcttGGGCGAAGCCAAATCGTGGGTGGTGGattatt        <  1:905331/33‑1 (MQ=37)
                            gCGTGGGCGAAGCCAAATCGTGGCTGGTGGgttatt        >  1:1771082/1‑36 (MQ=255)
                            gCGTGGGCGAAGCCAAATCGTGGCTGGTGGattatt        <  1:156689/36‑1 (MQ=255)
                             cGTGGGCGAAGCCAAATCGTGGCTGGTGGATTAttt       >  1:2663569/1‑36 (MQ=255)
                              gTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTc      <  1:5157608/36‑1 (MQ=255)
                               tGGGCGAAGCCAAATCGTGGCTGGTGGTTTAttttt     >  1:1095143/1‑34 (MQ=37)
                               tGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCa     <  1:2607763/36‑1 (MQ=255)
                                 ggCGAAGCCAAATCGTGGCTGGTGGATTATTTCaaa   >  1:3412057/1‑36 (MQ=255)
                                 ggCGAAGCCAAATCGTGGCTGGTGGATTATTTCaaa   >  1:6801431/1‑36 (MQ=255)
                                 gacGAAGCCAAATCGTGGCTGGTGGATTATTTCaa  <  1:2115748/34‑2 (MQ=255)
                                  gCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAAc  >  1:5466148/1‑36 (MQ=255)
                                   |                                  
CATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAAC  >  REL606/2174593‑2174662

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