Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,201,958 A→C 46.0% L201R (CTC→CGC)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,201,9580A→C46.0% 211.3 ‑0.3 ‑0.0133L201R (CTC→CGC) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (32/29):  ref base (41/31):  total (73/60)
Fisher's exact test strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.50e-01

Only 100 of 152 total aligned reads displayed.

GTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGT  >  REL606/4201924‑4201992
                                  |                                  
 tCAAATGAATAACCCCGTTTGCGCGTTTGGGCgag                                   <  1:5684603/35‑1 (MQ=255)
 tCAAATGAATAACCCCGTTTGCGCGTTTGGGCgag                                   >  1:15673927/1‑35 (MQ=255)
  cAAATGAATAACCCCGTTTGCGCGTTTGGGCgcga                                  <  1:3160130/35‑1 (MQ=255)
  cAAATGAATAACCCCGTTTGCGCGTTTGGGCgcga                                  >  1:5526516/1‑35 (MQ=255)
  cAAATGAATAACCCCGTTTGCGCGTTTGGGCgaga                                  <  1:947416/35‑1 (MQ=255)
   aaaTGAATAACCCCGTTTGCGCGTTTGGGCGCGAt                                 <  1:8044341/35‑1 (MQ=255)
   aaaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAt                                 >  1:2187013/1‑35 (MQ=255)
   aaaTGAATAACACCGTTTGCGCGTTTGGGCGAGAt                                 <  1:15101514/35‑1 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtc                                <  1:10689155/35‑1 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtc                                <  1:1034932/35‑1 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtc                                <  1:16070638/35‑1 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtc                                >  1:10048042/1‑35 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtc                                >  1:9113656/1‑35 (MQ=255)
    aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtc                                >  1:17696101/1‑35 (MQ=255)
     aTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtct                               >  1:8435685/1‑35 (MQ=255)
     aTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtct                               <  1:11171509/35‑1 (MQ=255)
     aTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtct                               >  1:3177339/1‑35 (MQ=255)
     aTGAATAACCCCGTTTGCGCGTTTGGGCGCGAtct                               <  1:2434270/35‑1 (MQ=255)
     aTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtct                               <  1:12790254/35‑1 (MQ=255)
      tGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                              >  1:7074421/1‑35 (MQ=255)
      tGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                              >  1:15128884/1‑35 (MQ=255)
      tGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                              <  1:6786283/35‑1 (MQ=255)
       gAATAACCCCGTTTGCGCGTTTGGGCGCGAtcttc                             >  1:18379802/1‑35 (MQ=255)
       gAATAACCCCGTTTGCGCGTTTGGGCGCGAtcttc                             >  1:13846931/1‑35 (MQ=255)
       gAATAACCCCGTTTGCGCGTTTGGGCGCGAtcttc                             >  1:14597167/1‑35 (MQ=255)
       gAATAACCCCGTTTGCGCGTTTGGGCGCGAtcttc                             >  1:11061163/1‑35 (MQ=255)
       gAATAACCCCGTTTGCGCGTTTGGGCGCGAtcttc                             >  1:10036781/1‑35 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGCGAtcttct                            <  1:13079143/35‑1 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGCGAtcttct                            >  1:5172734/1‑35 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                            <  1:15836529/35‑1 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                            >  1:15522906/1‑35 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                            >  1:5116902/1‑35 (MQ=255)
        aaTAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                            >  1:8599398/1‑35 (MQ=255)
         aTAACCCCGTTTGCGCGTTTGGGCGCGATcttctt                           <  1:16020841/35‑1 (MQ=255)
         aTAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                           <  1:12111116/35‑1 (MQ=255)
         aTAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                           <  1:6286701/35‑1 (MQ=255)
          tAACCCCGTTTGCGCGTTTGGGCGCGATCTTCttt                          >  1:13244049/1‑35 (MQ=255)
          tAACCCCGTTTGCGCGTTTGGGCGCGATCTTCttt                          >  1:12355215/1‑35 (MQ=255)
          tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCttt                          >  1:1218911/1‑35 (MQ=255)
           aaCCCCGTTTGCGCGTTTGGGCGCGATCTTCtttt                         >  1:7638852/1‑35 (MQ=255)
            aCCCCGTTTGCGCGTTTGGGCGCGATCTTCTTTTa                        <  1:11948566/35‑1 (MQ=255)
            aCCCCGTTTGCGCGTTTGGGCGCGATCTTCTTTTa                        >  1:19057303/1‑35 (MQ=255)
            aCCCCGTTTGCGCGTTTGGGCGCGATCTTCTTTTa                        <  1:15426810/35‑1 (MQ=255)
            aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTa                        <  1:2809617/35‑1 (MQ=255)
            aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTa                        <  1:14633335/35‑1 (MQ=255)
             ccccGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa                       <  1:19222562/35‑1 (MQ=255)
             ccccGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa                       <  1:13214902/35‑1 (MQ=255)
             ccccGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa                       <  1:10122924/35‑1 (MQ=255)
              cccGTTTGCGCGTTTGGGCGCGATCTTCTTTTaaa                      <  1:15561512/35‑1 (MQ=255)
              cccGTTTGCGCGTTTGGGCGCGATCTTCTTTTaaa                      >  1:2117423/1‑35 (MQ=255)
              cccGTTTGCGCGTTTGGGCGCGATCTTCTTTTaaa                      >  1:8160230/1‑35 (MQ=255)
              cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTaaa                      >  1:12687738/1‑35 (MQ=255)
               ccGTTTGCGCGTTTGGGCGCGATCTTCTTTTAAAt                     <  1:7815583/35‑1 (MQ=255)
               ccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAAt                     >  1:3234412/1‑35 (MQ=255)
                cGTTTGCGCGTTTGTGCGAGATCTTCTTTTAAATg                    >  1:3459739/1‑35 (MQ=255)
                cGTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATg                    >  1:12351479/1‑35 (MQ=255)
                cGTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATg                    >  1:2255679/1‑35 (MQ=255)
                 gTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATGc                   <  1:5405813/35‑1 (MQ=255)
                 gTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATGc                   >  1:9917047/1‑35 (MQ=255)
                 gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc                   <  1:1042310/35‑1 (MQ=255)
                 gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc                   >  1:874126/1‑35 (MQ=255)
                 gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc                   >  1:7049551/1‑35 (MQ=255)
                  tttGCGCGTTTGGGCGCGATCTTCTTTTAAATGca                  >  1:2735986/1‑35 (MQ=255)
                  tttGCGCGTTTGGGCGCGATCTTCTTTTAAATGca                  >  1:17824136/1‑35 (MQ=255)
                  tttGCGCGTTTGGGCGATATCTTCTTTTAAATGca                  >  1:15457166/1‑35 (MQ=255)
                  tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGca                  <  1:16647893/35‑1 (MQ=255)
                  tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGca                  <  1:6968567/35‑1 (MQ=255)
                  tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGca                  >  1:11733822/1‑35 (MQ=255)
                  tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGca                  <  1:10613577/35‑1 (MQ=255)
                   ttGCGCGTTTGGGCGCGATCTTCTTTTAAATGcac                 >  1:5374833/1‑35 (MQ=255)
                   ttGCGCGTTTGGGCGCGATCTTCTTTTAAATGcac                 <  1:17160965/35‑1 (MQ=255)
                   ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac                 >  1:18890490/1‑35 (MQ=255)
                    tGCGCGTTTGGGCGCGATCTTCTTTTAAATGcaca                <  1:20125693/35‑1 (MQ=255)
                    tGCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca                <  1:1979289/35‑1 (MQ=255)
                     gcgcGTTTGGGCGCGATCTTCTTTTAAATGCACAg               <  1:11984852/35‑1 (MQ=255)
                      cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAgg              <  1:7978130/35‑1 (MQ=255)
                      cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAgg              <  1:3019349/35‑1 (MQ=255)
                       gcgTTTGGGCGCGATCTTCTTTTAAATGCACAGga             >  1:787036/1‑35 (MQ=255)
                       gcgTTTGGGCGAGATCTTCTTTTAAATGCACAGga             >  1:10918727/1‑35 (MQ=255)
                        cgTTTGGGCGAGATCTTCTTTTAAATGCACAGgag            >  1:102531/1‑35 (MQ=255)
                         gTTTGGGCGCGATCTTCTTTTAAATGCACAGgaga           <  1:8280885/35‑1 (MQ=255)
                         gTTTGGGCGCGATCTTCTTTTAAATGCACAGgaga           >  1:4042418/1‑35 (MQ=255)
                         gTTTGGGCGAGATCTTCTTTTAAATGCACAGgaga           >  1:9644467/1‑35 (MQ=255)
                         gTTTGGGCGAGATCTTCTTTTAAATGCACAGgaga           >  1:8907933/1‑35 (MQ=255)
                         gTTTGGGCGAGATCTTCTTTTAAATGCACAGgaga           >  1:6533552/1‑35 (MQ=255)
                          tttGGGCGCGATCTTCTTTTAAATGCACAGTAGac          >  1:163663/1‑35 (MQ=38)
                          tttGGGCGCGATCTTCTTTTAAATGCACAGGAGac          >  1:9405700/1‑35 (MQ=255)
                          tttGGGCGAGATCTTCTTTTAAACGCACAGGAGac          >  1:3166400/1‑35 (MQ=255)
                           ttGGGCGCGATCTTCTTTTAAATGCACAGGAGaca         <  1:7962023/35‑1 (MQ=255)
                             gggCGAGATCTTCTTTTAAATGCACAGGAGACAcc       >  1:16906553/1‑35 (MQ=255)
                             gggCGAGACCTTCTTTTAAATGCACAGGAGACAcc       <  1:4680940/35‑1 (MQ=255)
                               gCGCGATCTTCTTTTAAATGCACAGGAGACACCAg     >  1:10912628/1‑35 (MQ=255)
                               gCGCGATCTTCTTTTAAATGCACAGGAGACACCAg     <  1:9213498/35‑1 (MQ=255)
                               gCGAGATCTTCTTTTAAATGCACAGGAGACACCAg     >  1:11728469/1‑35 (MQ=255)
                               gCGAGATCTTCTTTTAAATGCACAGGAGACACCAg     <  1:2961272/35‑1 (MQ=255)
                               gCGAGATCTTCTTTTAAATGCACAGGAGACACCAg     <  1:5997777/35‑1 (MQ=255)
                                cGAGATCTTCTTTTAAATGCACAGGAGACACCAGc    <  1:3493636/35‑1 (MQ=255)
                                 gagaTCTTCTTTTAAATGCACAGGAGACACCAGCg   <  1:13751806/35‑1 (MQ=255)
                                 gagaTCTTCTTTTAAATGCACAGGAGACACCAGCg   <  1:7369989/35‑1 (MQ=255)
                                  agaTCTTCTTTTAAATGCACAGGAGACACCAGCGt  <  1:18314952/35‑1 (MQ=255)
                                  |                                  
GTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGT  >  REL606/4201924‑4201992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.