New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	minE	19620 =	64 (1.040)	7 (0.130) +GTGACT	3/86	NT	11.6%	coding (536/942 nt)	ribF	bifunctional riboflavin kinase and FAD synthetase
	?	minE	= 2548513	58 (0.940)					coding (722/858 nt)	htrL	hypothetical protein
Rejected: Coverage evenness skew score above cutoff.

CCAGTAAACTCTCTGCCAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/19640‑19620 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CACATGACAAAAA  <  minE/2548513‑2548501                      ||||||              tgAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  <  1:3852700/38‑1  gAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  >  1:2325791/2‑39  gAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  >  1:3497478/2‑39   AGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  >  1:1293542/1‑38   AGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  <  1:1612156/38‑1   AGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  <  1:2651763/38‑1   AGTAAACTCTCAGTCAGAGGTGACTCACATGACAAA    <  1:5357655/36‑1                      ||||||              CCAGTAAACTCTCTGCCAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/19640‑19620 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CACATGACAAAAA  <  minE/2548513‑2548501

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