New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	69 (1.050)	4 (0.060)	4/280	11.9	5.8%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	= 4542986	66 (1.070)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

CGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542532‑4542690 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAA  <  NC_000913/4542986‑4542838                                                                                                                                                                                                                                                                                                                      CGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                                                                                               <  2:997961/151‑1   CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                                                                                             <  1:212390/151‑1                                 TGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATC                                                                                                                               >  1:663956/1‑151                                                             AATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTT                                                                                                   >  1:390562/1‑151                                                                           CGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCT                                                                                      >  2:160593/1‑150                                                                                                                                                    CAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCG            <  1:170764/151‑1                                                                                                                                                       TTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAG          <  2:386059/150‑1                                                                                                                                                         GGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCA       <  2:370860/151‑1                                                                                                                                                           GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCAC      >  1:475033/1‑150                                                                                                                                                              CATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAA  <  1:68208/151‑1                                                                                                                                                                                                                                                                                                                      CGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542532‑4542690 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAA  <  NC_000913/4542986‑4542838

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