Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,201,911 C→T 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% 180.3 59649iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (29/28):  ref base (0/1):  total (29/30)

AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943
                                 |                                
aaaTACGCGCGGCGAGGCAACGTAGCCCCAGCGTAt                                <  1:1178323/36‑1 (MQ=19)
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                >  2:1891493/1‑36 (MQ=255)
aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                <  1:2252761/36‑1 (MQ=255)
aaaAACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                <  1:1738620/36‑1 (MQ=37)
 aaTACACGCTGCAAGGCAACGTAGCCCCCGCGTATg                               >  2:225789/1‑36 (MQ=38)
 aaTACACGCTGCAAGGCAACGTAGCCCCAGCGTATg                               >  2:2695111/1‑36 (MQ=255)
  aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              >  2:440666/1‑36 (MQ=255)
  aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              >  1:2156437/1‑36 (MQ=255)
   tACACGCTGCAAGGCAACGTAGCCCCAGCGTATGtt                             <  2:388340/36‑1 (MQ=255)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGtta                            >  1:539173/1‑35 (MQ=255)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                            >  1:3308928/1‑36 (MQ=255)
    acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                            <  1:2285619/36‑1 (MQ=255)
     cacGCTGCAAGGCAACGTAGCTCCAGCGTATGTTcc                           <  1:2838817/36‑1 (MQ=38)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           >  1:322309/1‑36 (MQ=255)
     cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  2:1194294/36‑1 (MQ=255)
      acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                          >  1:2833591/1‑36 (MQ=255)
       cGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTc                         <  1:3124998/36‑1 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                        >  2:1333110/1‑36 (MQ=255)
        gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                        >  2:1371855/1‑36 (MQ=255)
         cTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCAt                       >  1:863483/1‑36 (MQ=255)
          tGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATc                      <  1:2353411/36‑1 (MQ=255)
           gaaaGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     <  1:118627/34‑1 (MQ=38)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     >  1:813299/1‑36 (MQ=255)
           gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     >  2:3486458/1‑36 (MQ=255)
            cAAGGCAGCGTAGCCCCAGCGTATGTTCCTCATCGt                    <  2:1640978/36‑1 (MQ=38)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCCTCGTc                   >  1:3460714/1‑36 (MQ=37)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  2:3242586/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  2:331397/1‑36 (MQ=255)
             aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   <  1:837876/36‑1 (MQ=255)
              tggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  >  1:2198307/2‑36 (MQ=38)
              atgCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  <  2:773282/34‑1 (MQ=38)
              aGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  <  2:622822/36‑1 (MQ=255)
               ggCTACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                 <  2:628270/36‑1 (MQ=37)
               ggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                 <  1:2205778/36‑1 (MQ=255)
               gcgaaCGTAGCCCCAGCGCATGTTCCTCATCGTCaa                 <  2:3589437/33‑1 (MQ=255)
                gaaaCGAAGCCCCAGCGAATGTTCCTCATCGTCaaa                <  2:604152/34‑1 (MQ=37)
                gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa                >  1:1837897/1‑36 (MQ=255)
                 cAACGTAGCCCCAGCGTATGTTCCTCATCGTCAAAt               <  1:3220794/36‑1 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCATATg              >  2:689957/1‑36 (MQ=37)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              >  2:2755616/1‑36 (MQ=255)
                  aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              <  1:134193/36‑1 (MQ=255)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa             <  2:2289631/36‑1 (MQ=255)
                   aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa             <  2:2080645/36‑1 (MQ=255)
                     gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt           >  2:3636888/1‑36 (MQ=255)
                     gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt           >  2:2715654/1‑36 (MQ=255)
                      tAGTCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  2:1677817/36‑1 (MQ=37)
                      tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  2:2231972/36‑1 (MQ=255)
                      tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  2:829144/36‑1 (MQ=255)
                       aGCCCCAGCGTATGTTCCTCATCGTCAAATGAATaa         >  1:1322666/1‑36 (MQ=255)
                       aGCCCCAGCGTATGTTCCTCATCGTAAAAGGACTaa         >  2:2535832/1‑36 (MQ=19)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc        >  2:904714/1‑36 (MQ=255)
                        gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc        <  1:393560/36‑1 (MQ=255)
                         ccctAGCGTATGTTCCTCATCGTCAAATGAATAAcc       <  2:1954486/36‑1 (MQ=37)
                         ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       >  1:3286873/1‑36 (MQ=255)
                          cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc      <  1:2957970/36‑1 (MQ=255)
                           ccAGCGTATGTTCCTCATCGTCAAATGAATAAcccc     <  2:1967980/36‑1 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    >  2:3614840/1‑36 (MQ=255)
                            cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    >  2:1246480/1‑36 (MQ=255)
                             aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt   <  1:1333136/36‑1 (MQ=255)
                              gcgtATGTTCCTCATCGTCAAATGAATAACCCCGtt  >  2:1511297/1‑36 (MQ=255)
                                 |                                
AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG