Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,201,911 C→T 81.5% A217T (GCG→ACG)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,201,9110C→T81.5% 35.1 ‑0.9 0.081649iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (35/31):  ref base (6/9):  total (41/40)
Fisher's exact test strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGT  >  REL606/4201878‑4201942
                                 |                               
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                               >  1:5104468/1‑36 (MQ=255)
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                               >  1:1753948/1‑36 (MQ=255)
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCAt                               <  1:1330707/36‑1 (MQ=255)
 aaTACACGCTGCAAGGCAACGTAGCCCCAGCGTATg                              <  1:3267921/36‑1 (MQ=255)
  aTACACGCTGCAAGGCCACGTAGCCCCAGCGTATGt                             <  1:4752902/36‑1 (MQ=37)
   tACACGCTGCAAGGCAACGTAGCCCCAGCGTATGtt                            >  1:6102819/1‑36 (MQ=255)
   tACACGCTGCAAGGCAACGTAGCCCCAGCGTATGtt                            <  1:1147803/36‑1 (MQ=255)
   tACACGCTGCAAGGCAACGTAGCCCCAGCGTATGtt                            <  1:5694930/36‑1 (MQ=255)
   tACACGCTGCAAGGCAACGTAGCCCCAGCGCATGtt                            >  1:1118712/1‑36 (MQ=255)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                           >  1:6551366/1‑36 (MQ=255)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                           <  1:3244109/36‑1 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                         >  1:3376813/1‑36 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                         >  1:8086659/1‑36 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                         >  1:4945401/1‑36 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCt                         <  1:835024/36‑1 (MQ=255)
       cGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTc                        <  1:4997746/36‑1 (MQ=255)
       cGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTc                        >  1:6365403/1‑36 (MQ=255)
       cGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTc                        <  1:4705743/36‑1 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                       <  1:2409863/36‑1 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                       >  1:2545078/1‑36 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                       >  1:3168752/1‑36 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCa                       >  1:3075251/1‑36 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCa                       >  1:6581017/1‑36 (MQ=255)
         cTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCAt                      >  1:3526010/1‑36 (MQ=255)
         cTGCAAGGCAACGCAGCCCCAGCGTATGTTCCTCAt                      >  1:2167058/1‑36 (MQ=37)
          tGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATc                     <  1:8021233/36‑1 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                    >  1:518742/1‑36 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                    >  1:4218514/1‑36 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                    >  1:8357080/1‑36 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                    <  1:2500543/36‑1 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                    <  1:142956/36‑1 (MQ=255)
            cgaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGt                   <  1:6992210/34‑1 (MQ=38)
            cAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGt                   <  1:754767/36‑1 (MQ=255)
            cAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGt                   >  1:6272934/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTTCTCATCGta                  >  1:105096/1‑35 (MQ=37)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                  >  1:3458759/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                  <  1:7422441/36‑1 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                  >  1:1148709/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                  <  1:7839541/36‑1 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                  <  1:8007405/36‑1 (MQ=255)
             aaGGCAACGGAGACCCAGCGTATGTTCCTCATCGTc                  >  1:6549895/1‑36 (MQ=37)
              aGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                 <  1:542060/36‑1 (MQ=255)
              aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCa                 <  1:1281563/36‑1 (MQ=255)
               ggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                >  1:5479574/1‑36 (MQ=255)
               ctcAACGTAGCCCCAGCGCATGTTCCTCATCGTCaa                <  1:241225/34‑1 (MQ=255)
                gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa               <  1:7047585/36‑1 (MQ=255)
                gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa               >  1:4429415/1‑36 (MQ=255)
                gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa               >  1:7362273/1‑36 (MQ=255)
                gCAACGTAGCCCCAGCGCATGTTCCTCATCGTCaaa               <  1:5307385/36‑1 (MQ=255)
                 cAACGTAGCCCCAGCGTATGTTCCTCATCGTCAAAt              <  1:6564000/36‑1 (MQ=255)
                 cAACGTAGCCCCAGCGTATGTTCCTCATCGTCAAAt              <  1:7822635/36‑1 (MQ=255)
                 cAACGTAGCCCCAGCGTATGTTCCTCATCGTCAAAt              <  1:3631501/36‑1 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg             >  1:4489198/1‑36 (MQ=255)
                  aaCCTAGCCCCAGCGTATGTTCCTCATCGTCAAATg             <  1:654721/36‑1 (MQ=37)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa            <  1:1545347/36‑1 (MQ=255)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa            <  1:469799/36‑1 (MQ=255)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTAAAATGa            >  1:8168280/1‑36 (MQ=37)
                    cGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGaa           >  1:1043286/1‑36 (MQ=255)
                    cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGaa           >  1:899001/1‑36 (MQ=255)
                    cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGaa           >  1:5647539/1‑36 (MQ=255)
                     gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt          >  1:8355718/1‑36 (MQ=255)
                     gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt          <  1:6462229/36‑1 (MQ=255)
                      tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa         >  1:2510742/1‑36 (MQ=255)
                      tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa         <  1:5322704/36‑1 (MQ=255)
                       aGCCCCAGCGTATGTTCCTCATCGTCAAATGAATaa        <  1:7694224/36‑1 (MQ=255)
                       aGCCCCAGCGTATGTTCCTCATCGTCAAATGAATaa        <  1:1283086/36‑1 (MQ=255)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc       <  1:436208/36‑1 (MQ=255)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc       <  1:8420838/36‑1 (MQ=255)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc       <  1:5717970/36‑1 (MQ=255)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc       >  1:4659558/1‑36 (MQ=255)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc      >  1:5842272/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc     >  1:404500/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc     >  1:7475813/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc     <  1:1598129/36‑1 (MQ=255)
                           ccAGCGTATGTTCCTCATCGTCAAATGAATAAcccc    >  1:3190825/1‑36 (MQ=255)
                           ccAGCGTATGTTCCTCATCGTCAAATGAATAAcccc    >  1:2970908/1‑36 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg   >  1:304609/1‑36 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg   <  1:3353032/36‑1 (MQ=255)
                            cAGCGCATGTTCCTCATCGTCAAATGAATAACCCCg   <  1:7972921/36‑1 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt  >  1:4441772/1‑36 (MQ=255)
                             aGCGCATGTTCCTCATCGTCAAATGAATAACCCCGt  <  1:3915052/36‑1 (MQ=255)
                                 |                               
AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGT  >  REL606/4201878‑4201942

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG