New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	232 (0.820)	5 (0.020)	5/260	NT	2.1%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	= 4542986	252 (0.950)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Frequency below/above cutoff threshold.

CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542534‑4542690 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGA  <  NC_000913/4542986‑4542839                                                                                                                                                                                                                                                                                                                   CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                                                                                             >  1:3227291/1‑150   ATACTGTGCGTTATACCGCCATTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                                                                                           <  1:1451469/150‑1   ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGTGAAATAAATAATCTAATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                                                                                           <  1:2944699/150‑1   ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                                                                                           >  2:4581289/1‑150   ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                                                                                           <  1:107224/150‑1      CTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                                                                                        <  2:3980898/150‑1       TGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                                                                                       >  1:4078478/1‑150         TGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAT                                                                                                                                                     >  1:3559754/1‑150                                                   GGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGT                                                                                                           <  2:4404928/150‑1                                                              AATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGG                                                                                                <  2:3174289/150‑1                                                                                                           TTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCA                                                   <  1:343632/150‑1                                                                                                                    ATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGT                                          >  2:2459335/1‑150                                                                                                                                   CAAATAAATCCAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAA                           <  2:3559754/150‑1                                                                                                                                                     TTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAG         <  1:1088515/150‑1                                                                                                                                                        GGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCA      <  2:1950746/150‑1                                                                                                                                                        GGGCCATTTGGACTCAGAGAGGATAGCATCGCGGACAATGTTTTTCAGTTTATTGGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCA      <  1:3245703/150‑1                                                                                                                                                         GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCAC     <  2:2362520/150‑1                                                                                                                                                          GCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACA    <  1:2412287/150‑1                                                                                                                                                          GCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACA    >  1:4164696/1‑150                                                                                                                                                          GCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACA    >  1:4214510/1‑150                                                                                                                                                          GCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACA    >  1:1989388/1‑150                                                                                                                                                           CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAG   >  2:1374540/1‑150                                                                                                                                                            CATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGA  <  1:2024450/150‑1                                                                                                                                                                                                                                                                                                                   CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542534‑4542690 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGA  <  NC_000913/4542986‑4542839

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