New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	minE	1258762 =	57 (1.100)	4 (0.090) +CAACGT	3/90	NT	8.0%	coding (2162/2430 nt)	torZ	trimethylamine N‑oxide reductase system III, catalytic subunit
	?	minE	= 1519712	47 (0.910)					coding (282/660 nt)	dedA	conserved inner membrane protein

CGAATGGCGATCTGGTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/1258780‑1258762 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTTGTTCGGTGAAAAG  <  minE/1519712‑1519697                    ||||||                 CGAATGGCGATCTGGTGCGCAACGTGTTGTTCGGT        >  1:1844455/1‑35     TGGCGATCTGGTGCGCAACGTGTTGTTCGGTGAAAA   <  1:3438834/36‑1     TGGCGATCTGGTGCGCAACGTGTTGTTCGGTGAAAA   >  1:3737543/1‑36       GCGATCTGGTGCGCAACGTGTTGTTCGGTGAAAAc  >  1:2647962/1‑34                    ||||||                 CGAATGGCGATCTGGTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/1258780‑1258762 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTTGTTCGGTGAAAAG  <  minE/1519712‑1519697

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