New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 4571590	79 (1.260)	3 (0.050)	3/158	3.9	3.7%	coding (1276/1590 nt)	hsdM	type I restriction‑modification system methyltransferase
	?	NZ_CP009273	4571594 =	81 (1.330)					coding (1272/1590 nt)	hsdM	type I restriction‑modification system methyltransferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑TC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4571541‑4571590 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGT  >  NZ_CP009273/4571594‑4571682                                                                                                                                              AGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAA                                                    >  1:855057/1‑90             GTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑TCGGCGTTAAAACTCC                                                                             >  1:1454874/1‑52             GTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑TCGGCGTTAAAACTCC                                                                             <  2:1454875/52‑1                                                   CGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGT  <  2:836938/90‑1                                                                                                                                              AGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑TC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4571541‑4571590 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGT  >  NZ_CP009273/4571594‑4571682

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