Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,561,448 G→T 7.6% L1097M (CTG→ATG)  hsdR ← endonuclease R

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,561,4480G→T 7.6% 5.5 0.0 ‑0.4663289hsdRendonuclease R
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (28/33):  total (30/36)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.70e-01

GCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCG  >  REL606/4561413‑4561483
                                   |                                   
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:5320756/1‑36 (MQ=255)
     tCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGt                                >  1:4431396/1‑36 (MQ=255)
      ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                               >  1:491102/1‑36 (MQ=255)
      ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                               <  1:395748/36‑1 (MQ=255)
      ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                               <  1:6270163/36‑1 (MQ=255)
      ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                               >  1:7146532/1‑36 (MQ=255)
      ccAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTg                               <  1:7893576/36‑1 (MQ=255)
       cAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGa                              >  1:1905940/1‑36 (MQ=255)
           cATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGa                          >  1:5239554/1‑36 (MQ=255)
           cATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGa                          <  1:4012540/36‑1 (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:7643453/1‑36 (MQ=255)
             ttttCGCCCTTAATGCGCGTCAGCGCGTGATCGaca                        >  1:3831337/1‑36 (MQ=255)
              tttCGCCCTTAATGCGCGTCAGCGCGTGATCGacac                       >  1:8044522/1‑36 (MQ=255)
               ttCGCCCTTAATGCGCGTGATCGCGTGATCGACACg                      <  1:2244275/36‑1 (MQ=37)
               ttCGCCCTTAATGCGCGTGATCGCGTGATCGACACg                      <  1:2505139/36‑1 (MQ=37)
               ttCGCCCTTAATGCGCGTCAGCGCGTGATCGACACg                      <  1:5594812/36‑1 (MQ=255)
               ttCGCCCTTAATGCGCGTCAGCGCGTGATCGACACg                      >  1:4372829/1‑36 (MQ=255)
                 cGCCCTTAATGCGCGTCAGCGCGTGATCGACACGtt                    >  1:4117719/1‑36 (MQ=255)
                 cGCCCTTAATGCGCGTCAGCGCGTGATCGACACGtt                    <  1:5388656/36‑1 (MQ=255)
                  gCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                   <  1:2904798/36‑1 (MQ=255)
                  gCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTc                   <  1:8435857/36‑1 (MQ=255)
                   cccTTAATGCGCGTCGGCGCGTGATCGACACGTTcc                  <  1:2810603/36‑1 (MQ=255)
                   cccTTAATGCGCGTCAGCGCGTGATCGACACGTTcc                  <  1:8268271/36‑1 (MQ=255)
                   cccTTAATGCGCGTCAGCGCGTGATCGACACGTTcc                  <  1:1351931/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:160457/1‑36 (MQ=255)
                    ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt                 <  1:7650540/36‑1 (MQ=255)
                    ccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCt                 <  1:6895331/36‑1 (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:8176637/1‑36 (MQ=255)
                        aaTGCGCGTCAGCGCGTGATCGACACGTTCCTCaaa             >  1:7481378/1‑36 (MQ=255)
                         aTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAAc            <  1:8107985/36‑1 (MQ=255)
                         aTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAAc            >  1:4493378/1‑36 (MQ=255)
                          tGCGCGTCAGCGCGTGATCGCCACGTTCCTCAAACg           >  1:4367467/1‑36 (MQ=255)
                          tGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACg           >  1:72727/1‑36 (MQ=255)
                           gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg          <  1:6268156/36‑1 (MQ=255)
                           gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg          <  1:5425726/36‑1 (MQ=255)
                           gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg          >  1:7105084/1‑36 (MQ=255)
                           gcgcGTCAGCGCGTGATCGACACGTTCCTCAAACgg          <  1:3773324/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)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   >  1:181281/1‑36 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   >  1:1712038/1‑36 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   <  1:1694059/36‑1 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   <  1:3745415/36‑1 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   <  1:6581994/36‑1 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   <  1:531029/36‑1 (MQ=255)
                                  aGCGCGTGATCGACACGTTCCTCAAACGGTTTCAgc   >  1:6647575/1‑36 (MQ=255)
                                   gcgcGTGATCGACACGTTCCTCAAACGGTTTCAgcg  <  1:5159448/36‑1 (MQ=255)
                                   gcgcGTCATCAACACGTTCCTGAAACGGTTTCAgcg  <  1:1940082/36‑1 (MQ=37)
                                   |                                   
GCTGCTCCAGTCATTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCG  >  REL606/4561413‑4561483

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