New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2589475	199 (0.910)	3 (0.010)	3/148	21.4	1.5%	coding (121/891 nt)	hcaR	DNA‑binding transcriptional activator of 3‑phenylpropionic acid catabolism
	?	REL606	4277048 =	190 (0.880)					coding (1141/2148 nt)	fdhF	formate dehydrogenase‑H, selenopolypeptide subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGAGAAAACATTCTCCCGCCGCCGTCAACGCGACTTTGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2589400‑2589475 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGT  >  REL606/4277048‑4277122                                                                                                                                                         GGAGAAAACATTCTCCCGCCGCCGTCAACGCGACTTTGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTT                                                                             >  1:5405202/1‑76                          CAACGCGACTTTGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTTAGCCGGTATGCGCTGGCAGGCTTTC                                                    >  1:3506430/1‑76                            ACGCGACTTTGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTTAGCCGGTATGCGCTGGCAGGCTTTCCA                                                  <  2:3506430/76‑1                                                                          GTTAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCAC    <  1:5197346/76‑1                                                                            TAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGT  <  2:5661716/76‑1                                                                            TAGCCGGTATGCGCTGGCAGGCTTTAAACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGCATCTTTCACGc  >  2:1761117/1‑75                                                                                                                                                         GGAGAAAACATTCTCCCGCCGCCGTCAACGCGACTTTGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2589400‑2589475 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGT  >  REL606/4277048‑4277122

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