New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3544651	31 (0.710)	5 (0.120)	5/280	2.8	14.5%	coding (194/771 nt)	bioH	pimeloyl‑ACP methyl ester carboxylesterase
	?	NC_000913	= 3544692	30 (0.730)					coding (153/771 nt)	bioH	pimeloyl‑ACP methyl ester carboxylesterase
Rejected: Frequency below cutoff threshold.

CTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAGCGCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3544497‑3544651 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcttcggGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTTAACCCATCCCGAGCGTGTTCAGG  <  NC_000913/3544692‑3544544                                                                                                                                                                                                                                                                                                                  CTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAGCGCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGG                                                                                                                                                             >  1:620347/1‑149                                                              GGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCCCCCCGACCCCAGCCTAACCAAATGGCTTTATCAAGGGCCCGTTGCAGCACGGCTTCCGGCGGAGCCGGGGGTTTTGTGCGCTGTCCCTTGCTGATATGGCCGAAGCCGGGCTGCA                                                                                              >  2:533/1‑151                                                                                 GGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGC                                                                            >  2:334628/1‑150                                                                                                       ACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGG                                                     >  2:488275/1‑151                                                                                                                         GGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGC                                    >  2:148541/1‑150                                                                                                                                     TGCCTGTTGCAGCACGGCTTCGGGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTTAA                       >  2:310424/1‑151                                                                                                                                                           GGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTTAACCCATCCCGAGCagggtaagg  <  2:539762/150‑10                                                                                                                                                                                                                                                                                                                  CTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAGCGCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3544497‑3544651 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcttcggGCGTAGCCGGGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTTAACCCATCCCGAGCGTGTTCAGG  <  NC_000913/3544692‑3544544

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