New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	360104 =	NA (NA)	16 (0.950)	10/168	0.4	44.4%	intergenic (‑35/+251)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
	?	NZ_CP009273	= 360286	20 (1.190)					intergenic (‑217/+69)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360016‑360193                                                                                         |                                                                                          aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCaa                                                                                          <  1:326421/90‑1 (MQ=255)        gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTgg                                                                                   <  1:345033/90‑1 (MQ=255)                                     ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAc                                                      <  1:357628/90‑1 (MQ=255)                                             ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCcgg                                              >  1:73481/1‑90 (MQ=255)                                             ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGGccccccccccccccc                                              >  1:238111/1‑87 (MQ=255)                                                                gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGGCTGCCCCCCCCCCCCCCGGGCGGTGCTTCGTAAccgtt                           >  2:462055/1‑86 (MQ=19)                                                                                    tCCCCAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  2:268365/90‑1 (MQ=255)                                                                                     ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGcccc                                                      >  2:14504/1‑40 (MQ=255)                                                                                     ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAc                                                      <  1:14504/42‑1 (MQ=255)                                                                                     ccccAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTTCTTCGCAAAGTTTAAAACCCCCCCccgcggggtt      >  2:394235/1‑86 (MQ=9)                                                                                       ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt    <  1:33591/90‑1 (MQ=255)                                                                                         aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:462055/90‑1 (MQ=25)                                                                                         aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:368911/90‑1 (MQ=35)                                                                                         aaTTCCTGGCAGTTTATGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:427499/90‑1 (MQ=25)                                                                                         aaTCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:309320/90‑1 (MQ=9)                                                                                         |                                                                                          ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360016‑360193

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