New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3388268	32 (0.650)	4 (0.090)	4/238	NT	11.5%	coding (74/480 nt)	PP_2990	MerR family transcriptional regulator
	?	NC_002947	= 3388288	31 (0.660)					coding (94/480 nt)	PP_2990	MerR family transcriptional regulator

ACACGTTGTGCTTGGCCGCCCGATGACTTCAGGTTTAAAGTCAACTTTAAGGTCAATAGGTTTTTCAGGGTGAGGTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3388119‑3388268 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcctTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAATGC  <  NC_002947/3388288‑3388241                                                                                                                                                                                                        ACACGTTGTGCTTGGCCGCCCGATGACTTCAGGTTTAAAGTCAACTTTAAGGTCAATAGGTTTTTCAGGGTGAGGTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCACTGGCCT                                                  >  2:935016/1‑150      TTGTGCTTGGCCGCCCGATGACTTCAGGTTTAAAGTCAACTTTAAGGTCAATAGGTTTTTCAGGGTGAGGTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCACTGGCCTcatag                                             <  1:935016/150‑6                                                                           GTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCTTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAATGC  >  1:185986/1‑124                                                                           GTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCTTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAATGC  <  2:185986/124‑1                                                                                                 AAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCTTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAAT    <  1:1162527/100‑1                                                                                                 AAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCTTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAAT    >  2:1162527/1‑100                                                                                                                                                                                                        ACACGTTGTGCTTGGCCGCCCGATGACTTCAGGTTTAAAGTCAACTTTAAGGTCAATAGGTTTTTCAGGGTGAGGTATGAGAATAGGTGAACTGGCAAAAATCAGTGGACTGGCACCGTCACGCATTCGATTCTATGAGGCCAGTGGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3388119‑3388268 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcctTGCGCTCGACGGAGCTGATCAGGCCACTGGCCTCATAGAATCGAATGC  <  NC_002947/3388288‑3388241

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