New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3828180 =	5 (0.070)	10 (0.220) +CTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCA	9/178	4.9	34.3%	intergenic (‑204/+155)	PP_3381/PP_3382	transposase/gluconate 2‑dehydrogenase cytochrome c subunit
	?	NC_002947	5545626 =	54 (0.770)					intergenic (‑766/+111)	rpsF/rlmB	30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Coverage evenness skew score above cutoff.

CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAA  >  NC_002947/3828169‑3828328            |                                                                                                                                                     ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGc               >  3:308131/1‑147 (MQ=11)   tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAgg           >  6:3095/1‑149 (MQ=33)   tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAgg           >  6:9202/1‑149 (MQ=33)     aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGgc         <  3:54963/149‑1 (MQ=35)        aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCt      >  6:408796/1‑149 (MQ=38)          agaaaaGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTCCACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTg     >  5:301117/3‑148 (MQ=11)            aaaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGaaa  <  3:160916/149‑1 (MQ=40)            |                                                                                                                                                     CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAA  >  NC_002947/3828169‑3828328

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.