New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1190141	43 (0.820)	4 (0.080)	4/236	NT	8.2%	coding (199/234 nt)	tatA‑I	twin arginine protein export complex subunit
	?	NC_002947	= 1190161	48 (0.950)					coding (179/234 nt)	tatA‑I	twin arginine protein export complex subunit

CGCGCCTGGCCCAGGCCGCGGCCGAGGGTACGGGCTGCCACCGGCAGGCGCTCGGGGCCCAGCACCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1189988‑1190141 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaacAGTGCCCGCTGAACGGCCAAGTTGCCCAGCA  <  NC_002947/1190161‑1190131                                                                                                                                                                                           CGCGTCTGGCCCAGGCCGCGGCCGAGGGTACGGGCTGCCACCGGCAGGCGCTCGGGGCCCAGCACCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGG                                     <  2:1324693/150‑1                                                                ACCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGC                                    >  1:306735/1‑88                                                                ACCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGC                                    <  2:306735/88‑1                                                                 CCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGGTTGCTGGGCAACAGTGCCCGCTGAACGGCCAAGTTGCCCAGCA  >  1:816729/1‑121                                                                 CCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGGTTGCTGGGCAACAGTGCCCGCTGAACGGCCAAGTTGCCCAGCA  <  2:816729/121‑1                                                                                                                TCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGCAACAGTGCCCGCTGAACGGCCAA             >  1:79818/1‑63                                                                                                                TCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGCAACAGTGCCCGCTGAACGGCCAA             <  2:79818/63‑1                                                                                                                                                                                           CGCGCCTGGCCCAGGCCGCGGCCGAGGGTACGGGCTGCCACCGGCAGGCGCTCGGGGCCCAGCACCAACAGCGCGACGACGCCGATCAGCAGCAACTCGCTGAAGCCTACCTCGAACATCAGGCCTGGCGGTCCGCTTGAGATTGCTGGGCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1189988‑1190141 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaacAGTGCCCGCTGAACGGCCAAGTTGCCCAGCA  <  NC_002947/1190161‑1190131

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