New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3653230	NA (NA)	5 (0.050)	5/494	7.0	17.2%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	4580649 =	24 (0.250)					coding (814/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3652992‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCG  >  NC_000913/4580649‑4580810                                                                                                                                                                                                                                                                                                                                                                                                                   CAGCGGAGAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACATCCACCATGTTTTGCCATGGCAGAATCTGCTCCAGTCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACCCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAACCGAA                                                                                                                                                        <  2:351274/251‑1                                       AAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGCTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAACCGAATATCGTTTTGATCTACATGGGCAGCACGTACAGAAat                                                                                                                   >  1:993955/1‑248                                                                                                                                                       GAAGGTAAGATGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAACCGAATATCGTTTTGATCTACATGGGCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAATACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACGTCTTAAATACAGAATGTTGCG  <  2:993955/251‑1                                                                                                                                                                                              GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAAT                                                   >  1:291243/1‑163                                                                                                                                                                                              GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAAT                                                   <  2:291243/163‑1                                                                                                                                                                                                                                                                                                                                                                                                                   CAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3652992‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCG  >  NC_000913/4580649‑4580810

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