New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 652785	108 (0.730)	57 (0.460)	39/222	3.4	55.5%	coding (1961/1962 nt)	ctpL	methyl‑accepting chemotaxis protein CtpL
	?	NC_002947	= 4808204	0 (0.000)					intergenic (‑268/‑287)	PP_4234/dsbD‑II	bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Rejected: Coverage evenness skew score above cutoff.

GACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCG  >  NC_002947/4808055‑4808236                                                                                                                                                      |                                 gACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGCGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTgg                                  >  2:3335549/1‑150 (MQ=255)   gaaTGTGGCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTAGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGga                                <  2:2374061/149‑1 (MQ=255)    aaTGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgag                               >  1:2626065/1‑150 (MQ=255)           cgcgGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGGCGGGAGAGTTTCAGCGCCTGGGAGATCGAgc                        <  2:2316511/150‑1 (MQ=255)           cgcgGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAgc                        <  2:2654361/150‑1 (MQ=255)            gcgGTGTGTAGAGGGCAATGGGCGGAATTCGTGCTTTGGCAAGTAGTTGTTTGGGGGGAGGGTTTCAGCGCCTGGGAGGAGTTCGCCGGGGGATTTTCGGCGCCTGGGAGGTGTTCGCCGGGAGGCTTTCAGCGCCTGGGAGATCGAgcg                       >  2:2310197/1‑150 (MQ=37)             cgGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAgcgc                      <  2:1499276/150‑1 (MQ=255)                  tgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGAGGTTTGGGGATCGAGCGccgccc                 >  2:868027/1‑150 (MQ=38)                      taGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgc             <  1:3428420/150‑2 (MQ=255)                        gggCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGATTTCAGCGCCTGGGAGGTGTTCGCCCGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGCGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcgg           >  2:2277731/1‑147 (MQ=255)                        gggCAATGGGCGGAATTCGCGCTTTGGCGGGAAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTACGCCGGGAGAGTTTCCTCGCCTGGGAGATCGCGCGCCGCCcgcgcgg           >  2:1428574/1‑147 (MQ=38)                         ggCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgcgcgggcggc          >  2:896826/1‑144 (MQ=255)                                 gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTGAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg  >  2:1592518/1‑150 (MQ=38)                                 gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgaga                              >  1:2542847/1‑122 (MQ=255)                                 gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgaga                              <  2:2542847/122‑1 (MQ=255)                                                  gCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGga                                <  2:325807/103‑1 (MQ=255)                                                  gCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGga                                >  1:325807/1‑103 (MQ=255)                                                            gtTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATc                            <  2:2298687/97‑1 (MQ=255)                                                            gtTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATc                            >  1:2298687/1‑97 (MQ=255)                                                                                                                                                      |                                 GACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCG  >  NC_002947/4808055‑4808236

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