New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	286345 =	NA (NA)	12 (0.980)	7/166	0.5	100%	intergenic (+2/+15)	BW25113_RS23270/BW25113_RS01405	YfjI family protein/IS1‑like element IS1B family transposase
	?	NZ_CP009273	3512441 =	0 (0.000)					intergenic (‑18/+383)	aroK/hofQ	shikimate kinase AroK/DNA uptake porin HofQ

GAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/286381‑286345 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCG  >  NZ_CP009273/3512441‑3512517                                                                                                                    GAACATAAAACACTATCAATAAGTTGGAGTCATTACCG                                                                              >  1:367608/1‑38 GAACATAAAACACTATCAATAAGTTGGAGTCATTACCG                                                                              <  2:367608/38‑1      TAAAACACTATCAATAAGTTGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGA                     >  2:307560/1‑90         AACACTATCAATAAGTTGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAC                  <  1:307560/90‑1         AACACTATCAATAAGTTGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAC                  <  1:76159/90‑1            ACTATCAATAAGTTGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTC               >  2:359454/1‑90                TCAATAAGTTGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAA           <  1:150923/90‑1                         TGGAGTCATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCG  >  1:113215/1‑90                               CATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATG                           <  1:267785/59‑1                               CATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATG                           <  1:93281/59‑1                               CATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATG                           >  2:267785/1‑59                               CATTACCGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATG                           >  2:93281/1‑59                                                                                                                    GAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/286381‑286345 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCG  >  NZ_CP009273/3512441‑3512517

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