New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5397806	52 (1.230)	3 (0.080)	3/240	NT	5.9%	coding (818/1731 nt)	hsdS	type I restriction modification system specificity protein
	?	NC_002947	= 5397858	47 (1.200)					coding (870/1731 nt)	hsdS	type I restriction modification system specificity protein

AATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397658‑5397806 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaactggtCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGG  <  NC_002947/5397858‑5397798                                                                                                                                                                                                                    AATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGT                                                               <  2:56737/149‑1                                                              TCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGG  <  1:241861/149‑1                                                                   GCCAAAGTCGAAGAACTAATGGCCCTACGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCAC         <  1:166076/137‑1                                                                   GCCAAAGTCGAAGAACTAATGGCCCTACGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCAAACTGGTCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCAC         >  2:166076/1‑137                                                                                                                                                                                                                    AATGTTGGCCTGATTCGTAGTGCCGCGGCACCTCTGCCTCCCCTCGCCGAACAACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397658‑5397806 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaactggtCGCGAAGTCGGCGGCGTCGATGGCTTGGGTCAGGCTGTCGAGCATCGCCTGCACCAGTTGG  <  NC_002947/5397858‑5397798

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.