New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 6002217	40 (0.980)	3 (0.070)	3/330	NT	6.8%	coding (84/1299 nt)	amaA	L‑pipecolate oxidase
	?	AE015451	= 6002230	43 (1.070)					coding (71/1299 nt)	amaA	L‑pipecolate oxidase

GGGCGGATCCACGTGCCGGCGTTGACCAGGCCGACGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/6002044‑6002217 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgCCCAGGCCCTGGCCGGTGAGCACAAGGCCGATGTCTGTGTGATCGGCGGCGGTATTACCGGCCTGTCGGCGGCCATCCACCTGCTCGAACAGGGCAAGTCGGTGATCGTGCTGGAGGCCTGGAAGATCGGCCA  <  AE015451/6002230‑6002098                                                                                                                                                                                                                                                                                                                     GGGCGGATCCACGGGCCGGCGTGGACCAGGCCGACGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACC                                                                                                                                        <  2:633222/173‑1                                                                                                                          ACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGAGCACAAGGCCGATGTCTGTGTGATCGGCGGCGGTATTACCGGCCTGTCGGCGGCCATCCACCTGCTCGAACAGCGCAGGTCGGTGATCGTGCTGGAGGCCTGGc            >  1:403235/1‑175                                                                                                                                             GCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGAGCACAAGGCCGATGTCTGTGTGATCGGCGGCGGTATTACCGGCCTGTCGGCGGCCATCCACCTGCTCGAACAGGGCAAGTCGGTGATCGTGCTGGAGGCCTGGAAGATCGGCCA  >  1:164014/1‑167                                                                                                                                              CCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGAGCACAAGGCCGAT                                                                                                       >  2:319635/1‑65                                                                                                                                                                                                                                                                                                                     GGGCGGATCCACGTGCCGGCGTTGACCAGGCCGACGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/6002044‑6002217 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgCCCAGGCCCTGGCCGGTGAGCACAAGGCCGATGTCTGTGTGATCGGCGGCGGTATTACCGGCCTGTCGGCGGCCATCCACCTGCTCGAACAGGGCAAGTCGGTGATCGTGCTGGAGGCCTGGAAGATCGGCCA  <  AE015451/6002230‑6002098

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