New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	360104 =	NA (NA)	33 (3.240)	23/168	0.0	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

AGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTT  >  NZ_CP009273/360199‑360373                                                                                        |                                                                                        aGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                                                                       <  1:54940/90‑1 (MQ=40)   gagaGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                                                                                     <  2:242979/90‑1 (MQ=255)      agCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaatt                                                                                  <  1:94177/90‑1 (MQ=255)       gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                 >  2:80268/1‑90 (MQ=255)       gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                 >  2:193792/1‑90 (MQ=255)       gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                 >  1:242664/1‑90 (MQ=255)       gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                 >  1:180119/1‑90 (MQ=255)                                                 ttCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTtat                                       >  1:46306/1‑90 (MQ=255)                                                 ttCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTtat                                       >  2:251125/1‑90 (MQ=255)                                                       tGAGCCTTTCGTTTTATTTGATGCCCGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTa                                 >  1:71236/1‑90 (MQ=255)                                                       tGAGCCATTCGCTTTGTTTGATGCCCGGGGAATTGTGTATTAATCCTTGAAGACGAAAGGGCCTAGTGATAGGGCGATTTTTATAGGtga                                 >  2:35501/1‑88 (MQ=255)                                                        gAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTaa                                <  2:161914/90‑1 (MQ=255)                                                        gAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTaa                                <  2:286677/90‑1 (MQ=255)                                                                       ccagATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATc                 >  1:276258/4‑90 (MQ=255)                                                                                    gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc    <  1:251125/90‑1 (MQ=255)                                                                                    gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc    <  1:193792/90‑1 (MQ=255)                                                                                    gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc    <  1:106793/90‑1 (MQ=255)                                                                                      attgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCtt  >  2:219439/1‑90 (MQ=255)                                                                                        |                                                                                        AGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTT  >  NZ_CP009273/360199‑360373

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