Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	1,293,015	IS1 (–) +8 bp	8.6%	intergenic (‑93/‑505)	hns ← / → tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 291401	NA (NA)	4 (0.040)	4/176	NT	5.0%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	1293015 =	69 (0.810)					intergenic (‑93/‑512)	hns/tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase
Rejected: Frequency below/above cutoff threshold.
*	?	NC_000913	= 1050545	NA (NA)	9 (0.100)	8/176	NT	10.6%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 1293022	69 (0.810)					intergenic (‑100/‑505)	hns/tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase

CACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/291363‑291401 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGATGATTTCG  >  NC_000913/1293015‑1293111                                                                                                                                          CACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAA                                       >  2:234543/1‑99         TACACTAAATCAGTAAGTTGGCAGCATCACCCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGC                               >  2:1908935/1‑99                           TGGCAGCATCACCCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACAC             <  1:1362323/99‑1                               AGCATCACCCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGAT         >  1:674260/1‑99                                      CCCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGATGATTTCG  >  2:2315652/1‑99                                                                                                                                          CACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/291363‑291401 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGATGATTTCG  >  NC_000913/1293015‑1293111

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