New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4585312 =	121 (1.080)	4 (0.040)	4/264	NT	3.4%	coding (1450/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
	?	NC_000913	4585357 =	112 (1.050)					coding (1405/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
Rejected: Frequency below/above cutoff threshold.

TTGCCGATAACCCTGGCATGAGCGAGCTGGGCCTGCGTTATTATCAGGAAGATGCCGTCCGCGCGGTTGAAAAGGCAATCGTCAAGGGGCAACAAGAGATCCTGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4585461‑4585312 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGC  >  NC_000913/4585357‑4585505                                                                                                                                                                                                                                                                                                                   TTGCCGATAACCCTGGCATGAGCGAGCTGGGCCTGCGTTATTATCAGGAAGATGCCGTCCGCGCGGTTGAAAAGGCAATCGTCAAGGGGCAACAAGAGATCCTGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCA                                                                                                                                                       <  1:1328399/150‑1   GCCGATAACCCTGGCATGAGCGAGCTGGGCCTGCGTTATTATCAGGAAGATGCCGTCCGCGCGGTTGAAAAGGCAATCGTCAAGGGGCAACAAGAGATCCTGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCAGg                                                                                                                                                     >  1:390014/1‑149                                                                                                      CTGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCG                                                  >  1:1225676/1‑150                                                                                                       TGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGCCGGCATCTTCCTGATAATAACGCAGGCCCAGCTCCCTCATGCCAGGGTTATCGG                                                 >  2:785956/1‑150                                                                                                         CTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCAGCAGGATCTCTTGTTGCCCCTGGACGAGTGCCTTTTCCACCGCGCGGACGGCAGCTTCCTGATAATACCGCCGGACCCGCTCGCTCGTGCCAGGGTTATCGGCA                                               >  1:1832150/1‑150                                                                                                          TGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTAT                                                    >  2:253850/1‑150                                                                                                                                                         GCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCC     >  2:1688760/1‑150                                                                                                                                                         GCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCATCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGAGTCTGTTTTTGCGGGTCGCTGCCGAGCATTTCC     >  2:1616973/1‑150                                                                                                                                                            AGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGC  <  1:182350/150‑1                                                                                                                                                                                                                                                                                                                   TTGCCGATAACCCTGGCATGAGCGAGCTGGGCCTGCGTTATTATCAGGAAGATGCCGTCCGCGCGGTTGAAAAGGCAATCGTCAAGGGGCAACAAGAGATCCTGCTGGCGATGGCGACCGGTACCGGTAAAACCCGTACGGCA‑‑‑‑‑‑ATCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4585461‑4585312 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGC  >  NC_000913/4585357‑4585505

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