New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4584201 =	26 (0.640)	3 (0.080)	3/272	6.0	10.0%	coding (2561/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
	?	NC_000913	4584215 =	29 (0.750)					coding (2547/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTGTAATATCGTGTTCCTGCGTAAAGTACGCAGCCGCATTCTGTACGAACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4584347‑4584201 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgtGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAGGAACACGATATTACAGATCGAC  >  NC_000913/4584215‑4584353                                                                                                                                                                                                                                                                                                CTGTAATATCGTGTTCCTGCGTAAAGTACGCAGCCGCATTCTGTACGAACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGA                                                                                                                                                 >  1:421895/1‑143                                                                                                                                TACAGCACGCTGGAGAGCGTGCTGTAGATATCGACACAGTCAAGAAT                                                                                                                  >  1:612214/1‑47                                                                                                                                TACAGCACGCTGGAGAGCGTGCTGTAGATATCGACACAGTCAAGAAT                                                                                                                  <  2:612214/47‑1                                                                                                                                       ctgtcgcGAGCGTGCTGTAGCTATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAGGAACACGATATTA           <  1:645553/136‑1                                                                                                                                                GCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTTGTACAGAATGCGGCTGCGTACTTTACGCAGGAACACGATATTACAGATCGAC  <  1:419336/143‑1                                                                                                                                                                                                                                                                                                CTGTAATATCGTGTTCCTGCGTAAAGTACGCAGCCGCATTCTGTACGAACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4584347‑4584201 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgtGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAGGAACACGATATTACAGATCGAC  >  NC_000913/4584215‑4584353

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