New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542682 =	131 (1.340)	7 (0.080)	7/280	15.5	5.6%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	4542996 =	112 (1.220)					intergenic (+363/‑119)	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.

GAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4542840‑4542682 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTC  >  NC_000913/4542996‑4543145                                                                                                                                                                                                                                                                                                                       GAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTT                                                                                                                                                                >  2:613173/1‑151    ACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGC                                                                                                                                                             >  2:89086/1‑151         AGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCA                                                                                                                                                         <  1:1372837/150‑1                                    AAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGT                                                                                                                             <  1:516921/151‑1                                                AATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCA                                                                                                                 <  1:69092/151‑1                                                      AATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAA                                                                                                           <  1:661583/151‑1                                                                               ATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAC                                                                                    >  1:651667/1‑149                                                                                                     ATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTT                                                             <  1:1227949/150‑1                                                                                                                TGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTT                                                  <  2:651640/150‑1                                                                                                                           ACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAA                                       <  2:549982/150‑1                                                                                                                                                       TGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT          >  1:1010794/1‑151                                                                                                                                                            cccAATGATTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAATATTAAAACTCTGGCAATCGTTG     >  2:416427/4‑151                                                                                                                                                            CCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTT      >  1:163975/1‑150                                                                                                                                                              AAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT   <  1:311922/151‑1                                                                                                                                                              AAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT   >  2:1501278/1‑151                                                                                                                                                               AATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTC  <  1:426608/151‑1                                                                                                                                                                                                                                                                                                                       GAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4542840‑4542682 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTC  >  NC_000913/4542996‑4543145

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