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.

ATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTA  >  REL606/498969‑499068
                                                  |                                                 
atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATc                                                   >  1:842606/1‑51 (MQ=255)
atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATc                                                   >  1:6425966/1‑51 (MQ=255)
atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATc                                                   >  1:505173/1‑51 (MQ=255)
atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATc                                                   <  1:4945642/51‑1 (MQ=255)
  gatgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAg                                                 <  1:3662882/51‑1 (MQ=255)
  gatgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAg                                                 >  1:2292102/1‑51 (MQ=255)
   atgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                <  1:1895101/51‑1 (MQ=255)
   atgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                >  1:797662/1‑51 (MQ=255)
   atgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                >  1:3750021/1‑51 (MQ=255)
   atgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                <  1:1076145/51‑1 (MQ=255)
    tgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAgtc                                               <  1:5594080/51‑3 (MQ=255)
    tgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                               >  1:2162560/1‑51 (MQ=255)
    tgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                               <  1:706813/51‑1 (MQ=255)
    tgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                               >  1:2284161/1‑51 (MQ=255)
     gaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCg                                              >  1:2570838/1‑51 (MQ=255)
      aaTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGt                                             >  1:3489946/1‑51 (MQ=255)
        tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                           <  1:5874282/51‑1 (MQ=255)
            gCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTc                                       >  1:4787537/1‑51 (MQ=255)
            gCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTc                                       >  1:6984426/1‑51 (MQ=255)
            gCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCAGCCGTACCGTc                                       >  1:5339636/1‑51 (MQ=255)
              aaCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                      >  1:7046342/1‑50 (MQ=255)
              aaCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                     <  1:6764722/51‑1 (MQ=255)
              aaCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                     >  1:1131672/1‑51 (MQ=255)
              aaCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                     <  1:3315711/51‑1 (MQ=255)
                  aGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                 <  1:4950763/51‑1 (MQ=255)
                    cTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                               <  1:3667808/51‑1 (MQ=255)
                    cTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                               <  1:7271343/51‑1 (MQ=255)
                     ttAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                              >  1:7034836/1‑51 (MQ=255)
                     ttAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                              >  1:3023079/1‑51 (MQ=255)
                     ttAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                              >  1:3076950/1‑51 (MQ=255)
                      tAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca                             >  1:411174/1‑51 (MQ=255)
                       aaTGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                            <  1:3487506/51‑1 (MQ=255)
                       aaTGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                            >  1:6211428/1‑51 (MQ=255)
                       aaTGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                            >  1:2196390/1‑51 (MQ=255)
                        aTGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagc                           >  1:3517115/1‑51 (MQ=255)
                         tGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                          >  1:2201027/1‑51 (MQ=255)
                          gaggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                         <  1:410841/51‑1 (MQ=255)
                          gaggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                         >  1:1752858/1‑51 (MQ=255)
                          gaggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                         <  1:4594408/51‑1 (MQ=255)
                           aggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagc                        <  1:5491265/51‑1 (MQ=255)
                           aggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagc                        <  1:6471048/51‑1 (MQ=255)
                            ggagAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagca                       >  1:3082261/1‑51 (MQ=255)
                             gagaACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCat                      <  1:3155883/51‑1 (MQ=255)
                              agaACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                     <  1:480622/51‑1 (MQ=255)
                              agaACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                     >  1:6509572/1‑51 (MQ=255)
                              agaACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                     >  1:4258183/1‑51 (MQ=255)
                               gaACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                    >  1:751317/1‑51 (MQ=255)
                                aaCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                   >  1:1587054/1‑51 (MQ=255)
                                aaCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                   >  1:5177974/1‑51 (MQ=255)
                                 aCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                  >  1:5074052/1‑51 (MQ=255)
                                 aCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                  <  1:2419057/51‑1 (MQ=255)
                                  cccGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                    <  1:7155085/48‑1 (MQ=255)
                                   ccGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                              <  1:4817536/37‑1 (MQ=255)
                                   ccGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                <  1:6544899/51‑1 (MQ=255)
                                   ccGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                <  1:5011415/51‑1 (MQ=255)
                                   ccGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                >  1:4622258/1‑51 (MQ=255)
                                    cGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg               >  1:5054359/1‑51 (MQ=255)
                                    cGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg               >  1:2011372/1‑51 (MQ=255)
                                    cGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg               >  1:3893847/1‑51 (MQ=255)
                                      catcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa             <  1:1047800/51‑1 (MQ=255)
                                      catcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa             <  1:4645885/51‑1 (MQ=255)
                                      catcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa             >  1:7198066/1‑51 (MQ=255)
                                      catcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa             <  1:3587960/51‑1 (MQ=255)
                                       atcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa             <  1:4846423/50‑1 (MQ=255)
                                       atcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGANt            <  1:4582704/51‑1 (MQ=255)
                                       atcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt            >  1:388319/1‑51 (MQ=255)
                                        tcatGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc           <  1:1219834/51‑1 (MQ=255)
                                         catGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                    <  1:2279331/41‑1 (MQ=255)
                                         catGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa          >  1:2557839/1‑51 (MQ=255)
                                         catGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa          <  1:7175972/51‑1 (MQ=255)
                                         catGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa          >  1:822706/1‑51 (MQ=255)
                                          atGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg         >  1:2914491/1‑51 (MQ=255)
                                          atGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg         >  1:2649318/1‑51 (MQ=255)
                                          atGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg         <  1:3254989/51‑1 (MQ=255)
                                          atGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg         >  1:6976132/1‑51 (MQ=255)
                                          atGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg         >  1:6073357/1‑51 (MQ=255)
                                           tgtgTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg        >  1:5680662/1‑51 (MQ=255)
                                             tgtATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGc      >  1:498220/1‑51 (MQ=255)
                                                aTCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGt   >  1:1992370/1‑51 (MQ=255)
                                                 tCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa  >  1:2021459/1‑51 (MQ=255)
                                                  |                                                 
ATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTA  >  REL606/498969‑499068

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