New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 652785	37 (1.150)	7 (0.260)	5/224	1.8	30.8%	coding (1961/1962 nt)	ctpL	methyl‑accepting chemotaxis protein CtpL
	?	NC_002947	= 4810126	0 (0.000)					intergenic (‑268/‑287)	PP_4234/dsbD‑II	bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2

GTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTG  >  NC_002947/4809991‑4810184                                                                                                                                        |                                                           gtgtgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCCGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTgg                                                            >  2:701766/1‑136 (MQ=255)                 gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcc                                            <  1:295993‑M2/136‑17 (MQ=255)                 gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcc                                            <  1:330378‑M2/136‑17 (MQ=255)                 gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcc                                            <  1:800619‑M2/136‑17 (MQ=255)                 gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcc                                            <  1:531142‑M2/136‑17 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccg                                          <  1:791627‑M2/136‑19 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccg                                          <  1:800596‑M2/136‑19 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccg                                          <  1:75127‑M2/136‑19 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccc                                          <  1:115968‑M2/136‑19 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGATCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccc                                          <  1:439125‑M2/136‑19 (MQ=255)                   gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGCGAATTTCAGCGCCTGGGAGATGTTCGCCCGGAGCCATTCAGCTACTGGgagatcgaccgccgcacg                                          <  2:91642‑M2/136‑19 (MQ=255)                                                          tttttcAGCGCCTGGCTGGTGTGCGCCGGGGGATTTTCAGCGCCTGGTAGGTGTTCGCCGGGAGAGTTTCCGCGCCTGGgagatcgaccgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgttt   >  2:75127‑M2/2‑79 (MQ=255)                                                          gTTTTCTGCGCGTGGGTGGTGTTCGCCGGGGGATTTTCTGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgttt   >  1:393598‑M2/1‑79 (MQ=255)                                                           ttttCTGCGGGTGGTTGGTGTGCGCCGGGTGATTTTCAGCGCCTGGTAGGTGTTCGCCGGTTGTTTTTCAGCGCCTGGgatatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgtttc  >  1:91642‑M2/1‑78 (MQ=255)                                                                               ttCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgttt   >  1:52665‑M2/1‑58 (MQ=255)                                                                               ttCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGgagatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgttt   <  2:52665‑M2/115‑58 (MQ=255)                                                                                                                                    ccTGGgagatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgt     >  1:557759‑M2/1‑5 (MQ=255)                                                                                                                                    ccTGGgagatcgagcgccgcccgcgcggcgcttcgcgggtaaacccgctcctacgtttgt     <  2:557759‑M2/60‑56 (MQ=255)                                                                                                                                        |                                                           GTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTG  >  NC_002947/4809991‑4810184

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