New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 652782	71 (1.010)	14 (0.230)	10/240	5.6	30.5%	coding (1961/1962 nt)	ctpL	methyl‑accepting chemotaxis protein CtpL
	?	NC_002947	= 4808159	3 (0.050)					intergenic (‑268/‑287)	PP_4234/dsbD‑II	bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Rejected: Coverage evenness skew score above cutoff.

CAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGG  >  NC_002947/4808012‑4808193                                                                                                                                                       |                                   cattgttTCGCGGTGTGTATAGGGCATTGGGCGGAATTCGCGCGTTGGCGGGTAGGTCTTTGGCGGGAGGTTTTCAGC‑GGCTGGG‑AGG‑TGTTCGCCGGGGGAGTTTCAGCGCCTGGGAGGGGTTCGCCGGGAGAGTTTCAGCGCCTggg                                   <  1:327552/146‑1 (MQ=255)          gcgGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGGCTGGGAGATCGAgc                          >  1:260623/1‑149 (MQ=255)                     aGGTCAATGGGCTGAATTGGCGCTTTGCGGGGTTTGTGTTTGGCGGGATTTTTTCTCCATCCGGGGTAGGCTCATCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccc                   <  1:329268/148‑1 (MQ=16)                      tggCAATGGTCGGAATTAGCTCTTTGGCGGTTAGGTTTTAGTCGGGAGGTTTTCAGC‑GCCTTGT‑TGG‑TGTTCGGCGGGTTATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgc               <  5:295088/147‑2 (MQ=25)                        tctttGGGTGGATTTGGCGCTTTGGGGGTTAGGTTTTGGGAGTGAGTTTTTCAGT‑GCCTGGT‑AGG‑TGTTGGGCGGTTGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcgg             <  8:325562/144‑4 (MQ=18)                        gCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGTTCGCCGGGGGAGTTTCAGCGACTGGGAGGTGTTCGCCGGGAGAGTTGGATCGCCTGGGAGATCGAGCGCCGCCCGCGAGgc            >  3:52392/1‑149 (MQ=38)                          aaTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTGGCCGGGAGAGTGTCAGCGCCTCGGAGATCGAGCGCCGCCCGCGCGgcgc          >  7:30274/1‑149 (MQ=38)                              ggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGGTCGCCGGGGGATTTTCAGCGCCTGGGAGCTGTTCGCCCGGAGAGTTTCATCGCCAGGGATATAGAGGGACGCACGCGAGgcgcttcc      >  7:156950/1‑145 (MQ=255)                              ggCGCAATTCGCGCTGTCGCGGGTAGGCTGTTGGCGGCAGGTGTTCTGG‑GTCTGGG‑TGG‑TGTTCGGCGGGTGTTTTGCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcg      <  8:30274/149‑5 (MQ=16)                              gcccgAAGTCTCGCTTTGGCGGGTCGGTCTTGGGCGGCAGGCTCGCCGC‑GCCTGGG‑AGG‑TGTTCGCCGGCTGTTGTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttc       <  2:260623/144‑4 (MQ=21)                                  gAATTGGGGCGTTGGGGGTTAGGTTTTTGGCGGGAGGTTTGCTGG‑TCCGGGT‑AGG‑TGTTGGCCGGTTGTTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg  <  2:92761/149‑1 (MQ=16)                                                                                                                                                       |                                   CAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGC‑GCCTGGG‑AGG‑TGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGG  >  NC_002947/4808012‑4808193

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