New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1612861	11 (0.430)	3 (0.120)	3/298	NT	19.8%	coding (658/693 nt)	ung	Uracil‑DNA glycosylase
	?	AE015451	= 1612883	14 (0.580)					coding (680/693 nt)	ung	Uracil‑DNA glycosylase

CCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1612687‑1612861 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgggggcAGTGCCCAGTCGATGGGCCCGAGCCCCCGCTGCTCGAGGAAGCTGTTGGCACGGCTGAAATGCCCGCAACCGATGAACCCGCGATACGCCGACAGCGGCGAAGGATGGACCGACTTCAAAACCAGATGCTTGGTACCGTCGATCAGCTTCTGCTTGCTCTGTGCATGCGCCCCCCATAGCAGGAACACCACGTTCGGGCATTGCTCACTGACCACCTGGATAATCCTGTCGG  <  AE015451/1612883‑1612652                                                                                                                                                                                                                                                                                                                                                                                                                         CCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGGTTTGAAGTCGGTCAATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCCTCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCAGTGCCCAGTCGATGGGCCC                                                                                                                                                                                                                      >  1:158141/1‑195                                                                              TTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCCTTTCAGCCGTTCCACCAGCTTCCTCGCGCAGC                                                                                                                                                                                                                                                >  2:35451/1‑92                                                                                                   GCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCAGTGCCCAGTCGATGGGCCCGAGCCCCCGCT                                                                                                                                                                                                           <  2:158137/108‑1                                                                                                                                                        ACAGCTTCCTCGAGCAGCGGGGGCAGTGCCCAGTCGATGGGCCCGAGCCCCCGCTGCTCGA                                                                                                                                                                                                     >  2:202228/1‑61                                                                                                                                                                            GGGCAGTGCCCAGTCGATGGGCCCGAGCCCCCGCTGCTCGAGGAAGCTGTTGGCACGGCTGAAATGCCCGCAACCGATGAACCCGCGATACGCCGACAGCGGCGAAGGATGGACCGACTTCAAAACCAGATGCTTGGTACCGTCGATCAGCTTCTGCTTGCTCTGTGCATGCGCCCCCCATAGCAGGAACACCACGTTCGGGCATTGCTCACTGACCACCTGGATAATCCTGTCGG  <  1:16257/236‑1                                                                                                                                                                            GGGCAGTGCCCAGTCGATGGGCCCGAGCCCCCGCTGCTCGAGGAAGCTGTTGGCACGGCTGAAATGCCCGCAACCGATGAACCCGCGATACGCCGACAGCGGCGAAGGATGGACCGACTTCAAAACCAGATGCTTGGTACCGTCGATCAGCTTCTGCTTGCTCTGTGCATGCGCCCCCCATAGCA                                                     >  2:315983/1‑185                                                                                                                                                                                                                                                                                                                                                                                                                         CCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1612687‑1612861 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgggggcAGTGCCCAGTCGATGGGCCCGAGCCCCCGCTGCTCGAGGAAGCTGTTGGCACGGCTGAAATGCCCGCAACCGATGAACCCGCGATACGCCGACAGCGGCGAAGGATGGACCGACTTCAAAACCAGATGCTTGGTACCGTCGATCAGCTTCTGCTTGCTCTGTGCATGCGCCCCCCATAGCAGGAACACCACGTTCGGGCATTGCTCACTGACCACCTGGATAATCCTGTCGG  <  AE015451/1612883‑1612652

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