Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,201,934 T→G Y209S (TAT→TCT)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,201,9340TG95.5% 262.5 / 2.1 66Y209S (TAT→TCT) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (37/26);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

CACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGG  >  REL606/4201863‑4202006
                                                                       |                                                                        
cACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAAcc                                                                      >  2:757536/1‑76 (MQ=255)
  cGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAAcccc                                                                    <  1:113680/76‑1 (MQ=255)
   ggTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCg                                                                   <  1:1962241/76‑1 (MQ=255)
      gCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGttt                                                                >  1:411504/1‑76 (MQ=255)
      gCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGttt                                                                >  1:1583699/1‑76 (MQ=255)
      gCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGttt                                                                <  2:1964516/76‑1 (MQ=255)
       cTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTg                                                               >  2:327796/1‑76 (MQ=255)
        tCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTgc                                                              >  1:1397856/1‑76 (MQ=255)
        tCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTgc                                                              >  1:978611/1‑76 (MQ=255)
         cATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTgcg                                                             >  2:461282/1‑76 (MQ=255)
         cATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTgcg                                                             >  2:1691960/1‑76 (MQ=255)
          aTCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTgcgc                                                            >  2:1693232/1‑76 (MQ=255)
           tCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGcgcg                                                           >  2:1676058/1‑76 (MQ=255)
             gAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGtt                                                         >  1:434590/1‑76 (MQ=255)
               aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTg                                                       >  1:13158/1‑76 (MQ=255)
               aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTg                                                       <  2:389550/76‑1 (MQ=255)
                 aTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTggg                                                     >  2:434526/1‑76 (MQ=255)
                   tcgcGCTGCAGGGCAGCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCg                                                   <  2:1521128/73‑1 (MQ=255)
                   acacGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCg                                                   <  1:231175/76‑1 (MQ=255)
                       gCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAt                                               >  1:620724/1‑76 (MQ=255)
                        cTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtc                                              >  1:718495/1‑76 (MQ=255)
                         tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtct                                             >  1:339016/1‑76 (MQ=255)
                         tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtct                                             >  2:742289/1‑76 (MQ=255)
                          gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGCGAtctt                                            >  1:532650/1‑76 (MQ=255)
                          gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtctt                                            >  2:192696/1‑76 (MQ=255)
                          gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtctt                                            >  1:665748/1‑76 (MQ=255)
                           cAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtcttc                                           >  2:1519073/1‑76 (MQ=255)
                            aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                                          <  2:650442/76‑1 (MQ=255)
                             aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                                         <  2:852487/76‑1 (MQ=255)
                             aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                                         >  2:194708/1‑76 (MQ=255)
                             aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                                         <  2:1834499/76‑1 (MQ=255)
                                cAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTa                                      <  2:49825/76‑1 (MQ=255)
                                 aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa                                     <  1:1868445/76‑1 (MQ=255)
                                 aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa                                     >  2:413530/1‑76 (MQ=255)
                                   cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAAt                                   >  1:813960/1‑76 (MQ=255)
                                    gTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATg                                  >  2:410101/1‑76 (MQ=255)
                                     tAGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc                                 >  2:895986/1‑76 (MQ=255)
                                      aGCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGca                                >  1:1947902/1‑76 (MQ=255)
                                       gCCCCAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac                               <  1:270423/76‑1 (MQ=255)
                                        ccccAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATGcaca                              >  2:1536327/1‑76 (MQ=255)
                                          ccAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAgg                            <  2:1814157/76‑1 (MQ=255)
                                          ccAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAgg                            <  1:418428/76‑1 (MQ=255)
                                           cAGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGga                           >  1:75940/1‑76 (MQ=255)
                                            aGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGgag                          <  2:1505642/76‑1 (MQ=255)
                                            aGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGgag                          >  2:858706/1‑76 (MQ=255)
                                            aGCGCATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGgag                          <  1:1412143/76‑1 (MQ=255)
                                             gcgcATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGgaga                         <  2:560641/76‑1 (MQ=255)
                                               gcATGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGaca                       <  1:1829807/76‑1 (MQ=255)
                                                 aTGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACAcc                     <  2:270171/76‑1 (MQ=255)
                                                 aTGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACAcc                     >  1:1108081/1‑76 (MQ=255)
                                                 aTGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACAcc                     >  1:1838459/1‑76 (MQ=255)
                                                  tGTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCa                    >  2:1342420/1‑76 (MQ=255)
                                                   gTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAg                   >  1:506740/1‑76 (MQ=255)
                                                   gTTCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAg                   >  2:472849/1‑76 (MQ=255)
                                                     tCCTCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCg                 <  1:1587153/76‑1 (MQ=255)
                                                        tCATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTg              <  2:1024511/76‑1 (MQ=255)
                                                         cATCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGc             >  1:100571/1‑76 (MQ=255)
                                                          aTCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCg            >  1:443806/1‑76 (MQ=255)
                                                           tCGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGt           <  1:677390/76‑1 (MQ=255)
                                                            cGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTg          >  2:1403463/1‑76 (MQ=255)
                                                            cGTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTg          <  1:1068831/76‑1 (MQ=255)
                                                             gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg         <  1:355813/76‑1 (MQ=255)
                                                             gTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg         <  2:459903/76‑1 (MQ=255)
                                                             gTCAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg         <  1:1295501/76‑1 (MQ=255)
                                                               cAAATGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGt       <  2:339016/76‑1 (MQ=255)
                                                                aaaTGAAGAACCCCGTTTGCGGGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGCCACCAGCGTTGCGGGGGtt      >  2:497288/1‑75 (MQ=255)
                                                                 aaTGAAGAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtat     <  1:854074/76‑1 (MQ=255)
                                                                    gaagaaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAgg  <  1:973255/76‑1 (MQ=255)
                                                                       |                                                                        
CACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGG  >  REL606/4201863‑4202006

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: