New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4586194 =	70 (0.780)	4 (0.050)	4/266	NT	5.7%	coding (568/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
	?	NC_000913	4586241 =	64 (0.730)					coding (521/3513 nt)	hsdR	endonuclease R Type I restriction enzyme

AAAACCTCTATCACCAGGAAGTGCTGACGCTAAAACAACAGCTTGAACAGCAGGTGCGAGAAAAAGCGCAGACTCAGGCAGAAGTCGAAGCGCAACAGCAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4586343‑4586194 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAAT  >  NC_000913/4586241‑4586388                                                                                                                                                                                                                                                                                                            AAAACCTCTATCACCAGGAAGTGCTGACGCTAAAACAACAGCTTGAACAGCAGGTGCGAGAAAAAGCGCAGACTCAGGCAGAAGTCGAAGCGCAACAGCAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAG                                                                                                                                                       >  2:59307/1‑149                                                                                        AAGCGCAACAGCAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTC                                                                >  1:32512/1‑149                                                                                             CAACAGCAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGT                                                           <  2:32512/149‑1                                                                                                   CAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGA                                                     >  2:845266/1‑149                                                                                                                         GCTATATCGCCATTCTGGAAGGCAAACAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTTCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAA                               <  1:845266/149‑1                                                                                                                                                   CAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCAT     >  1:1031051/1‑149                                                                                                                                                    AGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATA    >  1:492977/1‑149                                                                                                                                                     GCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAA   <  1:1202895/149‑1                                                                                                                                                      CTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTATCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAAT  <  1:465053/149‑1                                                                                                                                                                                                                                                                                                            AAAACCTCTATCACCAGGAAGTGCTGACGCTAAAACAACAGCTTGAACAGCAGGTGCGAGAAAAAGCGCAGACTCAGGCAGAAGTCGAAGCGCAACAGCAGAAGCTGGTTGCCCTGAACGGCTATATCGCCATTCTGGAAGGCAAACAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4586343‑4586194 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAAT  >  NC_000913/4586241‑4586388

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.