Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 514,834 C→T 8.5% L261F (CTC→TTC)  glxK → glycerate kinase II

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606514,8340C→T 8.5% 2.0 ‑2.8 ‑0.158781glxKglycerate kinase II
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (28/25):  total (33/25)
Fisher's exact test strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.22e-01

GCGGGCGGTATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGATATTAAACCGGGTATTGAGATTGT  >  REL606/514801‑514868
                                 |                                  
gcgggcggTATGGGTATTGCGGCGGCGGTATTtctc                                  <  1:5774302/36‑1 (MQ=255)
 cgggcggTATGGGTATTGCGGCGGCGGTATTTCTCa                                 <  1:6370018/36‑1 (MQ=255)
   gggggTATTGGTATTGCGGCGGCGGTATTTCTCAAt                               <  1:3045726/33‑1 (MQ=37)
   ggcggTATGGGTATTGCGGCGGCGGTATTTCTCAAt                               <  1:3814013/36‑1 (MQ=255)
    gcggTATGGGTATTGCGGCGGCGGTATTTCTCAATg                              >  1:8275344/1‑36 (MQ=255)
     cggTATGGGTATTGCGGCGGCGGTATTTCTTAAttg                             >  1:7726243/1‑34 (MQ=38)
     cggTATGGGTATTGCGGCGGCGGTATTTCTCAATGc                             <  1:2648353/36‑1 (MQ=255)
     cggTATGGGTATTGCGGCGGCGGTATTTCTCAATGc                             <  1:2292540/36‑1 (MQ=255)
      ggTATGGGTATTGCGGCGGCGGTATTTCTCAATTCg                            >  1:5550432/1‑36 (MQ=255)
      ggTATGGGTATTGCGGCGGCGGTATTTCTCAATGCg                            >  1:2280981/1‑36 (MQ=255)
      ggTATGGGTATTGCGGCGGCGGTATTTCTCAATGCg                            <  1:1987455/36‑1 (MQ=255)
      ggTATGGGTATTGCGGCGGCGGGATTTCTCAATGCg                            >  1:3227747/1‑36 (MQ=255)
       gTATGGGTATTGCGGCGGCGGTATTTTTTAATGGgg                           >  1:4684232/1‑36 (MQ=25)
       gTATGGGTATTGCGGCGGCGGTATTTCTCAATGCgg                           <  1:4721525/36‑1 (MQ=255)
        tATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGa                          <  1:4819191/36‑1 (MQ=255)
        tATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGa                          <  1:320737/36‑1 (MQ=255)
         aTGGGTATTGCGTCGGCGGTATTTCTCAATGCGGat                         <  1:2782473/36‑1 (MQ=255)
         aTGGGTATTGCGGCGGCGGTATTTCTCAATGCGGat                         <  1:4076479/36‑1 (MQ=255)
         aTGGGTATTGCGGCGGCGGTATTTCTCAATGCGGat                         >  1:1581624/1‑36 (MQ=255)
           gggTATTGCGGCGGCGGTATTTTTCCATGCGGgtat                       >  1:2108731/1‑36 (MQ=25)
           gggTATAGCGGCGGCTGTATTTCTCAATGCGGatat                       <  1:7165100/36‑1 (MQ=37)
            ggTATTGCGGCGGCGGTATTTCTCAATTGGGGTAtt                      >  1:5890276/1‑36 (MQ=25)
            ggTATTGCGGCGGCGGTATTTCTCAATGCGGGTAtt                      >  1:6832945/1‑36 (MQ=255)
             gTATTGCGGCGGCGGTTGTTCTCAATGCGGATAttt                     >  1:276318/1‑35 (MQ=37)
              tATTGCGGCGGCGGTATTTCTCAATGCGGATATTaa                    <  1:2235671/36‑1 (MQ=255)
               aTTGCGGCGGCGGTATTTTTCAATGCGGATATTaaa                   >  1:4162156/1‑36 (MQ=255)
               aTTGCGGCGGCGGTATTTCTCAATGCGGATATTaaa                   <  1:8301130/36‑1 (MQ=255)
               aTTGCGGCGGCGGGTTTTTTCCATTCGGGTTTTaaa                   >  1:6149926/1‑36 (MQ=2)
                ttGCGGCGGCGGTATTTTTCCATGCGGATATTTAAc                  >  1:1399304/1‑36 (MQ=25)
                ttGCGGCGGCGGTATTTCTCAATTCGGATATTAAAc                  >  1:157324/1‑36 (MQ=255)
                 tGCGGCGGCGGGCTTTCTTTAAGCGGGTATTaaatc                 >  1:6626833/1‑34 (MQ=2)
                  gcggcggcgGTATTTCTCAATGCGGATATTAAACCg                >  1:6071001/1‑36 (MQ=255)
                   cggcggcggTATTTCTCAATGCGGATATTAAACCgg               <  1:2168325/36‑1 (MQ=255)
                    ggcggcggTATTTCTCAATGCGGATATTAAACCggg              >  1:3330622/1‑36 (MQ=255)
                    ggcggcggTATTTCTCAATGCGGATATTAAACCggg              <  1:6544796/36‑1 (MQ=255)
                     ccggcggTATTTCTCAATTCGGATATTAAACCGGGt             <  1:5445857/35‑1 (MQ=38)
                       ggcggTATTTCTGCAAGGGGGTATTAAATCCGGTAt           >  1:7597303/1‑36 (MQ=0)
                       ggcggTATTTCTCAATGCGGATATTAAACCGGGTAt           <  1:2966097/36‑1 (MQ=255)
                        gcggTATTTCTCAATGCGGATATTAAACCGGGTAtt          >  1:1347523/1‑36 (MQ=255)
                         cggTATTTCTCAATGCGGATATTAAACCGGGTATTg         >  1:7465894/1‑36 (MQ=255)
                         cggTATTTCTCAATGCGGATATTAAACCGGGTATTg         >  1:6784100/1‑36 (MQ=255)
                          ggTATTTCTCAATGCGGATATTAAACCGGGTATTga        >  1:3730570/1‑36 (MQ=255)
                          ggTATTTCTCAATGCGGATATTAAACCGGGTATTga        <  1:6099436/36‑1 (MQ=255)
                          ggTATTTCTCAATGCGGATATTAAACCGGGTATTga        <  1:227335/36‑1 (MQ=255)
                          ggTATTTCTCAATGCGGATATTAAACCGGGTATTga        >  1:678001/1‑36 (MQ=255)
                           gTATTTCTCAATGCGGATATTAAACCGGGTATTgag       <  1:3778993/36‑1 (MQ=255)
                           gTATTTCTCAATGCGGATATTAAACCGGGTATTgag       >  1:2282462/1‑36 (MQ=255)
                            tATTTCTCAATGCGGATATTAAACCGGGTATTgaga      >  1:4000279/1‑36 (MQ=255)
                            tATTTCTCAATGCGGATATTAAACCGGGTATTgaga      <  1:7552516/36‑1 (MQ=255)
                            tATTTCTCAATGCGGATATTAAACCGGGTATTgaga      >  1:7459460/1‑36 (MQ=255)
                            tATTTCTCAATGCGGATATTAAACCGGGTATTgaga      >  1:7122421/1‑36 (MQ=255)
                            tATTTCTCAATGCGGATATTAAACCGGGTATTgaga      <  1:2852090/36‑1 (MQ=255)
                             aTTTCTCAATGCGGATATTAAACCGGGTATTGAGAt     >  1:4929196/1‑36 (MQ=255)
                             aTTTCTCAATGCGGATATTAAACCGGGTATTGAGAt     <  1:8165132/36‑1 (MQ=255)
                             aTTTCTCAATGCGGATATTAAACCGGGCATTGCGAt     >  1:3547870/1‑36 (MQ=37)
                             aTTGCTCAATGCGGGTATAAAAAAAGGTATTGAGAt     >  1:6065355/1‑36 (MQ=9)
                              tttCTCAATGCGGATATTAAACCGGGTGTTGAGAtt    >  1:1891436/1‑36 (MQ=255)
                              tttCTCAATGCGGATATTAAACCGGGTATTGAGAtt    <  1:4994932/36‑1 (MQ=255)
                               ttCTCAATGCGGATATTAAACCGGGTATTGAGATtg   >  1:4637669/1‑36 (MQ=255)
                                tctcAATGCGGATATTAAACCGGGTATTGAGATtgt  <  1:3294901/36‑1 (MQ=255)
                                tctcAATGCGGATATTAAACCGGGTATTGAGATtgt  <  1:3130092/36‑1 (MQ=255)
                                 |                                  
GCGGGCGGTATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGATATTAAACCGGGTATTGAGATTGT  >  REL606/514801‑514868

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