New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3248600 =	22 (0.730)	4 (0.150)	4/148	1.7	19.3%	intergenic (+34/+100)	BW25113_RS16135/cyuA	hypothetical protein/cysteine desulfidase
	?	NZ_CP009273	3248653 =	14 (0.530)					intergenic (+87/+47)	BW25113_RS16135/cyuA	hypothetical protein/cysteine desulfidase
Rejected: Frequency below/above cutoff threshold.

GGTTTTGGGGATCGTGCCGGATGCGATGCTGGCGCATCTTGTCCGGCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3248693‑3248600 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtaggccggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  NZ_CP009273/3248653‑3248739                                                                                                                                                                                       GGTTTTGGGGATCGTGCCGGATGCGATGCTGGCGCATCTTGTCCGGCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGC                                                                                             <  2:240502/90‑1           ATCGTGCCGGATGCGATGCTGGCGCATCTTGTCCGGCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGCCGGACAAGAT                                                                                   <  2:531031/90‑1                                    ATCTTGTCCGGCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGAT                                                          <  1:716748/90‑1                                              GCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACC                                                <  2:645363/90‑1                                                             AATTCGTTAATGTTGTTGGATTTTGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACTATCCCCAAAACCTGGCGTTTATCTGGC                                 <  2:581561/90‑1                                                                                      TAGGCCAGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATC        >  2:443437/1‑90                                                                                            GGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  1:536738/1‑90                                                                                            GGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  1:86569/1‑90                                                                                            GGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  2:253702/1‑90                                                                                                                                                                                       GGTTTTGGGGATCGTGCCGGATGCGATGCTGGCGCATCTTGTCCGGCCTACGGGTTGTTTAATTCGTTAATGTTGTTGGATTTTGTAGGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3248693‑3248600 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtaggccggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  NZ_CP009273/3248653‑3248739

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