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

TCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360065‑360193                                        |                                                                                          tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcggg                                                                     >  1:72043/1‑62 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcggg                                                                     <  2:72043/62‑1 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATg                                                                          >  1:244464/1‑57 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATg                                                                          <  2:244464/57‑1 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGcccct                                         >  2:363403/1‑87 (MQ=37)                              tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc            <  1:289668/90‑1 (MQ=255)                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAAACCCCCCCccgcgggt        >  2:154871/1‑88 (MQ=14)                                   tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  1:314119/90‑1 (MQ=255)                                      ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt    <  2:179711/90‑1 (MQ=38)                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:32350/90‑1 (MQ=35)                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:152486/90‑1 (MQ=35)                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTCGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:47188/90‑1 (MQ=25)                                        |                                                                                          TCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360065‑360193

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