New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3652533	NA (NA)	17 (0.260)	13/200	4.5	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	= 3697156	0 (0.000)					coding (280/4134 nt)	rhsA	rhsA element core protein RshA
Rejected: Coverage evenness skew score above cutoff.

TGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAA  >  REL606/3697056‑3697166                                                                                                     |           tGAACCCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            <  2:354767/101‑1 (MQ=255) tGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            >  1:656294/1‑101 (MQ=40) tGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            >  2:720422/1‑101 (MQ=40)  gAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcct           >  1:1304378/1‑100 (MQ=255)  gAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcct           <  2:1581706/101‑2 (MQ=40)      ccGTCCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  1:515319/101‑6 (MQ=255)      ccGACCTTGCCCTGCCCGGCCCGCTGCCGTTCGTCCTCTCCCGCTCCGACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  2:934216/101‑6 (MQ=14)      ccGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  1:106502/101‑6 (MQ=255)       cGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtact      <  1:856289/101‑7 (MQ=255)        gACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactg     >  1:1466097/1‑94 (MQ=255)         aCCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactga    <  1:247422/101‑9 (MQ=14)          ccTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgac   >  1:834712/1‑92 (MQ=255)           cTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgacc  >  1:168788/1‑91 (MQ=255)           cTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgacc  >  2:638848/1‑91 (MQ=255)                                                                                                     |           TGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAA  >  REL606/3697056‑3697166

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