New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4568611 =	40 (0.800)	4 (0.080)	4/322	NT	9.4%	coding (541/2070 nt)	malQ	4‑alpha‑glucanotransferase
	?	AE015451	4568628 =	39 (0.820)					coding (558/2070 nt)	malQ	4‑alpha‑glucanotransferase

TGAGCAGCAGGCGACTGGAAGGTGAGTAGGGGCTGTAGTGTTCCTGGTCGATGGCCGACAGGGCGTGAATCGGGCTGATAGCCAAGGCATCGGCGCCGCGCTCGGCGGCCGAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4568725‑4568611 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgagcgCGGCGCCGATGCCTTGGCTATCAGCCCGATTCACGCCCTGTCGGCCATCGACCAGGAACACTACAGCC  >  AE015451/4568628‑4568695                                                                                                                                                                                         TGAGCAGCAGGCGACTGGAAGGTGAGTAGGGGCTGTAGTGTTCCTGGTCGATGGCCGACAGGGCGTGAATCGGGCTGATAGCCAAGGCATCGGCGCCGCGCTCGGCGGCCGAGCGCGGCGCCGATGCCTTGGCTATCAGCCCGATTCACGCCCTGTCGGCCATCGACCAGGAACACTACAGCC  >  1:914979/1‑183       GCAGGCGACTGGAAGGTGAGTAGGGGCTGTAGTGTTCCTGGTCGATGGCCGACAGGGCGTGAATCGGGCTGATAGCCAAGGCATCGGCGCCGCGCTCGGCGGCCGAGCGCGGCGCCGATGCCTTGGCTATCAGCCCGATTCACGCCCTGTCGGCCATCGACCAG               <  2:914976/164‑1                                                       CCGACAGGGCGTGAATCGGGCTGATAGCCAAGGCATCGGCGCCGCGCTCGGCGGCCGAGCGCGGCGCCGATGCCTTGGCTATCAGCCCGATTCACGCCCTGTCGGCCATCGACCAGGAAC           <  1:358039/120‑1                                                                                                     CTCGGCGGCCGAGCGCGGCGCCGATGCCTTGGCTATCA                                               <  1:536559/38‑1                                                                                                                                                                                         TGAGCAGCAGGCGACTGGAAGGTGAGTAGGGGCTGTAGTGTTCCTGGTCGATGGCCGACAGGGCGTGAATCGGGCTGATAGCCAAGGCATCGGCGCCGCGCTCGGCGGCCGAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4568725‑4568611 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgagcgCGGCGCCGATGCCTTGGCTATCAGCCCGATTCACGCCCTGTCGGCCATCGACCAGGAACACTACAGCC  >  AE015451/4568628‑4568695

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