Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,462,320 IS150 (–) +3 bp 90.4% intergenic (‑82/‑257) mokB ← / → trg regulatory peptide/methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 664688 =NA (NA)66 (0.700) 36/44 NT 89.1% noncoding (1/1443 nt) IS150 repeat region
?REL606 1462320 = 7 (0.090)intergenic (‑82/‑259) mokB/trg regulatory peptide/methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
* ? REL606 = 666130NA (NA)66 (0.640) 36/48 NT 89.2% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1462322 7 (0.080)intergenic (‑84/‑257) mokB/trg regulatory peptide/methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

TTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCG  >  REL606/666101‑666164
                              |                                  
ttAACTGTCCAACTTTTTGGTGTCAGTACATTTCa                               >  1:2745646/1‑35 (MQ=255)
ttAACTGTCCAACTTTTTGGGGTCAGTACATTTcg                               >  1:17977660/1‑34 (MQ=255)
ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                               >  1:6292258/1‑35 (MQ=255)
 tAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                              <  1:6953889/35‑1 (MQ=255)
 tAACTGTACAACTTTTTGGGGTCAGTACATTTCAg                              <  1:7191247/35‑1 (MQ=255)
  aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAgg                             >  1:10974702/1‑35 (MQ=255)
  aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAgg                             <  1:18658842/35‑1 (MQ=255)
   aCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                            <  1:14159931/35‑1 (MQ=255)
   aCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                            >  1:2794403/1‑35 (MQ=255)
   aCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                            <  1:2304427/35‑1 (MQ=255)
   aCTGTCCAACTTTTTGGGGTCAGTACAATCCAggg                            <  1:2435937/35‑1 (MQ=11)
    cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                           >  1:13239558/1‑35 (MQ=255)
    cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                           >  1:2896254/1‑35 (MQ=255)
    cTGTCCAACTTTTTGGGGTCAGTACAATCCAGGGc                           >  1:12173974/1‑35 (MQ=11)
     tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                          >  1:19145898/1‑35 (MQ=255)
     tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                          >  1:18819481/1‑35 (MQ=255)
     tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                          >  1:3487247/1‑35 (MQ=255)
     tGTCCAACTTTTTGGGGTCAGTACAATCCAGGGCt                          >  1:16227253/1‑35 (MQ=14)
      gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                         <  1:12367228/35‑1 (MQ=255)
      gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                         >  1:20274748/1‑35 (MQ=255)
       tCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                        >  1:12797761/1‑35 (MQ=255)
       tCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                        <  1:10500762/35‑1 (MQ=255)
        ccAACTTTTTGGGGTCAGTACATTTCAGGGTtttt                       <  1:17975374/35‑1 (MQ=255)
        ccAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                       >  1:15357225/1‑35 (MQ=255)
         cAACTTTTTGGGGTCAGTACATTTCAGGGCttttt                      <  1:7000530/35‑1 (MQ=255)
          aaCTTTTTGGGGTCAGTACAACTTCGGGGCttttc                      >  1:18506961/1‑34 (MQ=0)
          aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                     <  1:6421696/35‑1 (MQ=255)
          aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                     >  1:16609813/1‑35 (MQ=255)
          aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                     <  1:12886781/35‑1 (MQ=255)
            cTTTTTGGGGTCAGTACATTTCAGGGCtttttaga                   >  1:8615229/1‑32 (MQ=255)
             tttttGGGGTCAGTACATTTCAGGGCTTTTTCACg                  <  1:1032592/35‑1 (MQ=255)
             tttttGGGGTCAGTACATTTCAGGGCTTTTTCACg                  <  1:8881663/35‑1 (MQ=255)
             tttttGGGGTCAGTACATTTCAGGGCTTTTTCACg                  >  1:14646941/1‑35 (MQ=255)
              tcttGGGGTCAGTACATTTCAGGGCTTTTTCACGa                 <  1:59424/33‑1 (MQ=255)
               tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                >  1:6581808/1‑35 (MQ=255)
               tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                <  1:5393534/35‑1 (MQ=255)
               tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                >  1:2310492/1‑35 (MQ=255)
               tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                <  1:17692511/35‑1 (MQ=255)
                ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc               >  1:5323053/1‑35 (MQ=255)
                ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc               <  1:13494482/35‑1 (MQ=255)
                ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc               >  1:1959364/1‑35 (MQ=255)
                 tGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc              <  1:11717832/35‑1 (MQ=255)
                 tGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc              >  1:8607355/1‑35 (MQ=255)
                   gggTCAGTACATTTCAGGGCTTTTTCACGAGCCgg            <  1:6969402/35‑1 (MQ=255)
                   gggTCAGTACATTTCAGGGCTTTTTCACGAGCCgg            <  1:19255583/35‑1 (MQ=255)
                   gggTCAGTACATTTAAGGGCTTTTTCACGAGCCgg            <  1:11595856/35‑1 (MQ=255)
                    ggTCAGTACATTTCAGTGCTTTTTCACGAGCCGgc           >  1:279280/1‑35 (MQ=255)
                    ggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgc           >  1:11910738/1‑35 (MQ=255)
                     gTCAGTACATTTCAGGGTTTTTTCACGAGCCGgcg          <  1:13384539/35‑1 (MQ=255)
                     gTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcg          >  1:4135610/1‑35 (MQ=255)
                      tCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc         >  1:882691/1‑35 (MQ=255)
                      tCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc         >  1:6462808/1‑35 (MQ=255)
                      tCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc         >  1:19570185/1‑35 (MQ=255)
                      tCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc         <  1:12229038/35‑1 (MQ=255)
                      tCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc         >  1:10697752/1‑35 (MQ=255)
                       cAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCt        >  1:17727729/1‑35 (MQ=255)
                        aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCtt       >  1:20283801/1‑35 (MQ=255)
                         gTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt      >  1:5555404/1‑35 (MQ=255)
                         gTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt      >  1:9626258/1‑35 (MQ=255)
                              aTTTCAGGGTTTTTTCGCGAGCCGGCGCTTTATcg  <  1:15646362/35‑1 (MQ=38)
                              aTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg  <  1:4457560/35‑1 (MQ=255)
                              aTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg  <  1:8097324/35‑1 (MQ=255)
                              aTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg  >  1:10494273/1‑35 (MQ=255)
                              |                                  
TTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCG  >  REL606/666101‑666164

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