New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5821896	41 (0.940)	3 (0.070)	3/324	NT	7.1%	coding (517/639 nt)	metW	methionine biosynthesis protein MetW
	?	AE015451	= 5821897	40 (0.970)					coding (518/639 nt)	metW	methionine biosynthesis protein MetW

GGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5821658‑5821896 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgaccgccAGGCGGTCGAGTACCTTGGCCCGGCGTTCGTGCACAAGTTCTTCGAAGTCGGCGAAGGTGCAGAAGTGGATGTTCGGCGTGTTGTACCAGGTATACGGCATGAAGTCCGATACCGGCATGCGGCCTTTGGTCGCCAGGTACCAGC  <  AE015451/5821897‑5821753                                                                                                                                                                                                                                                                                                                                                                                                  GGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGTTACTCGAC                                                                                                                                                       >  1:133202/1‑235                                                                             TTCGGCCACTGGCGCTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGACGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCC                                                                                                                                                   <  2:430069/163‑1                                                                                                                                              GCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGAGTACCTTGGCCCGGCGTTCGTGCACAAGTTCTTCGAAGTCGGCGAAGGTGCAGAAGTGGATGTTCGGCGTGTTGTACCAGGTATACGGCATGAAGTCCGATACCGGCGTGCGGCCTTTGGTCGCCAGG         >  1:527431/1‑236                                                                                                                                                     ACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGAGTACCTTGGCCCGGCGTTCGTGCACAAGTTCTTCGAAGTCGGCGAAGGTGCAGAAGTGGATGTTCGGCGTGTTGTACCAGGTATACGGCATGAAGTCCGATACCGGCATGCGGCCTTTGGTCGCCAGGTACCAGC  >  1:685100/1‑236                                                                                                                                                                                        CTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGAGTACCTTGGCCCG                                                                                                                            >  2:803712/1‑79                                                                                                                                                                                                                                                                                                                                                                                                  GGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5821658‑5821896 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgaccgccAGGCGGTCGAGTACCTTGGCCCGGCGTTCGTGCACAAGTTCTTCGAAGTCGGCGAAGGTGCAGAAGTGGATGTTCGGCGTGTTGTACCAGGTATACGGCATGAAGTCCGATACCGGCATGCGGCCTTTGGTCGCCAGGTACCAGC  <  AE015451/5821897‑5821753

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