Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 804,332 A→G 32.2% E62G (GAA→GGA)  bioC → predicted methltransferase, enzyme of biotin synthesis

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606804,3320A→G32.2% 3.1 ‑4.0 ‑0.021185bioCpredicted methltransferase, enzyme of biotin synthesis
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (6/8):  total (13/8)
Fisher's exact test strand distribution p-value = 1.80e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.53e-01

TGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGC  >  REL606/804298‑804367
                                  |                                   
tgtgGACCTGGCTGGATGAGCCGCCACTGGCGGGaa                                    <  1:592105/36‑1 (MQ=255)
 gtgGACCTGGCTGGATGAGCCGCCACTGGCGGGAAc                                   >  1:5477400/1‑36 (MQ=255)
   ggACCTGGCTGGATGAGCCGCCACTGGCGGGGACgg                                 >  1:2838879/1‑35 (MQ=255)
    gACCTGGCTGGATGAGCCGCCACTGGCGGGGACggc                                >  1:695387/1‑34 (MQ=37)
     aCCTGGCTGGATGAGCCGCCACTGGCGGGGACggcc                               >  1:1106747/1‑33 (MQ=37)
     aCCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCa                               >  1:4852373/1‑36 (MQ=37)
          gctggATGAGCCGCCACTGGCGGGGACGTTAcgcgg                          >  1:2566373/1‑35 (MQ=37)
           ctggATGAGCCGCCACTGGCGGGAACGTCACGCGCa                         <  1:5022119/36‑1 (MQ=255)
              gATGAGCCGCCACTGGCGGGGACGTCACGCGCAggg                      >  1:3108037/1‑35 (MQ=255)
              gATGAGCCGCCACTGGCGGGGACGGCACGCGCCggg                      >  1:316283/1‑35 (MQ=25)
              gATGAGCCGCCACTGGCGGGAAAGTCCCGCGCAggg                      >  1:3922262/1‑35 (MQ=37)
                 gAGCCGCCCCTGGCGGGAACGTCACGCGCAGGTGAc                   <  1:3240707/36‑1 (MQ=255)
                   gccgccACTGGCGGGAACGTCACGCGGAGGGGACgg                 >  1:4420767/1‑36 (MQ=37)
                     cgccACTGGCGGGAACGTCACGCGCAGGTGACGGcc               <  1:4416190/36‑1 (MQ=255)
                       ccACTGGCGGTAACGTCACGCGCAGGTGACGGCCtt             <  1:2771689/36‑1 (MQ=255)
                        cACTGGCGGGAACGTCACCCGCAGGTGACGGCCTTa            <  1:5624744/36‑1 (MQ=255)
                           tGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAt         <  1:5610552/36‑1 (MQ=255)
                             gCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtct       >  1:2501174/1‑36 (MQ=255)
                             gCGGGAACGTCACGCGCAGGAGACGGCCTTAGAtct       >  1:1619256/1‑36 (MQ=255)
                               gggAACGTCACGCGCAGGTGACGGCCTTAGAtctct     >  1:4285032/1‑36 (MQ=255)
                                ggAACGTCACGCGCAGGTGACGGCCTTAGAtctctc    <  1:3559359/36‑1 (MQ=255)
                                 gAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcg   <  1:1670274/36‑1 (MQ=255)
                                  aaCGTCACGCGCCGTTGAGGGCCTTAGATCTCTcgc  <  1:2719936/36‑1 (MQ=37)
                                  aaCGTCACGCGCAGGTGACGGCCTTAGATCTCTcgc  <  1:5434026/36‑1 (MQ=255)
                                  aaCGTCACGCGCAGGTGACGGCCTTAGATCTCTcgc  >  1:1827022/1‑36 (MQ=255)
                                  |                                   
TGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGC  >  REL606/804298‑804367

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG