Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 990,797 C→G 15.3% Q57E (CAG→GAG)  smtA → predicted S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606990,7970C→G15.3% 3.8 ‑3.5 0.042169smtApredicted S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (18/17):  total (25/17)
Fisher's exact test strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAG  >  REL606/990764‑990832
                                 |                                   
cTGCGTGTGCTGGATGCTGGCGGTGGGGAAGGGCAg                                   >  1:7633904/1‑36 (MQ=255)
 tGCGTGTGCTGGATGCTGGCGGTGGGGAAGGGGAGa                                  >  1:6160961/1‑36 (MQ=37)
  gCGTGTGCTGGATGCTGGCGGTGGGGAAGGGGAggg                                 >  1:3980852/1‑34 (MQ=25)
  gCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGAc                                 <  1:2875748/36‑1 (MQ=255)
   cgcgtgCTGTTTGCTGGCGGTGGAGAAGGGCAGAcc                                <  1:2721405/33‑1 (MQ=21)
     tgtgCTGGATGCTGGCGGTGGAGAAGGGCAGACCGc                              <  1:4986971/36‑1 (MQ=255)
     tgtgCTGGATGCTGGCGGTGGAGAAGGGCAGACCGc                              <  1:7447256/36‑1 (MQ=255)
      gtgCTGGATGCTGGCGGTGGGGAAGGGGAGACCGCa                             >  1:7516275/1‑36 (MQ=37)
      gtgCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCa                             <  1:4992189/36‑1 (MQ=255)
       tgCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCaa                            <  1:2721698/36‑1 (MQ=255)
       tgCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCaa                            >  1:5328180/1‑36 (MQ=255)
        gCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAAt                           <  1:5443142/36‑1 (MQ=255)
        gCTGGATGCTGGCGGGGGGGGAGGGCAGACCGGAAt                           >  1:5565788/1‑36 (MQ=19)
         cTGGATGCTGGCGGAGGGGGAGGGGAGAGCGGaaaa                          >  1:4269523/1‑34 (MQ=9)
          tGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCa                         <  1:5349543/36‑1 (MQ=255)
           ggATGCTGGCGGTGGGGAAGGGGAGACCGCAATCaa                        >  1:5138059/1‑36 (MQ=37)
           ggATGCTGGCGGTGGGGAAGGGCAGACCGCAATCaa                        >  1:1527827/1‑36 (MQ=255)
            gATGCTGGCGGTGGGGAAGGGCAGACCGCAATCaaa                       >  1:2790967/1‑36 (MQ=255)
            gATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaaa                       <  1:953723/36‑1 (MQ=255)
             aTGCTGGCGGTGGAGAAGGGCAGACCGCAATCaaaa                      <  1:7379301/36‑1 (MQ=255)
              tGCTGGCGGTGGGGAAGGGGAGACCGCAAACAAAAt                     >  1:3766697/1‑36 (MQ=25)
              tGCTGGCGGTGGAGAAGGGGAGAACGCAATCAAAAt                     >  1:3681063/1‑36 (MQ=37)
                cTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATgg                   <  1:3763646/36‑1 (MQ=255)
                 tGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                  >  1:1626870/1‑36 (MQ=255)
                 tGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                  <  1:7872747/36‑1 (MQ=255)
                 tGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                  <  1:4768258/36‑1 (MQ=255)
                   gCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCg                >  1:2171491/1‑36 (MQ=255)
                   gCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCg                <  1:6682047/36‑1 (MQ=255)
                   gCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCg                >  1:4659477/1‑36 (MQ=255)
                     ggTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAg              >  1:7130143/1‑36 (MQ=255)
                        ggAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGt           >  1:2785463/1‑36 (MQ=255)
                          agaAGGGCAGACCGCAATCAAAATGGCCGAGCGTgg         >  1:5708186/1‑36 (MQ=255)
                          agaAGGGCAGACCGCAATCAAAATGGCCGAGCGTgg         >  1:1921327/1‑36 (MQ=255)
                            aaGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGc       >  1:751440/1‑36 (MQ=255)
                            aaGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGc       >  1:6452252/1‑36 (MQ=255)
                             aGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCa      >  1:7573888/1‑36 (MQ=255)
                             aGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCa      >  1:5880123/1‑36 (MQ=255)
                             aGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCa      <  1:7790685/36‑1 (MQ=255)
                              gggCAGACCGCAATCAAAATGGCCGAGCGTGGGCAt     >  1:6808308/1‑36 (MQ=255)
                              gggCAGACCGCAATCAAAATGGCCGAGCGTGGGCAt     >  1:4274286/1‑36 (MQ=255)
                               ggCAGACCGCAATAAAAATGGCCGAGCGTGGGCATc    <  1:5308408/36‑1 (MQ=255)
                                gCAGACCGCAATCAAAAGGGCCGAGCGTGGGCATCa   <  1:2663974/36‑1 (MQ=255)
                                 cAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAg  >  1:5767853/1‑36 (MQ=255)
                                 cAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAg  <  1:1255801/36‑1 (MQ=255)
                                 |                                   
CTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAG  >  REL606/990764‑990832

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG