Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,202,088 C→T 5.0% A158T (GCT→ACT)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,202,0880C→T 5.0% 5.1 ‑0.4 ‑0.1123A158T (GCT→ACT) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (69/47):  total (73/50)
Fisher's exact test strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.28e-01

Only 100 of 130 total aligned reads displayed.

AGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGAC  >  REL606/4202054‑4202122
                                  |                                  
aGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGc                                    >  1:19817348/1‑35 (MQ=255)
aGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGc                                    >  1:2832308/1‑35 (MQ=255)
 gTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCg                                   <  1:9976003/35‑1 (MQ=255)
 gTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCg                                   <  1:893080/35‑1 (MQ=255)
 gTTGGGCTAAAAAGGCTTACCTGCACCGGAAGCgg                                  <  1:2091891/35‑2 (MQ=25)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGTGt                                  <  1:14741836/35‑1 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  <  1:15153619/35‑1 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  <  1:2878228/35‑1 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  >  1:1215200/1‑35 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  >  1:6550253/1‑35 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  <  1:7790973/35‑1 (MQ=255)
  ttGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGt                                  <  1:19227405/35‑1 (MQ=255)
   tGGGCTAAAAAGGCTTTACCTGCACCGGAAGTGTg                                 <  1:1834325/35‑1 (MQ=255)
   tGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTg                                 <  1:7678052/35‑1 (MQ=255)
    gggCTAAAAAGGCTTTACCTGCACCGGAAGCGTGc                                <  1:784335/35‑1 (MQ=255)
    gggCTAAAAAGGCTTTACCTGCACCGGAAGCGTGc                                <  1:6036805/35‑1 (MQ=255)
     ggCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCa                               >  1:13239833/1‑35 (MQ=255)
     ggCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCa                               <  1:14939768/35‑1 (MQ=255)
      gCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCAt                              >  1:11050726/1‑35 (MQ=255)
      gCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCAt                              >  1:844746/1‑35 (MQ=255)
      gCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCAt                              >  1:8557804/1‑35 (MQ=255)
       cTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATc                             <  1:15901461/35‑1 (MQ=255)
       cTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATc                             <  1:4292402/35‑1 (MQ=255)
        tAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCg                            >  1:15189569/1‑35 (MQ=255)
        tAAAAAGGCTTTACCTGCACAGGAAGCGTGCATCg                            <  1:6532774/35‑1 (MQ=255)
         aaaaaGGCTTTACCTGCACCGGAAGCGTGCATCgg                           >  1:11951731/1‑35 (MQ=255)
         aaaaaGGCTTTACCTGCACCGGAAGCGTGCATCgg                           <  1:1023236/35‑1 (MQ=255)
         aaaaaGGCTTTACCTGCACCGGAAGCGTGCATCgg                           <  1:20021050/35‑1 (MQ=255)
          aaaaGGCTTTACCTGCACCGGAAGCGTGCATCGGc                          >  1:18934852/1‑35 (MQ=255)
          aaaaGGCTTTACCTGCACCGGAAGCGTGCATCGGc                          >  1:17419520/1‑35 (MQ=255)
          aaaaGGCTTTACCTGCACCGGAAGCGTGCATCGGc                          >  1:12367876/1‑35 (MQ=255)
          aaaaGGCTTTACCTGCACCGGAAGCGTGCATCGGc                          <  1:1205199/35‑1 (MQ=255)
           aagggCTTTCCCTGCACCGGAAGCGTGCATCGGCa                         <  1:18441453/32‑1 (MQ=25)
           aaaGGCTTTACCTGCACCGGAAGCGTGCATCGGCa                         <  1:14025693/35‑1 (MQ=255)
           aaaGGCTTTACCTGCACCGGAAGCGTGCATCGGCa                         <  1:7446706/35‑1 (MQ=255)
           aaaGGCTTTACCTGCACCGGAAGCGTGCATCGGCa                         >  1:1537377/1‑35 (MQ=255)
           aaaGGCTTTACCTGCACCGGAAGCGTGCATCGGCa                         >  1:10205011/1‑35 (MQ=255)
           aaaGGCTTTACCTGCACCGGAAGCGTACATCGGCa                         <  1:13718748/35‑1 (MQ=255)
            aaGGCTTTACCTGCACCGGAAGCGTGCATCGGCaa                        <  1:5367690/35‑1 (MQ=255)
            aaGGCTTTACCTGCACCGGAAGCGTGCATCGGCaa                        >  1:11319538/1‑35 (MQ=255)
            aaGGCTTTACCTGCACCGGAAGCGTGCATCGGCaa                        >  1:9253865/1‑35 (MQ=255)
             aGGCTTTACCTGCACCGGAAGCGTGCATCGGCAAt                       >  1:2477539/1‑35 (MQ=255)
             aGGCTTTACCTGCACCGGAAGCGTGCATCGGCAAt                       <  1:9731819/35‑1 (MQ=255)
             aGGCTTTACCTGCACCGGAAGCGTGCATCGGCAAt                       >  1:10068337/1‑35 (MQ=255)
             aGGCTTTACCTGCACCGGAAGCGTGCATCGGCAAt                       <  1:13472593/35‑1 (MQ=255)
              ggCTTTCCCTGCACCGGAAGCGTGCATCGGCAAtt                      <  1:12175917/35‑1 (MQ=255)
              ggCTTTACCTGCACCGGAAGCGTGCATCGGCAAtt                      >  1:10735690/1‑35 (MQ=255)
              ggCTTTACCTGCACCGGAAGCGTGCATCGGCAAtt                      <  1:475652/35‑1 (MQ=255)
              ggCTTTACCTGCACCGGAAGCGTGCATCGGCAAtt                      >  1:3149706/1‑35 (MQ=255)
               gCTTTACCTGCACCGGAAGCGTGCATCGGCAAttt                     >  1:11282041/1‑35 (MQ=255)
               gCTTTACCTGCACCGGAAGCGTGCATCGGCAAttt                     >  1:714650/1‑35 (MQ=255)
               gCTTTACCTGCACCGGAAGCGTGCATCGGCAAttt                     >  1:19782213/1‑35 (MQ=255)
               gCTTTACCTGCACCGGAAGCGTGCATCGGCAAtct                     >  1:10737886/1‑33 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    <  1:7894949/35‑1 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    <  1:10809/35‑1 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    <  1:7245189/35‑1 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    <  1:3701519/35‑1 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    >  1:12830551/1‑35 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    >  1:17474664/1‑35 (MQ=255)
                cTTTACCTGCACCGGAAGCGTGCATCGGCAATTTa                    >  1:19162648/1‑35 (MQ=255)
                 tttACCTGCACCGGAAGCGTGCATCGGCAATTTAc                   >  1:20411270/1‑35 (MQ=255)
                 tttACCTGCACCGGAAGCGTGCATCGGCAATTTAc                   >  1:10615945/1‑35 (MQ=255)
                 tttACCTGCACCGGAAGCGTGCATCGGCAATTTAc                   <  1:11452697/35‑1 (MQ=255)
                 tttACCTGCACCGGAAGCGTGCATCGGCAATTTAc                   >  1:12607929/1‑35 (MQ=255)
                 tttACCTGCACCGGAAGCGTGCATCGGCAATTTAc                   <  1:497554/35‑1 (MQ=255)
                  ttACCTGCACCGGAAGTGTGCATAGGCAATTTAcc                  <  1:2893328/35‑1 (MQ=25)
                  ttACCTGCACCGGAAGCGTGCATCGGCAATTTAcc                  <  1:1552415/35‑1 (MQ=255)
                  ttACCTGCACCGGAAGCGTGCATCGGCAATTTAcc                  >  1:14707069/1‑35 (MQ=255)
                  ttACCTGCACCGGAAGCGTGCATCGGCAATTTAcc                  <  1:7333392/35‑1 (MQ=255)
                  ttACCTGCACCGGAAGCGTGCATCGGCAATTTAcc                  >  1:6964803/1‑35 (MQ=255)
                  ttACCTGCACCGGAAGCGTGCATCGGCAATTTAcc                  <  1:642926/35‑1 (MQ=255)
                  ttACCTGCACCGGAAGCGTGCATCGGCAAATTAcc                  >  1:1389275/1‑35 (MQ=255)
                   tACCTGCACCGGAAGCGTGCATCGGCAATTTAccg                 >  1:4083122/1‑35 (MQ=255)
                   tACCTGCAACGGAAGCGTGCATCGGCAATTTAccg                 <  1:19146588/35‑1 (MQ=255)
                    aCCTGCACCGGAAGCGTGCATCGGCAATTTAccgc                <  1:12387830/35‑1 (MQ=255)
                    aCCTGCACCGGAAGCGTGCATCGGCAATTTAccgc                <  1:5529156/35‑1 (MQ=255)
                    aCCTGCACCGGAAGCGTGCATCGGCAATTTAccgc                >  1:4776531/1‑35 (MQ=255)
                     ccTGCACCGGAAGCGTGCATCTGCAATTTAccgcc               >  1:9416720/1‑35 (MQ=255)
                      cTGCACCGGAAGCGTGCATCGGCAATTTAccgccg              >  1:9496919/1‑35 (MQ=255)
                       tGCACCGGAAGCGTGCATCGGCAATTTACCGCCGa             >  1:2025927/1‑35 (MQ=255)
                       tGCACCGGAAGCGTGCATCGGCAATTTACCGCCGa             >  1:754437/1‑35 (MQ=255)
                       tGCACCGGAAGCGTGCATCGGCAATTTACCGCCGa             >  1:12074564/1‑35 (MQ=255)
                       tGCACCGGAAGCGTGCATCGGCAATTTACCGCCGa             >  1:5313986/1‑35 (MQ=255)
                        gCACCGGAAGCGTGCATCGGCAATTTACCGCCGAt            >  1:10580629/1‑35 (MQ=255)
                         cACCGGAAGCGTGCATCGGCAATTTACCGCCGATa           >  1:10561026/1‑35 (MQ=255)
                         cACCGGAAGCGTGCATCGGCAATTTACCGCCGATa           >  1:8227554/1‑35 (MQ=255)
                          aCCGGAAGCGTGCATCGGCAAATTACCGCCGATAg          >  1:15289984/1‑35 (MQ=255)
                           ccGGAAGTGTGCATCGGCAATTTACCGCCGATAgg         >  1:18886647/1‑35 (MQ=255)
                           ccGGAAGTGTGCATCGGCAATTTACCGCCGATAgg         >  1:12859433/1‑35 (MQ=255)
                           ccGGAAGCGTGCATCGGCAATTTACCGCCGATAgg         >  1:1802008/1‑35 (MQ=255)
                           ccGGAAGCGTGCATCGGCAATTTACCGCCGATAgg         >  1:7390613/1‑35 (MQ=255)
                            cGGAAGTGTGCATCGGCAATTTACCGCCGATAGgc        >  1:1028823/1‑35 (MQ=255)
                            cGGAAGCGTGCATCGGCAATTTACCGCCGATAGgc        >  1:7692835/1‑35 (MQ=255)
                            cGGAAGCGTGCATCGGCAATTTACCGCCGATAGgc        >  1:14491088/1‑35 (MQ=255)
                             ggAAGCGTGCATCGGCAATTTACCGCCGATAGgcg       >  1:9141815/1‑35 (MQ=255)
                             ggAAGCGTGCATCGGCAATTTACCGCCGATAGgcg       >  1:12158348/1‑35 (MQ=255)
                                aGCGTGCATCGGCAATTTACCGCCGATAGGCGCgg    >  1:15359784/1‑35 (MQ=255)
                                aGCGTGCATCGGCAATTTACCGCCGATAGGCGCgg    >  1:7034793/1‑35 (MQ=255)
                                 gCGTGCATCGGCAATTTACCGCCGATAGGGGCGGa   >  1:16804093/1‑35 (MQ=255)
                                  cGTGCATCGGCAATTTACCGCCGATAGGCGCGGAc  >  1:1800035/1‑35 (MQ=255)
                                  |                                  
AGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGAC  >  REL606/4202054‑4202122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.