New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	730478 =	11 (0.210)	3 (0.060)	3/146	3.9	23.7%	pseudogene (875/1506 nt)	BW25113_RS03635	RHS repeat protein
	?	NZ_CP009273	3613677 =	NA (NA)					coding (1126/4236 nt)	rhsB	rhs element protein RhsB
Rejected: Coverage evenness skew score above cutoff.

GAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATA  >  NZ_CP009273/730390‑730561                                                                                         |                                                                                    gAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGc                                                                                    >  2:1035832/1‑90 (MQ=35)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:1378174/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:1603246/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:229846/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:454173/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:775453/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  2:78283/90‑1 (MQ=31)                                                                             agcagcTTATCCGCCTGCCGGGGCAGCAGCATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGa        <  1:1463497/90‑1 (MQ=14)                                                                               cagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAgg      <  2:1048196/90‑1 (MQ=33)                                                                                   ttATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGtata  >  1:1214407/1‑90 (MQ=33)                                                                                   ttATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGgata  >  1:706949/1‑90 (MQ=18)                                                                                         |                                                                                    GAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATA  >  NZ_CP009273/730390‑730561

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