Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	1,293,006	IS1 (–) Δ11 bp	16.1%	intergenic (‑84/‑511)	hns ← / → tdk	global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	257908 =	NA (NA)	8 (0.030)	8/278	NT	4.9%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 1293005	167 (0.540)					intergenic (‑83/‑522)	hns/tdk	global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
Rejected: Frequency below/above cutoff threshold.
*	?	NC_000913	1293017 =	167 (0.540)	56 (0.200)	30/278	NT	26.6%	intergenic (‑95/‑510)	hns/tdk	global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
	?	NC_000913	= 1979270	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

GAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/258056‑257908 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGAGCGAAGCACTTAAAATTCTGAACAACATCCGTACTCTTCGTGCGCAGGCAA  <  NC_000913/1293005‑1292868                                                                                                                                                                                                                                                                                                 GAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                                                           >  1:145878/1‑150       AGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCT                                                                                                                                     >  1:4383255/1‑150               AACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCTGGCTATTG                                                                                                                             <  2:2820433/150‑1                                       GCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATT                                                                                                     <  1:3053375/150‑1                                               GGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGATCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTC                                                                                             <  1:2998377/150‑1                                                                                       CATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGA                                                     >  1:1934291/1‑150                                                                                                       GGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGAGCGAAGCACTTAAAAT                                     >  1:1478207/1‑150                                                                                                                                       GTTGGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGAGCGAAGCACTTAAAATTCTGAACAACATCCGTACTCTTCGTGCGCAGG     >  2:4635240/1‑150                                                                                                                                          GGAGTCATTACCCATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGAGCGAAGCACTTAAAATTCTGAACAACATCCGTACTCTTCGTGCGCAGGCAA  <  2:4383255/150‑1                                                                                                                                                                                                                                                                                                 GAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/258056‑257908 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CATTCCTGGCTATTGCACAACTGAATTTAAGGCTCTATTATTACCTCAACAAACCACCCCAATATAAGTTTGAGATTACTACAATGAGCGAAGCACTTAAAATTCTGAACAACATCCGTACTCTTCGTGCGCAGGCAA  <  NC_000913/1293005‑1292868

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