New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	3049681 =	5 (0.020)	55 (0.260)	34/242	10.8	37.3%	coding (841/1188 nt)	PP_2662	conserved exported protein of unknown function
	?	AE015451	3366553 =	180 (0.840)					intergenic (+219/‑142)	PP_2962/PP_2964	Alcohol dehydrogenase, zinc‑containing/putative transposase
Rejected: Coverage evenness skew score above cutoff.

GCAGATCTACAGTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTCGGCAA  >  AE015451/3049560‑3049805                                                                                                                          |                                                                                                                             gcAGATCTACAGTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGa                                                                                                                           <  1:2635126/125‑1 (MQ=255) gcAGATCTACAGTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGa                                                                                                                           <  1:2635127/125‑1 (MQ=255)      tCTACAGTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGACgggg                                                                                                                      <  1:2635128/125‑2 (MQ=255)            gTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTNCAA‑GGNGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGACGGGGTcatgc                                                                                                                <  1:2635129/124‑5 (MQ=255)                                                         cGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGa                                                                                                                           <  1:2635136/69‑1 (MQ=255)                                                                                                                    aacccccTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGtt        >  1:2635145/7‑125 (MQ=255)                                                                                                                    aacccccTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGtt        >  1:2635157/7‑125 (MQ=255)                                                                                                                     acccccTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTc       >  2:2635149/6‑125 (MQ=255)                                                                                                                       ccccTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTCgg     >  1:2635180/4‑125 (MQ=255)                                                                                                                        cccTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTCGGc    >  2:2635173/3‑125 (MQ=255)                                                                                                                          cTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTCGGCaa  >  2:2635156/1‑125 (MQ=255)                                                                                                                          |                                                                                                                             GCAGATCTACAGTTGGGTCAACGGCAGTTACCAGACCTCGGACAATACTGACTCCACGGTAGAAACGATCACTTCCAAGGGGGTGGGTTACGGGGTGCAGGCAAAGATGGGCGGCTGGTCGCTGACCGGCTCGGGGTTCCAGGCCAAAGGCATCAACCCGTTCTTCACCAACAATGCCGGCGAACCGGTTTTGCGCAATGTCGACAGTGATGGCTACCTGCTGCAGGGCTCGTACAAGTTCGGCAA  >  AE015451/3049560‑3049805

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