New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	131895 =	39 (0.980)	4 (0.110)	3/320	NT	10.2%	coding (87/618 nt)	cc	Cytochrome c4
	?	AE015451	131994 =	33 (0.870)					coding (186/618 nt)	cc	Cytochrome c4

CTGGTCGCTGAATGCGGCCAGCATGCCGGTCATTTCCAGCACGGTACGCTTGCCCGACTTGATATCGTGAAGTTGTTTTTCGAGATAGCGCTGGCCTTGGCCGGCCAGTTTCGGGAAGTTCGGTGCCAGGCTGTTGCCGTCGGGATTGTGGCAGGCGCCACACACGGCGGTCTTGGCCTGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/132081‑131895 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccagcgCTATCTCGAAAAACAACTTCACGATATCAAGTCGGGCAAGCG  >  AE015451/131994‑132035                                                                                                                                                                                                                                       CTGGTCGCTGAATGCGGCCAGCATGCCGGTCATTTCCAGCACGGTACGCTTGCCCGACTTGATATCGTGAAGTTGTTTTTCGAGATAGCGCTGGCCTTGGCCGGCCAGTTTCGGGAAGTTCGGTGCCAGGCTGTTGCCGTCGGGATTGTGGCAGGCGCCACACACGGCGGTCTTGGCCTGGCCAGCGCTATCTCGAAAAACAACTTCACGATATCA               <  2:553630/216‑1                                                              ATATCGTGAGGTTGTTTTTCGAGATAGCGCTGGCCTTGGCCGGCCAGTTTCGGGAAGTTCGGTGCCAGGCTGTTGCCGTCGGGATTGTGGCAGGCGCCACACACGGCGGTCTTGGCCTGGCCAGCGCTATCTCGAAAAACAACTTCACGATATCAAG             <  2:520808/157‑1                                                              ATATCGTGAAGTTGGTTTTCGAGATAGCGCTGGCCTTGGCCGGCCAGTTTCGGGAAGTTCGGTGCCAGGCTGTTGCCGTCGGGATTGTGGCAGGCGCCACACACGGCGGTCTTGGCCGGGCCAGCGCTATCTCGAAAAACAACTTCACGATATCAAGTCGGGCAAGCG  <  1:646407/168‑1                                                                                                                                             CGGGATTGTGGCAGGCGCGACACACGGCGGTCTTGGCCTGGCCAGCGCTATCTCGAAAAACAACTTCACGATATCA               <  2:441990/76‑1                                                                                                                                                                                                                                       CTGGTCGCTGAATGCGGCCAGCATGCCGGTCATTTCCAGCACGGTACGCTTGCCCGACTTGATATCGTGAAGTTGTTTTTCGAGATAGCGCTGGCCTTGGCCGGCCAGTTTCGGGAAGTTCGGTGCCAGGCTGTTGCCGTCGGGATTGTGGCAGGCGCCACACACGGCGGTCTTGGCCTGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/132081‑131895 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccagcgCTATCTCGAAAAACAACTTCACGATATCAAGTCGGGCAAGCG  >  AE015451/131994‑132035

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