Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,560,632 T→C 30.9% Y131C (TAC→TGC)  hsdM ← DNA methylase M

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,560,6320T→C30.9% 67.1 0.0 ‑0.781392hsdMDNA methylase M
Reads supporting (aligned to +/- strand):  new base (15/10):  ref base (32/24):  total (47/34)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.18e-01

TGGTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGTGACTTACCGTGC  >  REL606/4560597‑4560667
                                   |                                   
ttgTTCCGTTGGCGTTATTCTGCAACAGCCCTTCGt                                     <  1:6309433/34‑1 (MQ=37)
tGGTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGt                                     <  1:7489134/36‑1 (MQ=255)
tGGTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGt                                     >  1:5866569/1‑36 (MQ=255)
 ggTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTa                                    >  1:7728700/1‑36 (MQ=255)
 ggTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTa                                    >  1:3593119/1‑36 (MQ=255)
  gTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTAc                                   >  1:981736/1‑36 (MQ=255)
  gTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTAc                                   <  1:1145060/36‑1 (MQ=255)
  gTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTAc                                   >  1:7806300/1‑36 (MQ=255)
  gTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGCAc                                   >  1:4603328/1‑36 (MQ=255)
   tttCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACa                                  >  1:640712/1‑36 (MQ=255)
   tttCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACa                                  <  1:5115618/36‑1 (MQ=255)
   tttCGTTGGCGTTTTTCTGCAACAGCCCTTCGCACa                                  <  1:416474/36‑1 (MQ=255)
    ttCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACat                                 >  1:6510536/1‑36 (MQ=255)
    ttCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACat                                 <  1:2853600/36‑1 (MQ=255)
    ttCGTTGGCGTTTTTCTGCAACAGCCCTTCGCACat                                 <  1:1850091/36‑1 (MQ=255)
     tCGTTGGCGTTTTTCTGCAACAGCCCTTCGCACata                                <  1:6598894/36‑1 (MQ=255)
      cGTTGGCGTTTTTCTGCAACAGCCCTTCGTACatat                               <  1:1681444/36‑1 (MQ=255)
      cGTTGGCGTTTTTCTGCAACAGCCCTTCGCACatat                               >  1:4228358/1‑36 (MQ=255)
       gTTGGCGTTTTTCTGCAACAGCCCTTCGTACATATc                              <  1:5254009/36‑1 (MQ=255)
       gTTGGCGTTTTTCTGCAACAGCCCTTCGCACATATc                              <  1:7508329/36‑1 (MQ=255)
        ttGGCGTTTTTCTGCAACAGCCCTTCGTACATATCg                             <  1:3155416/36‑1 (MQ=255)
         tGGCGTTTTTCTGCAACAGCCCTTCGTACATATCGc                            >  1:273226/1‑36 (MQ=255)
         tGGCGTTTTTCTGCAACAGCCCTTCGCACATATCGc                            >  1:6505956/1‑36 (MQ=255)
          ggCGTTTTTCTGCAACAGCCCTTCGTACATATCGcc                           <  1:6077641/36‑1 (MQ=255)
          ggCGTTTTTCTGCAACAGCCCTTCGCACATATCGcc                           <  1:4348666/36‑1 (MQ=255)
          ggCGTTTTTCTGCAACAGCCCTTCGCACATATCGcc                           >  1:5589342/1‑36 (MQ=255)
           gCGTTTTTCTGCAACAGCCCTTCGTACATATCGCCg                          <  1:3078996/36‑1 (MQ=255)
           gCGGTTTTCTGCAACAGCCCTTCTTTCATATCGCCg                          >  1:4104693/1‑36 (MQ=25)
            cGTTTTTCTGCAACAGCCCTTCGCACATATCGCCGa                         >  1:1405095/1‑36 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGTACATATCGCCGaa                        >  1:3446173/1‑36 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGTACATATCGCCGaa                        >  1:376321/1‑36 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGTACATATCGCCGaa                        <  1:2998756/36‑1 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGTACATATCGCCGaa                        >  1:7420179/1‑36 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGTACATATCGCCGaa                        <  1:63707/36‑1 (MQ=255)
             gTTTTTCTGCAACAGCCCTTCGCACATATCGCCGaa                        >  1:3676726/1‑36 (MQ=255)
               ttttCTGCAACAGCCCTTCGTACATATCGCCGAAGt                      >  1:3062594/1‑36 (MQ=255)
                tttCTGCAACAGCCCTTCGTACATATCGCCGCAGTc                     >  1:6906866/1‑36 (MQ=255)
                tttCTGCAACAGCCCTTCGTACATATCGCCGAAGTc                     >  1:7679785/1‑36 (MQ=255)
                tttCTGCAACAGCCCTTCGCACATATCGCCGAAGTc                     >  1:3067341/1‑36 (MQ=255)
                 ttCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                    >  1:5937408/1‑36 (MQ=255)
                 ttCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                    >  1:863951/1‑36 (MQ=255)
                 ttCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                    <  1:3708155/36‑1 (MQ=255)
                 ttCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                    <  1:1810596/36‑1 (MQ=255)
                  tCTGCAACAGCCCTTCGTACATATCGCCGAAGTCAt                   >  1:5116499/1‑36 (MQ=255)
                  tCTGCAACAGCCCTTCGTACATATCGCCGAAGTCAt                   >  1:3698292/1‑36 (MQ=255)
                  tCTGCAACAGCCCTTCGCACATATCGCCGAAGTCAt                   <  1:7108804/36‑1 (MQ=255)
                   cTGCAACAGCCCTTCGTACATATCGCCGAAGTCATc                  <  1:6320537/36‑1 (MQ=255)
                   cTGCAACAGCCCTTCGTACATATCGCCGAAGTCATc                  <  1:4931692/36‑1 (MQ=255)
                   cTGCAACAGCCCTTCGCACATATCGCCGAAGTCATc                  >  1:4559179/1‑36 (MQ=255)
                   cTGCAACAGCCCTTCGCACATATCGCCGAAGTCATc                  >  1:4889712/1‑36 (MQ=255)
                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATcg                 >  1:3705254/1‑36 (MQ=255)
                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATcg                 <  1:2053357/36‑1 (MQ=255)
                    tGCAACAGCCCTTCGCACATATCGCCGAAGTCATcg                 <  1:5795564/36‑1 (MQ=255)
                     gCAACAGCCCTTCGTACATATCGCCGAAGTCATcgc                >  1:1214074/1‑36 (MQ=255)
                     gCAACAGCCCTTCGCACATATCGCCGAAGTCATcgc                >  1:6422723/1‑36 (MQ=255)
                       aaCAGCCCTTCGTACATAGCGCCGAAGTCATagcgg              >  1:2146432/1‑35 (MQ=37)
                       aaCAGCCCTTCGCACATATCGCCGAAGTCATCGCGt              <  1:5352568/36‑1 (MQ=255)
                       aaCAGCCCTTCGCACATATCGCCGAAGTCATCGCGt              <  1:3090493/36‑1 (MQ=255)
                        aCAGCCCTTCGCACATATCGCCGAAGTCATCGCGTg             >  1:901568/1‑36 (MQ=255)
                          aGCCCTTCGTACATATCGCCGAAGTCATCGCGTGAc           >  1:1899610/1‑36 (MQ=255)
                          aGCCCTTCGTACATATCGCCGAAGTCATCGCGTGAc           >  1:6542065/1‑36 (MQ=255)
                           gCCCTTCGTACATATCGCCGAAGTCATCGCGTGACt          <  1:1591729/36‑1 (MQ=255)
                           gCCCTTCGTACATATCGCCGAAGTCATCGCGTGACt          <  1:1291877/36‑1 (MQ=255)
                            cccTTCGTACATATCGCCGAAGTCATCGCGTGACtt         >  1:1967027/1‑36 (MQ=255)
                            cccTTCGTACATATCGCCGAAGTCATCGCGTGACtt         >  1:7922059/1‑36 (MQ=255)
                             ccTTCGTACATATCGCCGAAGTCATCGCGTGACTTa        >  1:1036570/1‑36 (MQ=255)
                             ccTTCGTACATATCGCCGAAGTCATCGCGTGACTTa        >  1:5948359/1‑36 (MQ=255)
                             ccTTCGTACATATCGCCGAAGTCATCGCGTGACTTa        <  1:5199158/36‑1 (MQ=255)
                             ccTTCGTACATATCGCCGAAGTCATCGCGTGACTTa        >  1:505061/1‑36 (MQ=255)
                             ccTTCGTACATAACACCGAAGTCATCGCGTCACTTa        >  1:1347106/1‑36 (MQ=37)
                              cTTCGTACATATCGCCGAAGTCATCGCGTGACTTAc       <  1:6581470/36‑1 (MQ=255)
                              cTTCGTACATATCGCCGAAGTCATCGCGTGACTTAc       <  1:928351/36‑1 (MQ=255)
                              cTTCGTACATATCGCCGAAGTCATCGCGTGACTTAc       >  1:1269322/1‑36 (MQ=255)
                              cTTCGCACATATCGCCGAAGTCATCGCGTGACTTAc       >  1:1813294/1‑36 (MQ=255)
                               ttCGTACATATCGCCGAAGTCATCGCGTGACTTAcc      <  1:7283581/36‑1 (MQ=255)
                               ttCGTACATATCGCCGAAGTCATCGCGTGACTTAcc      >  1:2113401/1‑36 (MQ=255)
                               ttCGCACATATCGCCGAAGTCATCGCGTGACTTAcc      >  1:7775665/1‑36 (MQ=255)
                                tCGCACATATCGCCGAAGTCATCGCGTGACTTACCg     >  1:2576582/1‑36 (MQ=255)
                                tCGCACATATCGCCGAAGTCATCGCGTGACTTACCg     >  1:7531826/1‑36 (MQ=255)
                                tCGCACATATCGCCGAAGTCATCGCGTGACTTACCg     <  1:3157810/36‑1 (MQ=255)
                                 cGTACATATCGCCGAAGTCATCGCGTGACTTACCGt    <  1:2703942/36‑1 (MQ=255)
                                 cGTACATATCGCCGAAGTCATCGCGTGACTTACCGt    <  1:4509183/36‑1 (MQ=255)
                                 cGTACATATCGCCGAAGTCATCGCGTGACTTACCGt    <  1:6627881/36‑1 (MQ=255)
                                  gTACATATCGCCGAAGTCATCGCGTGACTTACCGTg   >  1:5429090/1‑36 (MQ=255)
                                   tACATATCGCCGAAGTCATCGCGTGACTTACCGTgc  <  1:25930/36‑1 (MQ=255)
                                   tACATATCGCCGAAGTCATCGCGTGACTTACCGTgc  <  1:1782913/36‑1 (MQ=255)
                                   |                                   
TGGTTTCGTTGGCGTTTTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGTGACTTACCGTGC  >  REL606/4560597‑4560667

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG