New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 2412887	30 (0.750)	3 (0.080)	3/316	NT	10.7%	intergenic (+151/‑270)	rlmM/PP_2114	Ribosomal RNA large subunit methyltransferase M/transposase
	?	AE015451	= 4908014	22 (0.590)					coding (11/198 nt)	PP_4317	conserved protein of unknown function

GAAGTGACGTGCCACCTGCGGCGGCTCGATTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/2412692‑2412887 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgcccgtcaCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAA  <  AE015451/4908014‑4907930                                                                                                                                                                                                                                                                                           gacGTGACGTGCCACCTGCGGCGGCTCGATTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTC                                                 <  2:640945/231‑1                         CTCGATTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTCGC                                               <  2:693381/212‑1                                                                                                                  ATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGGAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCACCTGCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAA  <  1:25456/168‑1                                                                                                                                                                                                                                                                                           GAAGTGACGTGCCACCTGCGGCGGCTCGATTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/2412692‑2412887 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgcccgtcaCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAA  <  AE015451/4908014‑4907930

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