New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5397806	55 (0.910)	3 (0.050)	3/246	NT	5.3%	coding (818/1731 nt)	hsdS	type I restriction modification system specificity protein
	?	NC_002947	= 5397858	55 (0.980)					coding (870/1731 nt)	hsdS	type I restriction modification system specificity protein

CAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397650‑5397806 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaactggtCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  <  NC_002947/5397858‑5397717                                                                                                                                                                                                                                                                                                             CAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTC                                                                                                                                                        <  4:267413/149‑1      TTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACT                                                                                                                                                   >  5:150020/1‑149                                                         GAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGC                                                                                                >  4:36465/1‑149                                                                      TCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGG                                                                                   <  3:36465/149‑1                                                                           GCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTTGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCTCCTGCACCAGTTGGACGT                                                                               <  2:12492/148‑1                                                                                                                                                       AACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  <  2:228997/149‑1                                                                                                                                                       AACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  <  2:97270/149‑1                                                                                                                                                       AACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  >  3:315674/1‑149                                                                                                                                                       AACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  >  4:329177/1‑149                                                                                                                                                       AACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  >  7:3640/1‑149                                                                                                                                                                                                                                                                                                             CAACATTAAATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397650‑5397806 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaactggtCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGGACGTGGGCGGACTCGGCGTCGGCCTGCTGGGCTTCCAAGCGGTCGCATAGGGCCATTAGTTCTTCGACTTTGGCAACGATG  <  NC_002947/5397858‑5397717

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