New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4908303	7 (0.200)	7 (0.210)	5/284	3.8	50.5%	intergenic (‑279/‑3)	PP_4317/PP_4318	hypothetical protein/transposase
	?	NC_002947	= 4908336	NA (NA)					coding (31/1302 nt)	PP_4318	transposase
Rejected: Coverage evenness skew score above cutoff.

TTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTT  >  NC_002947/4908157‑4908322                                                                                                                                                   |                    tttttGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGCCCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCg                     >  1:198554/1‑147 (MQ=17)        cgcgTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGcc             <  6:180135/148‑1 (MQ=32)        cgcgTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGc              <  8:37950/147‑1 (MQ=32)            tttGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAAc         >  4:87781/1‑148 (MQ=32)             ttGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACt        >  3:172970/1‑148 (MQ=32)               gACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACtt       <  5:179638/147‑1 (MQ=32)                 cAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGcc    <  7:91954/148‑1 (MQ=32)                   aaaaGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCtt  <  4:172929/148‑1 (MQ=32)                                                                                                                                                   |                    TTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTT  >  NC_002947/4908157‑4908322

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