Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 507,302 T→A 9.1% F407L (TTT→TTA gcl → glyoxylate carboligase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606507,3020T→A 9.1% 2.5 ‑4.2 ‑0.0331221gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (24/6):  total (24/9)
Fisher's exact test strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01

TCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGG  >  REL606/507267‑507337
                                   |                                   
tCACAACTCTCTGCGGCTCACATGCTGCATTTCttt                                     <  1:85825/36‑1 (MQ=16)
tCACAAATCGCTGCGGCACAAATGCTGCATGTCttt                                     <  1:6035409/36‑1 (MQ=255)
tCACAAATCGCTGCGGCACAAATGCTGCATGTCttt                                     <  1:2022610/36‑1 (MQ=255)
 cacaAATCGCTGCGGCACAAATGCTGCATGTCTTTa                                    <  1:7027623/36‑1 (MQ=255)
   caAATCGCTGCGGCACAAATGCTGCATGTCTTTaaa                                  <  1:3221668/36‑1 (MQ=255)
     aaTCGCTGCGGCACAAATGCTGCATGTCTTTAAAGa                                >  1:1411958/1‑36 (MQ=255)
     aaTCACTGCGGCACAAATGCTGCATGTCTTTAAAGa                                <  1:7130274/36‑1 (MQ=255)
       tCGTTGCGGCACAAATGCTGCATGTCTTTAAAGAcc                              <  1:3750065/36‑1 (MQ=255)
       tCGCTGCGGCACAAATGCTGCATGTCTTTAAAGAcc                              >  1:8299596/1‑36 (MQ=255)
       tCGCTGCGGCACAAATGCTGCATGTCTTTAAAGAcc                              >  1:3401483/1‑36 (MQ=255)
        cGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCg                             <  1:5756581/36‑1 (MQ=255)
        cGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCg                             >  1:3010972/1‑36 (MQ=255)
        cGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCg                             <  1:3429317/36‑1 (MQ=255)
         gCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGc                            >  1:5658582/1‑36 (MQ=255)
         gCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGc                            >  1:3564760/1‑36 (MQ=255)
         gCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGc                            >  1:2993650/1‑36 (MQ=255)
          cTGCGGCACAAATGCTGCATGTCTTTAAAGACCGcc                           >  1:1644705/1‑36 (MQ=255)
            gCGGCACAAATGCTGCATGTCTTTAAAGACCGCCAc                         >  1:3073448/1‑36 (MQ=255)
             cGGCACAAATGCTGCATGTCTTTAAAGACCGCCACt                        >  1:2759572/1‑36 (MQ=255)
                cacaAATGCTGCATGTCTTTAAAGACCGCCACTGGa                     >  1:1633370/1‑36 (MQ=255)
                cacaAATGCTGCATGTCTTTAAAGACCGCCACTGGa                     >  1:4813963/1‑36 (MQ=255)
                 acaAATGCTGCATGTCTTTAAAGACCGCCACTGGAt                    >  1:6387819/1‑36 (MQ=255)
                 acaAATGCTGCATGTCTTTAAAGACCGCCACTGGAt                    >  1:7167072/1‑36 (MQ=255)
                 aaaaaTGCTGAATTTTTTAAAAGCCCGCCACTGGAt                    <  1:4916758/34‑1 (MQ=9)
                  caAATGCTGCATGTCTTTAAAGACCGCCACTGGATc                   >  1:6612158/1‑36 (MQ=255)
                  aaaaTGCTGCATTTATAAAAAGACCGCCACTGGATc                   <  1:7274019/35‑1 (MQ=16)
                  aaaaTCCTGCTTTTCTAAAAAGCCCGCCACTGGATc                   <  1:1193894/35‑1 (MQ=9)
                      tgctgcATGTCTTTAAAGACCGCCACTGGATCAACt               >  1:5129081/1‑36 (MQ=255)
                      tgctgcATGTCTTTAAAGACCGCCACTGGATCAACt               >  1:7132733/1‑36 (MQ=255)
                       gctgcATGTCTTTAAAGACCGCCACTGGATCAACtg              >  1:5164813/1‑36 (MQ=255)
                        ctgcATGTCTTTAAAGACCGCCACTGGATCAACtgt             >  1:8463763/1‑36 (MQ=255)
                          gcATGTCTTTAAAGACCGCCACTGGATCAACTGTgg           >  1:4781609/1‑36 (MQ=255)
                           cATGTCTTTAAAGACCGCCACTGGATCAACTGTGGt          >  1:6450716/1‑36 (MQ=255)
                            aTGTCTTTAAAGACCGTCACTGGATCAACTGTGGTc         >  1:186761/1‑36 (MQ=255)
                            aTGTCTTTAAAGACCGCCACTGGATCAACTGTGGTc         >  1:1662328/1‑36 (MQ=255)
                                cTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGc     >  1:1301937/1‑36 (MQ=255)
                                 tttAAAGACCGCAACTGGATCAACTGTGGTCAGGCt    >  1:3943747/1‑36 (MQ=255)
                                  ttAAAGACCGCCACTGGATCTACTGTGGTCAGGCTg   >  1:3257949/1‑36 (MQ=255)
                                   tAAAGACCGCCACTGGATCAACTGTGGTCAGGCTgg  >  1:2215936/1‑36 (MQ=255)
                                   |                                   
TCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGG  >  REL606/507267‑507337

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG