New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3828223 =	7 (0.100)	24 (0.380)	17/250	3.7	75.6%	intergenic (‑204/+155)	PP_3381/PP_3382	transposase/gluconate 2‑dehydrogenase cytochrome c subunit
	?	NC_002947	5224584 =	9 (0.140)					intergenic (‑240/+139)	PP_4603/PP_4604	transposase/DMT superfamily permease
Rejected: Coverage evenness skew score above cutoff.

CCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAA  >  NC_002947/3828211‑3828371             |                                                                                                                                                     cccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGa             <  2:438607/150‑1 (MQ=32)  ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAg            <  1:3765/150‑1 (MQ=32)  ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAg            <  2:452745/150‑1 (MQ=32)    tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGAAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAggg          >  1:936023/1‑150 (MQ=21)       cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTCTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGgcgc       <  1:103666/150‑1 (MQ=37)        aaaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCt      >  1:839071/1‑150 (MQ=38)            cAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGaaa  <  2:1064059/150‑1 (MQ=40)             |                                                                                                                                                     CCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAA  >  NC_002947/3828211‑3828371

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