New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3265438	132 (0.800)	6 (0.060) +GCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAA	5/174	NT	5.5%	coding (1074/1635 nt)	PP_2861	methyl‑accepting chemotaxis transducer
	?	NC_002947	= 4457258	188 (1.140)					coding (1693/2178 nt)	PP_3950	methyl‑accepting chemotaxis transducer

ATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3265391‑3265438 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAGACCAACCTGCTGGCGCTCAATGCCGCCATCGAGGCGGCGCGTGCCGGGGAGCATGGCCGT  <  NC_002947/4457258‑4457196                                                 |||||||||||||||||||||||||||||||||||||||||||||||                                                                ATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCGAGGCGGCACGGGCTGGCGAGC          >  1:5075155/1‑150         GCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCGAGGCGGCACGGGCTGGCGAGCAAGGCCGg  >  2:612371/1‑149         GCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCGAGGCGGCACGGGC                  <  1:2512988/134‑1         GCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCGAGGCGGCACGGGC                  >  2:2512988/1‑134                 CGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCG                               >  1:755311/1‑113                 CGGCCAGGTCGCGGATTACCTGTACGGCCACAGCACGCGCATCGGCGGCGTGCTGGACGTGATTCGTTCCATTGCCGAACAGACCAACCTGCTGGCGCTGAATGCCGCCATCG                               <  2:755311/113‑1                                                 |||||||||||||||||||||||||||||||||||||||||||||||                                                                ATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3265391‑3265438 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAGACCAACCTGCTGGCGCTCAATGCCGCCATCGAGGCGGCGCGTGCCGGGGAGCATGGCCGT  <  NC_002947/4457258‑4457196

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