New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 3538215	38 (1.470)	3 (0.120)	3/302	NT	7.8%	coding (37/1116 nt)	PP_5587	conserved protein of unknown function
	?	AE015451	= 3538223	35 (1.430)					coding (45/1116 nt)	PP_5587	conserved protein of unknown function

GCTGGCGATACGCTGGTCGTTCCTTGATCTTGCCTTTAATAAGAACATTTAATGCCTCTTCTTAGAACAACTCGCATTGCAACGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/3538128‑3538215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aacggcaaGACCAAGATTGCCGTTGCAATGCGAGTTGTTCTAAGAAGAGGCATTAAATGTTCTTATTAAAGGCAAGATCAAGGAACGACCAGCGTATCGCCAGCCTGAAGTTG  <  AE015451/3538223‑3538119                                                                                                                                                                                                   GCTGGCGATACGCTGGTCGTTCCTTGATCTTGCCTTTAATAAGAACATTTAATGCCTCTTCTTAGAACAACTCGCATTGCAACGGCAAGACCAAGCTTGCCGTTGCAATGCGAGTTGTTCTCAGAAGAGGCATTAAATGTTCTTATTAAAGGCAAGATCAAGGAACGACCAGCGTAac                 >  1:373956/1‑176    aGCGATACGCTGGTCGTTCCTTGATCTTGCCTTTAATAAGAACATTTAATGCCTCTTCTTAGAACAACTCGCATTGCAACGGC                                                                                                             <  2:322704/82‑1                         TGATCTTGCCTTTAATAAGAACATTTAATGGCTCTTCTTAGAACAACTCGCATTGCAACGGCAAGACCAAGATTGCCGTTGCAATGCGAGTTGTTCTAAGAAGAGGCATTAAATGTTCTTATTAAAGGCAAGATCAAGGAACGACCAGCGTATCGCCAGCCTGAAGGTG  <  2:470931/169‑1                         TGATCTTGCCTTTAATAAGAACATTTAATGCCTCTTCTTAGAACAACTCGCATTGCAACGGCAAGACCAAGATTGCCGTTGCAATGCGAGTTGTTCTAAGAAGAGGCATTAAATGTTCTTATTAAAGGCAAGATCAAGGAACGACCAGCGTATCGCCAGCCTGAAGTTG  >  1:470935/1‑169                                                                                                                                                                                                   GCTGGCGATACGCTGGTCGTTCCTTGATCTTGCCTTTAATAAGAACATTTAATGCCTCTTCTTAGAACAACTCGCATTGCAACGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/3538128‑3538215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aacggcaaGACCAAGATTGCCGTTGCAATGCGAGTTGTTCTAAGAAGAGGCATTAAATGTTCTTATTAAAGGCAAGATCAAGGAACGACCAGCGTATCGCCAGCCTGAAGTTG  <  AE015451/3538223‑3538119

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