New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1661078 =57 (0.550)10 (0.120) 4/40 NT 15.8% coding (1002/1302 nt) rstB sensory histidine kinase in two‑component regulatory system with RstA
?REL606 = 3625456 59 (0.690)pseudogene (7/738 nt) yhjQ b3534; cell division protein (chromosome partitioning ATPase)

CCGACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAG‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1661115‑1661078
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagGGGGGGCGGGGA  <  REL606/3625456‑3625445
                                                  
CCGACAGTAGCAGGCTGGTTTCAACCGTTGAATGG                 <  1:12173529/35‑1
CCGACAGTAGCAGGCTGGTTTCAACCGTTGAATGG                 >  1:18633453/1‑35
 CGACAGTAGCAGGCTTGTTTCAACCGTTGAATGGC                >  1:258731/1‑35
 CGACAGTAGCAGGCTGGTTTCAACCGTTGAATGGC                >  1:1753879/1‑35
  GACAGTAGCATGCTGGTTTCAACCGTTGAATGGCA               >  1:576126/1‑35
  GACAGTAGCAGGCTTGTTTCAACCGTTGAATGGCA               >  1:1640734/1‑35
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               >  1:4648859/1‑35
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:4095948/35‑1
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               >  1:3496051/1‑35
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:17953481/35‑1
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:8203768/35‑1
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:15133674/35‑1
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:14774035/35‑1
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               >  1:8297455/1‑35
  GACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCA               <  1:12977866/35‑1
   ACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAG              <  1:7832799/35‑1
    CAGTAGCAGGCTTGTTTCAACCGTTGAATGGCAGG             >  1:5768226/1‑35
    CAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAGG             >  1:20428485/1‑35
    CAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAGG             >  1:1300595/1‑35
    CAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAGG             >  1:971121/1‑35
     AGTAGCAGGCTGGTTTCAACCGTTGAATGGCAGGG            >  1:16212311/1‑35
     AGTAGCAGGCTGGTTTCAACCGTTGAATGGCAGGG            >  1:18087011/1‑35
      GTAGCAGGCTGGTTTCAACCGTTGAATGGCAGGGG           >  1:19127598/1‑35
      GTAGCAGGCTGGTTTCAACCGTTGAATGGCAGGGG           >  1:11724850/1‑35
       TAGCAGGCTGGTTTCAACCGTTGAATGGCAGGGGG          >  1:15460511/1‑35
       TAGCAGGCTGGTTTCAACCGTTGAATGGCAGGGGG          >  1:8604599/1‑35
          CAGGCTGGTTTCAACCGTTGAATGGCAGGGGGGca       >  1:10844572/1‑33
          CAGGCTGGTTTCAACCGTTGAATGGCAGGGGGGGg       >  1:15180992/1‑34
           AGGCTGGTTTCAACCGTTGAATGGCAGGGGGGcag      >  1:7980870/1‑32
           AGGCTGGTTTCAACCGTTGAATGGCAGGAGGGGCG      >  1:18222591/1‑35
           AGGCTGGTTTCAACCGTTGAATGGCAGGAGGGGCG      >  1:16017432/1‑35
            GGCTGGTTTCAACCGTTGAATGGCAGGGGCGGCGG     >  1:5701936/1‑35
             GCTGGTTTCAACCGTTGAATGGCAGGGGGGGGGGG    >  1:14419950/1‑35
             GCTGGTTTCAACCGTTGAATGGCAGGGGGGGCGGG    >  1:6484225/1‑35
              CTGGTTTCAACCGTTGAATGGCAGGGGGGGCGGGG   >  1:13467863/1‑35
              CTGGTTTCAACCGTTGAATGGCAGGGGGGGCGGGG   >  1:9786142/1‑35
               TGGTTTCAACCGTTGAATGGCAGGGGGGGGGGGGg  >  1:9685832/1‑34
               TGGTTTCAACCGTTGAATGGCAGGGGGGGCGGGGg  >  1:3913297/1‑34
                                                  
CCGACAGTAGCAGGCTGGTTTCAACCGTTGAATGGCAG‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1661115‑1661078
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagGGGGGGCGGGGA  <  REL606/3625456‑3625445

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.