New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3864155	46 (0.870)	4 (0.080)	4/242	NT	7.6%	coding (3078/3282 nt)	PP_3413	two‑component system sensor histidine kinase/response regulator
	?	NC_002947	= 3864199	53 (1.050)					coding (3122/3282 nt)	PP_3413	two‑component system sensor histidine kinase/response regulator

AGTGGGTTCGACCTGGCCGAACTGCGTCACCTGACCCAGGACGACGAACAGTTGACGCGGCATTTGCTGGAGCAATTGTCGATAAGCGTCGGTGAAGATCTGGCCGCCCTGCGTGCCTTGTCAACAGACGCCCCTGATGAAACCGTGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3864006‑3864155 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgcgcACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGGCAAGGGCGCGCA  <  NC_002947/3864199‑3864151                                                                                                                                                                                                         AGTGGGTTCGACCTGGCCGAACTGCGTCACCTGACCCAGGACGACGAACAGTTGACGCGGCATTTGCTGGAGCAATTGTCGATAAGCGTCGGTGAAGATCTGGCCGCCCTGCGTGCCTTGTCAACAGACGCCCCTGATGAAACCGTGCG                                                    >  4:231569/1‑149                                        GACGACGAACAGTTGACGCGGCATTTGCTGGAGCAATTGTCGATAAGCGTCGGTGAAGATCTGGCCGCCCTGCGTGCCTTGTCAACAGACGCCCCTGATGAAACCGTGCGCACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGG              >  2:147403/1‑148                                                   GTTGACGCGGCATTTGCTGGAGCAATTGTCGATAAGCGTCGGTGAAGATCTGGCCGCCCTGCGTGCCTTGTCAACAGACGCCCCTGATGAAACCGTGCGCACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGGCAAGGGCGCGCA  <  1:147403/149‑1                                                                                                                     CTTGTCAACAGACGCCCCTGATGAAACCGTGCGCACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGGCAAGGGCGC     <  1:45248/80‑1                                                                                                                     CTTGTCAACAGACGCCCCTGATGAAACCGTGCGCACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGGCAAGGGCGC     >  2:45248/1‑80                                                                                                                                                                                                         AGTGGGTTCGACCTGGCCGAACTGCGTCACCTGACCCAGGACGACGAACAGTTGACGCGGCATTTGCTGGAGCAATTGTCGATAAGCGTCGGTGAAGATCTGGCCGCCCTGCGTGCCTTGTCAACAGACGCCCCTGATGAAACCGTGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3864006‑3864155 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgcgcACCCTGATCATTTTGGCGCCACCCTTGATGCGGTGGGCAAGGGCGCGCA  <  NC_002947/3864199‑3864151

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