New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	2507756 =	NA (NA)	5 (0.020)	4/272	NT	1.7%	intergenic (+33/‑54)	nupC/insL1	nucleoside (except guanosine) transporter/IS186 transposase
	?	CP009273	= 4533433	290 (0.950)					coding (502/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
Rejected: Frequency below/above cutoff threshold.

GGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATggg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  CP009273/2507862‑2507756 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTAAATGTCGCACCATCCAGCGTCAGCGC  <  CP009273/4533433‑4533328                                                                                                                                                                                                                       GGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAAT                                                                 <  2:227863/150‑1                                      CGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTA                            >  1:3257333/1‑150                                         GAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTAAAT                         <  1:4545195/150‑1                                         GAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTAAAT                         <  1:4545468/150‑1                                                                AATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTAAATGTCGCACCATCCAGCGTCAGCGC  <  1:869558/150‑1                                                                                                                                                                                                                       GGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATggg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  CP009273/2507862‑2507756 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAACTAAATGTCGCACCATCCAGCGTCAGCGC  <  CP009273/4533433‑4533328

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.