New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	360104 =	NA (NA)	16 (1.760)	11/166	0.1	42.1%	intergenic (‑35/+251)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
	?	NZ_CP009273	= 360286	22 (2.410)					intergenic (‑217/+69)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

TGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360022‑360193                                                                                   |                                                                                          tGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTg                                                                                    >  2:274017/1‑90 (MQ=255)   ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGc                                                                                  <  1:251565/90‑1 (MQ=255)                                       ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCcgg                                              >  1:67922/1‑90 (MQ=255)                                            tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATAGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                         >  1:93640/1‑90 (MQ=255)                                                                              tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  1:224976/90‑1 (MQ=255)                                                                                 ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt    <  1:26697/90‑1 (MQ=38)                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:120898/90‑1 (MQ=25)                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:274017/90‑1 (MQ=35)                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:61802/90‑1 (MQ=35)                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCCCCCGGCGGATTTGTcc  <  1:42169/90‑1 (MQ=25)                                                                                   |                                                                                          TGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360022‑360193

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