Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	5,859,118	40 bp→GCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGCCG	7.4%	coding (829‑868/1026 nt)	PP_5135 ←	ABC transporter, periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5859117	62 (1.080)	4 (0.090) +GCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGCCG	4/262	NT	7.4%	coding (869/1026 nt)	PP_5135	ABC transporter, periplasmic binding protein
	?	AE015451	5859158 =	68 (1.180)					coding (828/1026 nt)	PP_5135	ABC transporter, periplasmic binding protein

GCCATTCGGCCAGAATCTGGTCACGGTTTTTAGAAGCCCAGACAAAGTCGTTCTTGATCAACCGCTGCTCATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5859004‑5859117 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCGGCCTTGGGCGAGCCTTTGATCACTGCGGTTGCCTCGATCTCCCAGCCCAAACCCTCCTTTGGCAGCACGATGTCCAACGGCGCGCCTTGGCGCTTG  >  AE015451/5859158‑5859296                                                                                                                   ||||||||||||||||||||||||||||||||||||||||                                                                                                                                            GCCATTCGGCCAGAATCGGGTCCCGGTTTTTAGACGCCCAGACAAAGTCGTTCTTGATCAACCGCTGCTCATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCGGCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGC‑‑CGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCG                                                                                                  <  2:297988/195‑1                                                           CAACCGCTGCTCATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCGGCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGCCGCGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCGGCCTTGGGCGAGCCTTTGATCACTGCGGTTGCCTCGATCTCCCAGCCCAAACCCTCCTTTGGCAGCACGATGTCCAACGGCGCGCCTTGGTGCTTG  >  1:444554/1‑235                                                                      CATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCGGCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGCCGCGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCGGCCTTGG                                                                                           <  1:655588/135‑1                                                                      CATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCGGCCATGGAACTGTACAAGGAAAACTTCGCCGTACTGGCCGCGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCGGCCTTGG                                                                                           >  2:655585/1‑135                                                                                                                   ||||||||||||||||||||||||||||||||||||||||                                                                                                                                            GCCATTCGGCCAGAATCTGGTCACGGTTTTTAGAAGCCCAGACAAAGTCGTTCTTGATCAACCGCTGCTCATAATCCGCCGGTAGCTCGGTCTGCGGCTTGGCGATGCCCGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5859004‑5859117 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCCGGGCTGGCCGAGAAATCGGCTAGGCGCTTGGCCACATCGGCCTTGGGCGAGCCTTTGATCACTGCGGTTGCCTCGATCTCCCAGCCCAAACCCTCCTTTGGCAGCACGATGTCCAACGGCGCGCCTTGGCGCTTG  >  AE015451/5859158‑5859296

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