New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1043472	66 (0.920)	4 (0.060)	4/256	NT	6.0%	coding (1282/1311 nt)	colS	integral membrane sensor signal transduction histidine kinase
	?	NC_002947	= 1043490	63 (0.920)					coding (1300/1311 nt)	colS	integral membrane sensor signal transduction histidine kinase

CGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1043321‑1043472 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgaccaGGGCGTCCGGGTCGCCTTTGGTCGCGG  <  NC_002947/1043490‑1043464                                                                                                                                                                                     CGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCG                                 <  8:410293/148‑1   TGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACC                              <  2:32797/149‑1                                                         AGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACCAGGGCGTCCGGGTCGCCTTTGGTCGC    <  3:36513/121‑1                                                         AGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACCAGGGCGTCCGGGTCGCCTTTGGTCGC    >  4:36513/1‑121                                                                    GACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACCAGGGCGTCCGGGTCGCCTTTGGTCGCGG  <  2:109365/112‑1                                                                    GACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGACGCTTCCAGGTGGACCTGAGCAACACCGCGACCAGGGCGTCCGGGTCGCCTTTGGTCGCGG  >  1:109365/1‑112                                                                                                                                                                                     CGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACATCCACCTTGCCGCATGGCTGCCGCTTCCAGGTGGACCTGAGCAACACCGCGACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1043321‑1043472 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgaccaGGGCGTCCGGGTCGCCTTTGGTCGCGG  <  NC_002947/1043490‑1043464

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.