New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5200909 =	32 (0.730)	3 (0.070)	3/332	NT	8.9%	coding (94/774 nt)	yafJ	putative amidotransferase
	?	AE015451	5200932 =	30 (0.710)					coding (71/774 nt)	yafJ	putative amidotransferase

ACCGTGACGACTTGCCGCCAGCCGCTGCGCGTGGCGCAGGCCGGGTCGAGCCCCGTCTGTAAACGGAGACACGCCTATGTGTGAACTGCTGGGCATGAGTGCCAACGTCCCTACCGACATCGTTTTCAGCTTCACCGGCCTGATGCAGCGGGGCGGGCGAACCGGC‑‑‑‑‑CCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5201078‑5200909 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑c‑‑‑‑‑ccgcTGCATCAGGCCGGTGAAGCTGAAAACGATGTCGGTAGGGACGTTGGCACTCATGC  >  AE015451/5200932‑5200986                                                                                                                                                                                                                                        ACCGTGACGACGTGCCGCCAGCCGCTGCGCGTGGCGCAGGCCGGGTCGAGCCCCGTCTGTAAACGGAGACACGCCTATGTGTGAACTGCTGGGCATGAGTGCCAACGTCCCTACCGACATCGTTTTGAGCTTCACCGTCCTGATGCAGCGGGGCGGGCGAACCGGC‑‑‑‑‑CCGCTGCATCAGGCCGGTGAAGCTGAAAACGATGTCGGTAGGGACGTTGGCACTCATGC  <  1:435691/225‑1                                    GCAGGCCGGGTCGAGCCCCGTCTGTAAACGGAGACACGCCTATGTGTGAACTGCTGGGCATGAGTGCCAACGTCCCTACCGACATCGTTTTCAGCTTCACCGGCCTGATGCAGCGGGGCGGGCGAACCGGCCCGCACCGCTGCATCAGGCCGGTGAAGCTGAAAACGATGTCGGT                      <  2:777425/175‑1                                                   CCCCGTCTGTAAACGGAGACACGCCTACGTGTGAACTGCTGGGCATGAGTGCCAACGTCCCTACCGACATCGTTTTCAGCTTCACCGGCCTGATGCAGCGGGGCGGGCGAACCGGC‑‑‑‑‑CCGCTGCATCAGGCCGGTGAAGCTGAAAACGATGTCGGTAGGGACGTTGG           <  1:572882/166‑1                                                                                                                                                                                                                                        ACCGTGACGACTTGCCGCCAGCCGCTGCGCGTGGCGCAGGCCGGGTCGAGCCCCGTCTGTAAACGGAGACACGCCTATGTGTGAACTGCTGGGCATGAGTGCCAACGTCCCTACCGACATCGTTTTCAGCTTCACCGGCCTGATGCAGCGGGGCGGGCGAACCGGC‑‑‑‑‑CCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5201078‑5200909 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑c‑‑‑‑‑ccgcTGCATCAGGCCGGTGAAGCTGAAAACGATGTCGGTAGGGACGTTGGCACTCATGC  >  AE015451/5200932‑5200986

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