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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606514,8340C→T30.3% 2.2 0.0 0.022781glxKglycerate kinase II
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (14/1):  total (21/1)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGGCGGTATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGATATTAAACCGGGTATTGAGATTG  >  REL606/514803‑514867
                               |                                 
gggcggTATGGGTATTGCGGCGGCGGGTTTtctcca                               >  1:4867774/1‑34 (MQ=25)
   cggTATGGGTATTGCGGCGGCGGTTTTTTTCAAttg                            >  1:5307831/1‑34 (MQ=25)
     gTATGGGTATTGCGGAGGGGGTATTTTTTAAAGCgg                          >  1:893531/1‑36 (MQ=16)
      tATGGGTATTGCGGCGGCGGTATTTTTCAATGCggg                         >  1:142614/1‑35 (MQ=255)
        tGGGTATTGCGGCGGCGGTATTTTTCAATGCGGgta                       >  1:1399657/1‑36 (MQ=37)
         gggTATTGCGGCGGCGGTATTTCTCAATGCGGatat                      >  1:26328/1‑36 (MQ=255)
          ggTATTGCGGCGGCGGTATTTCTCAATGCGGATAtt                     <  1:2289369/36‑1 (MQ=255)
            tATTGCGGCGGCGGCATTTCCTCATGCTGGTATTaa                   >  1:867941/1‑36 (MQ=9)
             aTTGCGGCGGCGGTATTTTTCCTTGCGGATATTata                  >  1:5607707/1‑34 (MQ=21)
             aTTGCGGCGGCGGTATTTCTCAATGCGGATATTaaa                  >  1:857785/1‑36 (MQ=255)
             aTTGCGGCGGCGGTATTTCTCAATGCGGATATTaaa                  >  1:5811101/1‑36 (MQ=255)
             aTTGCGGCGGCGGTATTTCTCAATGCGGATATTaaa                  >  1:5291966/1‑36 (MQ=255)
               tGCGGCGGCGGGATTTTTCCAAGCGGGTATTAAAcc                >  1:3654961/1‑36 (MQ=12)
                  ggcggcggTATTTTTCAATGCGGATATTAAACCggg             >  1:3235430/1‑36 (MQ=255)
                      gcggTATTTCTCAATGCGGATATTAAACCGGGTAtt         >  1:3799413/1‑36 (MQ=255)
                      gcggTATTTCTCAATGCGGATATTAAACCGGGGAtt         >  1:2437456/1‑36 (MQ=255)
                       cggTATTTCTCAATGCGGATATTAAACCGGGTATTg        >  1:1417053/1‑36 (MQ=255)
                        ggTATTTCTCAATGCGGATATTAAACCGGGTATTga       >  1:2979662/1‑36 (MQ=255)
                        ggTATTTCTCAATGCGGATATTAAACCGGGTATTga       >  1:1798797/1‑36 (MQ=255)
                          tATTTCTCAATGCGGATATTAAACCGGGTCTTgcga     >  1:1030740/1‑36 (MQ=38)
                            tttCTCAATGCGGATATTAAACCGGGTATTGAGAtt   >  1:2535569/1‑36 (MQ=255)
                             ttCTCAATGCGGATATTAAACCGGGTATTGAGATtg  >  1:3896372/1‑36 (MQ=255)
                             ttCTCAATGCGGATATTAAACCGGGTATTGAGATtg  >  1:4624531/1‑36 (MQ=255)
                               |                                 
GGGCGGTATGGGTATTGCGGCGGCGGTATTTCTCAATGCGGATATTAAACCGGGTATTGAGATTG  >  REL606/514803‑514867

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG