New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	360104 =	NA (NA)	17 (0.940)	16/170	0.2	NA	intergenic (‑35/+251)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
	?	NZ_CP009273	= 360286	NA (NA)					intergenic (‑217/+69)	lacZ/BW25113_RS23320	beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGC  >  NZ_CP009273/360197‑360371                                                                                          |                                                                                      tCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                       <  1:406133/90‑1 (MQ=38) tCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                       >  2:390758/1‑90 (MQ=38)         gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                               >  1:354029/1‑90 (MQ=255)         gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                               >  1:35891/1‑90 (MQ=255)                          acaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                                                              <  1:442160/90‑1 (MQ=255)                             gATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc                                                           >  2:163832/1‑90 (MQ=255)                                      aaaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa                                                  <  1:335899/90‑1 (MQ=255)                                      aaaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa                                                  <  1:450568/90‑1 (MQ=255)                                       aaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                       >  1:274973/1‑52 (MQ=255)                                       aaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                       <  2:274973/52‑1 (MQ=255)                                                         tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGACTTAATTCTTGAAGACGAAAGGGCCTCGTGAt                                                   <  1:385775/70‑1 (MQ=255)                                                         tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt                                                   >  2:385775/1‑70 (MQ=255)                                                         tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTa                               >  1:28087/1‑90 (MQ=255)                                                               tttCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCa                         <  1:196755/90‑1 (MQ=255)                                                                    ttttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGa                    <  1:163832/90‑1 (MQ=255)                                                                                      gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGCGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc  <  2:69978/90‑1 (MQ=255)                                                                                          |                                                                                      TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGC  >  NZ_CP009273/360197‑360371

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