New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3512449	0 (0.000)	10 (0.810)	10/166	0.3	100%	intergenic (‑26/+375)	aroK/hofQ	shikimate kinase AroK/DNA uptake porin HofQ
	?	NZ_CP009273	3576788 =	NA (NA)					coding (391/417 nt)	BW25113_RS17825	YrhA family protein

AAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCT  >  NZ_CP009273/3576718‑3576864                                                                       |                                                                             aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACt                                                           <  1:177980/90‑1 (MQ=255)                                ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGGTTTAGTGTATGATGGTGTTTTTGAGGGGc                            >  2:295614/1‑90 (MQ=255)                                            aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                >  2:146487/1‑90 (MQ=255)                                                       tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa     >  2:125257/1‑90 (MQ=255)                                                            aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt  >  1:178002/1‑88 (MQ=255)                                                            aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt  <  2:178002/88‑1 (MQ=255)                                                                       |                                                                             AAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCT  >  NZ_CP009273/3576718‑3576864

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