Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,561,446 C→G 9.6% L1097L (CTG→CTC hsdR ← endonuclease R

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,561,4460C→G 9.6% 10.1 ‑0.4 ‑0.0623291hsdRendonuclease R
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (27/29):  total (29/33)
Fisher's exact test strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

TCGCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCA  >  REL606/4561411‑4561480
                                   |                                  
tCGCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTc                                    <  1:448500/36‑1 (MQ=255)
 cGCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTCa                                   <  1:1695093/36‑1 (MQ=255)
  gctgctCCAGTCATTTTCGCCCTTAATGCGCGTCAg                                  >  1:1713265/1‑36 (MQ=255)
  gctgctCCAGTCATTTTCGCCCTTAATGCGCGTCAg                                  >  1:6506802/1‑36 (MQ=255)
     gctCCAGTCATTTTCGCCCTTAATGCGCGTGAtcgc                               <  1:3450983/36‑1 (MQ=37)
      ctCCAGTCATTTTCGCCCTTAATGCGCGTCAGcgcg                              <  1:6279002/36‑1 (MQ=255)
      ctCCAGTCATTTTCGCCCTTAATGCGCGTCAGcgcg                              <  1:2166563/36‑1 (MQ=255)
       tCCAGTCATTTTCGCCCTTAATGCGCGTGATCGCGt                             >  1:7875757/1‑36 (MQ=37)
       tCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGt                             >  1:4933195/1‑36 (MQ=255)
       tCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGt                             >  1:4431396/1‑36 (MQ=255)
       tCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGt                             >  1:5320756/1‑36 (MQ=255)
        ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                            <  1:6270163/36‑1 (MQ=255)
        ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                            <  1:395748/36‑1 (MQ=255)
        ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                            >  1:7146532/1‑36 (MQ=255)
        ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                            >  1:491102/1‑36 (MQ=255)
        ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                            <  1:7893576/36‑1 (MQ=255)
         cAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGa                           >  1:1905940/1‑36 (MQ=255)
             cATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGa                       <  1:4012540/36‑1 (MQ=255)
             cATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGa                       >  1:5239554/1‑36 (MQ=255)
              aTTTTCGCCCTTAATGCGCGTCAGCTCGTGATCGac                      >  1:7915342/1‑36 (MQ=255)
              aTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGac                      >  1:7324201/1‑36 (MQ=255)
              aTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGac                      >  1:6249620/1‑36 (MQ=255)
               ttttCGCCCTTAATGCGCGTCAGCGCGTGATCGaca                     >  1:3831337/1‑36 (MQ=255)
               ttttCGCCCTTAATGCGCGTCAGCGCGTGATCGaca                     >  1:7643453/1‑36 (MQ=255)
                tttCGCCCTTAATGCGCGTCAGCGCGTGATCGacac                    >  1:8044522/1‑36 (MQ=255)
                 ttCGCCCTTAATGCGCGTGATCGCGTGATCGACACg                   <  1:2505139/36‑1 (MQ=37)
                 ttCGCCCTTAATGCGCGTGATCGCGTGATCGACACg                   <  1:2244275/36‑1 (MQ=37)
                 ttCGCCCTTAATGCGCGTCAGCGCGTGATCGACACg                   >  1:4372829/1‑36 (MQ=255)
                 ttCGCCCTTAATGCGCGTCAGCGCGTGATCGACACg                   <  1:5594812/36‑1 (MQ=255)
                   cGCCCTTAATGCGCGTCAGCGCGTGATCGACACGtt                 >  1:4117719/1‑36 (MQ=255)
                   cGCCCTTAATGCGCGTCAGCGCGTGATCGACACGtt                 <  1:5388656/36‑1 (MQ=255)
                    gCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                <  1:8435857/36‑1 (MQ=255)
                    gCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                <  1:2904798/36‑1 (MQ=255)
                     cccTTAATGCGCGTCGGCGCGTGATCGACACGTTcc               <  1:2810603/36‑1 (MQ=255)
                     cccTTAATGCGCGTCAGCGCGTGATCGACACGTTcc               <  1:1351931/36‑1 (MQ=255)
                     cccTTAATGCGCGTCAGCGCGTGATCGACACGTTcc               <  1:8268271/36‑1 (MQ=255)
                     cccTTAATGCGCGTCAGCGCGTGATCGACACGTTcc               <  1:3754129/36‑1 (MQ=255)
                     cccTTAAGGCGCGTCAGCGCGTGATCGACACGTTcc               <  1:2730420/36‑1 (MQ=255)
                      ccTTAATGCGCGTGATCGCGTGATCGACACGTTCCt              >  1:6975227/1‑36 (MQ=37)
                      ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt              <  1:7650540/36‑1 (MQ=255)
                      ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt              <  1:6895331/36‑1 (MQ=255)
                      ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt              >  1:160457/1‑36 (MQ=255)
                      ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt              >  1:4300546/1‑36 (MQ=255)
                      ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt              <  1:4577408/36‑1 (MQ=255)
                       cTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTc             <  1:1359489/36‑1 (MQ=255)
                       cTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTc             >  1:6969242/1‑36 (MQ=255)
                         tAATGCGCGTGATCGCGTGATCGACACGTTCCTCaa           <  1:4877620/36‑1 (MQ=37)
                         tAATGCGCGTCAGCGCGTGATCGACACGTTCCTCaa           >  1:2036390/1‑36 (MQ=255)
                          anTGCGCGTCAGCGCGTGATCGACACGTTCCTCaaa          <  1:3941819/36‑1 (MQ=255)
                          aaTGCGCGTCAGCGCGTGATCGACACGTTCCTCaaa          >  1:7481378/1‑36 (MQ=255)
                          aaTGCGCGTCAGCGCGTGATCGACACGTTCCTCaaa          >  1:8176637/1‑36 (MQ=255)
                           aTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAAc         >  1:4493378/1‑36 (MQ=255)
                           aTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAAc         <  1:8107985/36‑1 (MQ=255)
                            tGCGCGTCAGCGCGTGATCGCCACGTTCCTCAAACg        >  1:4367467/1‑36 (MQ=255)
                            tGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACg        >  1:72727/1‑36 (MQ=255)
                             gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg       >  1:7105084/1‑36 (MQ=255)
                             gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg       <  1:3773324/36‑1 (MQ=255)
                             gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg       <  1:6268156/36‑1 (MQ=255)
                             gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg       <  1:5425726/36‑1 (MQ=255)
                              cgcgTCAGCGCGTGATCGACACGTTCCTCAAACGGt      <  1:8066221/36‑1 (MQ=255)
                              cgcgTCAGCGCGTGATCGACACGTTCCTCAAACGGt      <  1:5031593/36‑1 (MQ=255)
                               gcgTCAGCGCGTGATCGACACGTTCCTCAAACGGtt     <  1:666236/36‑1 (MQ=255)
                                 gccAGCGCGGGATCGACACGTTCCTCAAACGGTTTc   <  1:5291999/34‑1 (MQ=37)
                                  tCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCa  <  1:52910/36‑1 (MQ=255)
                                   |                                  
TCGCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCA  >  REL606/4561411‑4561480

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG