New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 499019 =80 (1.000)10 (0.140) 6/92 8.6 13.8% coding (3592/4281 nt) rhsD rhsD element protein
?REL606 = 3099737 51 (0.700)coding (409/2259 nt) parC DNA topoisomerase IV subunit A
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGT  >  REL606/3099687‑3099787
                                                  |                                                  
cTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGccc                                                    <  1:7235007/51‑1 (MQ=255)
cTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGccc                                                    >  1:6263903/1‑51 (MQ=255)
 tCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGcccc                                                   >  1:7244295/1‑51 (MQ=255)
 tCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGcccc                                                   >  1:5267574/1‑51 (MQ=255)
    tGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgc                                                >  1:3923834/1‑51 (MQ=255)
    tGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgc                                                >  1:5249288/1‑51 (MQ=255)
      cAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctc                                              <  1:4805553/51‑1 (MQ=255)
      cAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctc                                              <  1:7043927/51‑1 (MQ=255)
      cAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctc                                              >  1:3352160/1‑51 (MQ=255)
       aaaGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcg                                             >  1:6368425/1‑51 (MQ=255)
       aaaGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcg                                             >  1:2011573/1‑51 (MQ=255)
       aaaGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcg                                             <  1:1201604/51‑1 (MQ=255)
       aaaGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcg                                             <  1:2280067/51‑1 (MQ=255)
         aGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcgct                                           >  1:1120136/1‑51 (MQ=255)
           tGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCa                                         <  1:1201955/51‑1 (MQ=255)
            gCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCaa                                        <  1:1701782/51‑1 (MQ=255)
             ccGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAAt                                       <  1:2539994/51‑1 (MQ=255)
               gTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATag                                     <  1:2372061/51‑1 (MQ=255)
                    aGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCt                                <  1:7225351/51‑1 (MQ=255)
                    aGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCt                                <  1:5887945/51‑1 (MQ=255)
                       ttGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCgg                             >  1:4259458/1‑51 (MQ=255)
                        tGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGa                            <  1:4503589/51‑1 (MQ=255)
                        tGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGa                            >  1:3673665/1‑51 (MQ=255)
                          gCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAat                          >  1:5977015/1‑51 (MQ=255)
                           cACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAata                         >  1:6235412/1‑51 (MQ=255)
                           cACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAata                         >  1:7173830/1‑51 (MQ=255)
                            aCCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAatat                        <  1:2243179/51‑1 (MQ=255)
                              ccAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATAttt                      <  1:1139814/51‑1 (MQ=255)
                               cAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTc                     <  1:6641865/51‑1 (MQ=255)
                                aGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCg                    >  1:155733/1‑51 (MQ=255)
                                aGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCg                    >  1:3397163/1‑51 (MQ=255)
                                 gTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGa                   <  1:3704107/51‑1 (MQ=255)
                                  tCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTC‑‑AATATTTCGACaa                <  1:1571390/51‑1 (MQ=255)
                                    aGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACaa                <  1:6972437/51‑1 (MQ=255)
                                     gCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAc               <  1:6223238/51‑1 (MQ=255)
                                     gCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAc               >  1:4305221/1‑51 (MQ=255)
                                      ccGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAcc              <  1:2243142/51‑1 (MQ=255)
                                       cGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTTGGAATATTTCGACAACCg             <  1:4395622/51‑1 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            >  1:5539960/1‑51 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            <  1:1415291/51‑1 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            <  1:3628145/51‑1 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            <  1:59994/51‑1 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            <  1:4081473/51‑1 (MQ=255)
                                        gTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg            >  1:6852609/1‑51 (MQ=255)
                                         tCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCggg           >  1:268828/1‑51 (MQ=255)
                                         tCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCggg           >  1:164793/1‑51 (MQ=255)
                                         tCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCggg           >  1:1236671/1‑51 (MQ=255)
                                          ccccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGa          <  1:3248219/51‑1 (MQ=255)
                                          ccccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGa          <  1:2625472/51‑1 (MQ=255)
                                          ccccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGa          >  1:3128896/1‑51 (MQ=255)
                                           cccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAt         <  1:38443/51‑1 (MQ=255)
                                            ccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAtt        >  1:7040245/1‑51 (MQ=255)
                                            ccTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAtt        >  1:2274500/1‑51 (MQ=255)
                                              tGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCg      >  1:7071738/1‑51 (MQ=255)
                                                ccccAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGgt    <  1:6432197/51‑1 (MQ=255)
                                                 cccAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGgtg   <  1:4094020/51‑1 (MQ=255)
                                                  ccAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGgtgt  >  1:2403436/1‑51 (MQ=255)
                                                  |                                                  
CTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGT  >  REL606/3099687‑3099787

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