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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,201,9110C→T100.0% 176.2 49649iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (22/27):  ref base (0/0):  total (22/27)

AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943
                                 |                                
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                <  1:4379217/36‑1 (MQ=255)
 aaTACACGCTGCAAGGCAACGTAGCCCCAGCGTATg                               >  1:382865/1‑36 (MQ=255)
 aaTACACGCTGCAAGGCAACGTAGCCCCAGCGTATg                               <  1:1911637/36‑1 (MQ=255)
  aTACACGCTGCAAGGCAATGTAGCCCCAGCGTATGt                              >  1:133927/1‑36 (MQ=38)
  aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              >  1:2585635/1‑36 (MQ=255)
  aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              <  1:1103633/36‑1 (MQ=255)
   tACACGCTGCAAGGCAACGAAACCCCAGCGTATGtt                             >  1:5840535/1‑36 (MQ=25)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                            <  1:1145912/36‑1 (MQ=255)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  1:4995532/36‑1 (MQ=255)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  1:5646090/36‑1 (MQ=255)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  1:2375952/36‑1 (MQ=255)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  1:2559549/36‑1 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                          >  1:616175/1‑36 (MQ=255)
       cGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTc                         <  1:3724855/36‑1 (MQ=255)
         cTGCAAGGCAACGTAACCCCAGCGTATGTTCCTCAt                       <  1:1699156/36‑1 (MQ=37)
          tGCAAGGCAACTTAGCCCCAGCGTATGTTCCTCATc                      <  1:1191397/36‑1 (MQ=37)
          tGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATc                      <  1:1613500/36‑1 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     >  1:323304/1‑36 (MQ=255)
            cAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCGt                    <  1:4642556/36‑1 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   <  1:677689/36‑1 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  1:1577418/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   <  1:896167/36‑1 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  1:1347974/1‑36 (MQ=255)
               ggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                 <  1:3953546/36‑1 (MQ=255)
                gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa                >  1:5268030/1‑36 (MQ=255)
                 cAACGTAGCCCCATCGTATGTTCCTCATCGTCAAAt               <  1:5339061/36‑1 (MQ=37)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              >  1:3047263/1‑36 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              >  1:5275500/1‑36 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              >  1:4772807/1‑36 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              <  1:944477/36‑1 (MQ=255)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa             <  1:4270392/36‑1 (MQ=255)
                    cGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGaa            <  1:4327335/36‑1 (MQ=255)
                    cGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGaa            >  1:3781519/1‑36 (MQ=255)
                     gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt           <  1:1698220/36‑1 (MQ=255)
                       aGCCCCAGCGTATGTTCCTCATCGTCAAATGAATaa         >  1:5573482/1‑36 (MQ=255)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       >  1:997253/1‑36 (MQ=255)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       <  1:62642/36‑1 (MQ=255)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       >  1:3079517/1‑36 (MQ=255)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       <  1:4276760/36‑1 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc      >  1:3344735/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc      >  1:5408993/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCAGCGTCAAATGAATAAccc      >  1:3649986/1‑36 (MQ=37)
                           ccAGCGTATGTTCCTCATCGTCAAATGAATAAcccc     <  1:4995729/36‑1 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    >  1:3595349/1‑36 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    <  1:4800774/36‑1 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt   <  1:3506069/36‑1 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt   <  1:5132375/36‑1 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt   >  1:1507026/1‑36 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATAAATAACACCg  >  1:1107219/1‑35 (MQ=21)
                              gcgtATGTTCCTCATCGTCAAATGAATAACCCCGtt  <  1:32421/36‑1 (MQ=255)
                              gcgtATGTTCCTCATCGTCAAATGAATAACCCCGtt  >  1:1868032/1‑36 (MQ=255)
                                 |                                
AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG