New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4962552 =	31 (1.150)	6 (0.250)	5/264	2.4	16.9%	intergenic (‑41/+122)	fliH/fliG	flagellar assembly protein FliH/flagellar motor switch protein
	?	NC_002947	4962642 =	31 (1.280)					intergenic (‑131/+32)	fliH/fliG	flagellar assembly protein FliH/flagellar motor switch protein
Rejected: Frequency below/above cutoff threshold.

GCTCGGCGGCAAGGGCGCCGAGGAAATGATCTAAGGTATTCGCGATCCCATAGCCGACGCCGCCTTTTGTAGGAGCGGGTTTACCCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4962707‑4962552 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aacccgctcctacaaAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCG  >  NC_002947/4962642‑4962789                                                                                                                                                                                                                                                                                                                  GCTCGGCGGCAAGGGCGCCGAGTAAATGATCTAAGGTATTCGCGATCCCATAGCCGACGCCGCCTTTTGTAGGAGCGGGTTTACCCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTC                                                                                                                                                            <  1:245912/150‑1   TCGGCGGCAAGGGCGCCGAGGAAATGATCTAAGGTATTCGCGATCCCATAGCCGACGCCGCCTTTTGTAGGAGCGGGTTTACCCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCT                                                                                                                                                          >  2:377511/1‑150                                                                                     CCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAAAAGGCGGCGTCGGCTCTGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCG                                                                        >  2:539733/1‑150                                                                                     CCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCG                                                                        >  2:185359/1‑150                                                                                         GAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCC                                                                         >  1:219839/1‑145                                                                                         GAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCC                                                                         <  2:219839/145‑1                                                                                                        GGCAATCGAGGTTGTCGCTACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGC                                                           <  1:516964/144‑1                                                                                                        GGCAATCGAGGTTGTCGCTACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGC                                                           >  2:516964/1‑144                                                                                                                                                   GCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGAC          >  1:425752/1‑150                                                                                                                                                           AAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCG  >  1:70185/1‑150                                                                                                                                                                                                                                                                                                                  GCTCGGCGGCAAGGGCGCCGAGGAAATGATCTAAGGTATTCGCGATCCCATAGCCGACGCCGCCTTTTGTAGGAGCGGGTTTACCCGCGAAGGCGTCGGTCCTGGCAATCGAGGTTGTCGCTGCCGACGCATTCGCGGGTGAACCCGCTCCTACAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4962707‑4962552 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aacccgctcctacaaAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCG  >  NC_002947/4962642‑4962789

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.