New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 2018583	11 (0.230)	52 (1.120)	7/106	0.5	83.1%	coding (579/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	= 2521059	NA (NA)					intergenic (+170/+30)	insL/yfeA	predicted transposase/predicted diguanylate cyclase

TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGA  >  W3110S.gb/2521028‑2521113                                |                                                       tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:1406493/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:2023757/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:2663729/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:3819601/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:4100238/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:4142797/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:469623/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:487236/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:739596/1‑61 (MQ=255) tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:817178/1‑61 (MQ=255) tcAGATGGTTCAAACCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                           >  1:1504994/1‑61 (MQ=255)     aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                       <  1:2648150/61‑1 (MQ=255)     aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                       <  1:3150665/61‑1 (MQ=255)      tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                     >  1:1651708/1‑62 (MQ=255)              cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATg             >  1:4214549/1‑62 (MQ=255)                          cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGa  <  1:1039082/61‑1 (MQ=255)                                |                                                       TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGA  >  W3110S.gb/2521028‑2521113

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