Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,138,473 G→T 10.3% Q52K (CAA→AAA)  grxB ← glutaredoxin 2 (Grx2)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,138,4730G→T10.3% 2.9 ‑3.8 ‑0.039154grxBglutaredoxin 2 (Grx2)
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (12/23):  total (16/23)
Fisher's exact test strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.77e-01

CATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGG  >  REL606/1138438‑1138508
                                   |                                   
cATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTg                                     <  1:6583219/36‑1 (MQ=255)
    cTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGa                                 <  1:2897320/36‑1 (MQ=255)
      ttCTCGCATATAGCGGCTGTCGTCTTTTTGCAGAAt                               <  1:5180421/36‑1 (MQ=37)
       tCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATg                              <  1:5695368/36‑1 (MQ=255)
         tGGCATATAGCGGCTGTCATCTTTTTGCAGAATggg                            <  1:7496377/36‑1 (MQ=255)
          ggCATATAGCGGCTGTGATTTTTTTTTAGAATtggg                           >  1:667261/1‑35 (MQ=9)
           gCATATAGCGGCTGTCATCTTTTTGCAGAATGGGaa                          <  1:1152715/36‑1 (MQ=255)
            cATATAGCGGCTGNCATCTTTTTTCAGAATGGGGAc                         >  1:4339112/1‑36 (MQ=37)
             atatAGCGGCTGTCATCTTTTTTCAGAATTGGAAcc                        >  1:3807710/1‑36 (MQ=37)
              tataGCGGCTGTCATCTTTTTGCAGAATGGGAACCt                       <  1:4439439/36‑1 (MQ=255)
               ataGCGGCTGTCATCTTTTTTCAGAATGGGAACCTg                      >  1:5330862/1‑36 (MQ=255)
                taGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGt                     <  1:6771404/36‑1 (MQ=255)
                taGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGt                     <  1:4292991/36‑1 (MQ=255)
                 aGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGtt                    <  1:4926611/36‑1 (MQ=255)
                  gCGGCTGTCATCTTTTTGGATAATGGGAACCTTttt                   <  1:1977377/36‑1 (MQ=37)
                  gCGGCTGTCATCTTTTTGCAGAATGGGAACCTGttt                   <  1:1691066/36‑1 (MQ=255)
                    ggCTGTCATCTTTTTGCAGAATGGGAACCTGttttt                 <  1:1467110/36‑1 (MQ=255)
                     gCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTg                >  1:2720215/1‑36 (MQ=255)
                     gCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTg                <  1:5142925/36‑1 (MQ=255)
                      cTGTCATCTTTTTGCTGAATGGGAACCTGTTTTTga               >  1:5147992/1‑36 (MQ=255)
                      cTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTga               <  1:5858629/36‑1 (MQ=255)
                       tGTCATCTTTTTGCAGAATGGGAACCTGTTTTTgac              <  1:4142029/36‑1 (MQ=255)
                       tGTCATCTTTTTGCAGAATGGGAACCTGTTTTTgac              <  1:6301628/36‑1 (MQ=255)
                        ttCATCTTTTTGCAGAATGGGAACCTGTTTTTgacc             <  1:1289649/35‑1 (MQ=255)
                        gTCATCTTTTTGCAGAATGGGAATCTGTTTTTgacc             >  1:6710429/1‑36 (MQ=255)
                        gTCATCTTTTTGCAGAATGGGAACCTGTTTTTgacc             >  1:3980983/1‑36 (MQ=255)
                         tCATCTTTTTGCAGAATGGGAACCTGTTTTTgaccg            <  1:6812280/36‑1 (MQ=255)
                           aTCTTTTTGCAGAATGGGAACCTGTTTTTgaccgac          <  1:1706440/36‑1 (MQ=255)
                            tCTTTTTGCAGAATGGGAACCTGTTTTTgaccgacc         <  1:925523/36‑1 (MQ=255)
                            tCTTTTTGCAGAATGGGAACCTGTTTTTgaccgacc         >  1:3671628/1‑36 (MQ=255)
                              tttttGCAGAATGGGAACCTGTTTTTGACCGACCAt       >  1:2259372/1‑36 (MQ=255)
                               ttttGCAGAATGGGATCCTGTTTTTGACCGACCATc      <  1:5200196/36‑1 (MQ=255)
                               ttttGCAGAATGGGAACCTGTTTTTGACCGACCATc      >  1:1587526/1‑36 (MQ=255)
                               ttttGCAGAATGGGAACCTGTTTTTGACCGACCATc      >  1:2746254/1‑36 (MQ=255)
                               ttttGCAGAATGGGAACCTGTTTTTGACCGACCATc      >  1:6882758/1‑36 (MQ=255)
                                tttGCAGAATGGGAACCTGTTTTTGACCGACCATcc     <  1:6701353/36‑1 (MQ=255)
                                tttGCAGAATGGGAACCTGTTTTTGACCGACCATcc     >  1:673843/1‑36 (MQ=255)
                                  tGCAGAATGGGAACCTGTTTTTGACCGACCATCCgg   >  1:6970270/1‑36 (MQ=255)
                                  tGCAGAATGGGAACCTGTTTTTGACCGACCATCCgg   <  1:3245282/36‑1 (MQ=255)
                                  tGCAGAATGGGAACCTGTTTTTGACCGACCATCCgg   >  1:878262/1‑36 (MQ=255)
                                   gCAGAATGGGAACCTGTTTTTGACCGACCATCCggg  >  1:1716768/1‑36 (MQ=255)
                                   gCAGAATGGGAACCTGTTTTTGACCGACCATCCggg  >  1:2807163/1‑36 (MQ=255)
                                   |                                   
CATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGG  >  REL606/1138438‑1138508

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG