New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	97 (1.210)	4 (0.050)	4/280	14.0	4.4%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	= 4542986	81 (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:1226986/1‑151  GCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                                                                                              >  2:854819/1‑151   CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                                                                                             >  1:51541/1‑151   CCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                                                                                             >  2:230780/1‑151     ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGG                                                                                                                                                           >  2:475276/1‑151                                        TTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGAC                                                                                                                         <  1:825331/150‑1                                                   ATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAAATTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTC                                                                                                               >  1:1245234/1‑149                                                                                                             TTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCAT                                                   >  2:879433/1‑151                                                                                                              TATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATG                                                  >  2:1042090/1‑151                                                                                                                                                       TTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAG          <  2:144403/150‑1                                                                                                                                                       TTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTAGAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAG          <  1:1154196/150‑1                                                                                                                                                         GGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCC        >  1:871335/1‑150                                                                                                                                                           GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACA     <  1:1153863/151‑1                                                                                                                                                              CATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGA   >  2:852294/1‑150                                                                                                                                                               ATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAA  >  2:311345/1‑150                                                                                                                                                                                                                                                                                                                      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 < 29 ≤ ATCG/ATCG < 34 ≤ 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.