Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	5,225,843	164 bp→TTGT	5.6%	coding (81‑244/876 nt)	PP_4605 →	Transcriptional regulator, AraC family

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5225842	72 (1.160)	4 (0.070) +TTGT	3/336	NT	5.6%	coding (80/876 nt)	PP_4605	Transcriptional regulator, AraC family
	?	AE015451	5226007 =	67 (1.080)					coding (245/876 nt)	PP_4605	Transcriptional regulator, AraC family

AGGGCGCAGTAGGATGGCAGGCAAGAGGACCCGCCCATGACGCCACTCACCCAGTTGCAAGTGTTCAACGCCATGCATGCCTCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5225727‑5225842 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGGCACACCTGTATGTCAGCGAGGCACAGTTCGCCCTGGGCTGTGTACGCCTGCTCGACCGGGAACCGCGCGAGCTGCAACTGCAGGAAGCGACCTTTCTGGACGACCCGCAGCAGGCTGGCCGCTTTCGT  >  AE015451/5226007‑5226182                                                                                                                     ||||                                                                                                                                                                                 AGGGCGCAGTAGGATGGCAGGCAAGAGGACCCGCCCATGACGCCACTCACCCAGTTGCAAGTGTTCAACGCCATGCATGCCTCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCTTTGTCGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGGCACACCTGTATGTCAGCGAGGCACAGTTCGCCCTGGGCTGTGTACGCCTGCTC                                                                             >  1:938180/1‑221                                                                             ATGCCTCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCTTTGTCGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGG                                                                                                                                  <  1:679391/92‑1                                                                             ATGCCTCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCTTTGTCGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGG                                                                                                                                  >  2:679390/1‑92                                                                                  TCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCTT                                                                                                                                                                                     >  2:3863/1‑36                                                                                                      ttgtcggGCGCACCTTTGTCGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGGCACACCTGTATGTCAGCGAGGCACAGTTCGCCCTGGGCTGTGTACGCCTGCTCGACCGGGAACCGCGCGAGCTGCAACTGCAGGAAGCGACCTTTCTGGACGACCCGCAGCAGGCTGGCCGCTTTCGT  >  1:874356/8‑195                                                                                                                     ||||                                                                                                                                                                                 AGGGCGCAGTAGGATGGCAGGCAAGAGGACCCGCCCATGACGCCACTCACCCAGTTGCAAGTGTTCAACGCCATGCATGCCTCGCCCAATGCCAGGCTGGAGCTGAGCGCGCACCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5225727‑5225842 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGGCCGGGCAGGAATCGAACTGGGTCGTCAATGGCGCGATCCGCCTGGCACACCTGTATGTCAGCGAGGCACAGTTCGCCCTGGGCTGTGTACGCCTGCTCGACCGGGAACCGCGCGAGCTGCAACTGCAGGAAGCGACCTTTCTGGACGACCCGCAGCAGGCTGGCCGCTTTCGT  >  AE015451/5226007‑5226182

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