Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	867,950	45 bp→TTCGGTA	9.9%	intergenic (+48/+85)	PP_0750 → / ← mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 867949	31 (1.200)	3 (0.120) +TTCGGTA	3/304	NT	9.9%	intergenic (+47/+130)	PP_0750/mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase
	?	AE015451	867995 =	26 (1.010)					intergenic (+93/+84)	PP_0750/mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase

GATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/867878‑867949 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAGATTTCCTTGAG  >  AE015451/867995‑868210                                                                         |||||||                                                                                                                                                                                                                         GATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGCCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTCGGTAATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCC                                                                                                                                                                                    <  1:469116/117‑1 GATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGCCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTCGGTAATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCC                                                                                                                                                                                    >  2:469112/1‑117                                                            CCGCTCCCACAGGTTCGGTAATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAGATTTCCTTGAG  >  1:452481/1‑236                                                                         |||||||                                                                                                                                                                                                                         GATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/867878‑867949 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAGATTTCCTTGAG  >  AE015451/867995‑868210

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