New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	253338 =	36 (1.020)	5 (0.180) +CCATTTATTATCTTTCTGCGCACTTCACG	3/196	3.9	16.4%	intergenic (+177/‑141)	yafP/ykfJ	GNAT family putative N‑acetyltransferase/pseudogene
	?	NC_000913	= 803472	30 (0.850)					intergenic (+152/‑31)	ybhI/ybhJ	putative DASS family tricarboxylate or dicarboxylate transporter/aconitase family protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/253380‑253338 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGCGCAGAAAGATAATAAATGGATATACTCTAAATAATTCGAGTTGCAGGAAGGCGACAAGCGAGTGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGG  <  NC_000913/803472‑803370                                            |||||||||||||||||||||||||||||                                                                                                        CTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGAAAGATAATAAATGGATATA                                                                             >  5:76655/1‑100 CTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGAAAGATAATAAATGGATATA                                                                             <  6:76655/100‑1                  TGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGAAAGATAATAAATGGATATACTCTAAATAATTCGAGTTGCAGGAAGGCGACAAGCGAGTGAATCCCCAGGAGCTTACATAAGTAAG           >  6:67710/1‑149                             TTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACCGTGAAGTGCGCAGAAAGATAATAAATGGATATACTCTAAATAATTCGAGTTGCAGGAAGGCGACAAGCGAGTGAATCCCCAGTCACTTACTTATGTAAGCTCCTGGGG  >  5:173040/1‑147                             TTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACCGTGAAGTGCGCAGAAAGATAATAAATGGATATACTCTAAATAATTCGAGTTGCAGGAAGGCGACAAGCGAGTGAATCCCCAGTCACTTACTTATGTAAGCTCCTGGGG  <  6:173040/147‑1                                            |||||||||||||||||||||||||||||                                                                                                        CTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/253380‑253338 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGCGCAGAAAGATAATAAATGGATATACTCTAAATAATTCGAGTTGCAGGAAGGCGACAAGCGAGTGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGG  <  NC_000913/803472‑803370

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