New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 652785	8 (0.610)	4 (0.360)	4/234	1.0	53.9%	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

GGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCC  >  NC_002947/652639‑652902                                                                                                                                                   |                                                                                                                      acaaCCGCATTAAAAACAAACGGCGCATGGCGGACGCGCGAGGGGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTCGGAGCCGGTTTACCCGGAAAGCGCCGCGCGGGCGGCGCTCGATCTCAGa                                                                                                                    <  1:212229/148‑1 (MQ=25)              aaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACACTGAACGCtgg                                                                                                       >  1:84118/1‑148 (MQ=255)                                               aaCGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcag                                                                      >  1:170639/1‑150 (MQ=255)                                                                                         cGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatccc                                                                     <  1:270725‑M1/109‑51 (MQ=255)                                                                                         cGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatccc                                                                     >  2:270725‑M1/1‑59 (MQ=255)                                                                                          gAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCATGTTTACGCCCCTCGCCAATGCGCTGCAgg                           >  1:151756/1‑150 (MQ=255)                                                                                                 acGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctaccc                    <  1:91966‑M1/150‑100 (MQ=255)                                                                                                           gggTTTTCCCGCGAAGCGCCGCGCGGGCGGTGCTCGATCTCccaggcgctgaaactcgcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgc          <  1:43897‑M1/150‑110 (MQ=255)                                                                                                              tttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAGCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCACTGCGCGGCAGGTTGTCTTCACCCAGTTCGTg       >  1:285493/1‑150 (MQ=255)                                                                                                                   ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGc                              >  1:155136/1‑122 (MQ=255)                                                                                                                   ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGc                              <  2:155136/122‑1 (MQ=255)                                                                                                                   ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACACTCCCGTATATTCCCATTGCATTCAGCTCACCCTCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATcc  >  1:85225/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                      GGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCC  >  NC_002947/652639‑652902

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