New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	4542683 =	38 (0.670)	9 (0.170)	7/260	6.1	18.6%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_bme_pgi	4542997 =	43 (0.810)					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/above cutoff threshold.

CGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/4542837‑4542683 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCG  >  NC_000913_3_bme_pgi/4542997‑4543140                                                                                                                                                                                                                                                                                                             CGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTATTATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGG                                                                                                                                                        <  2:997867/149‑1  gttcgaTTTACATATAGCGGATGTAATTTGAACTTAAGTTATGATAGAAATAATTTACATATCAAATAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGC                                                                                                                                                       <  2:485475/143‑1  GTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAGATACGTAACTTATTTATGATATGGACAGTTTGac                                                                                                                                                       >  1:833765/1‑147       CTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGCAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA                                                                                                                                                  >  1:985787/1‑149       CTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA                                                                                                                                                  >  1:596031/1‑149                 aaCGGAGTTGATGTTAAATTAATTTACAATTGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATC                                                                                                                                        <  2:290864/147‑1                                            AATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGG                                                                                                             <  2:897392/149‑1                                                                   ACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTG                                                                                                                                   <  1:876864/104‑1                                                                   ACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTG                                                                                                                                   >  2:876864/1‑104                                                                                                           TTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGT                                              >  1:181733/1‑149                                                                                                               GAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTT                                          <  1:305078/149‑1                                                                                                               GAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTT                                          <  1:305042/149‑1                                                                                                               GAAAAAATAAGTTACGTATTTATGATATGGACATATTGGTCTCAAATGTTTAATATTTTGGGGGAAAACTGTGCAGTGTTGGTAGTCAAACTCGTTTATAAAATAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTT                                          <  1:305077/149‑1                                                                                                                    aagACGTAACTTATTTAAGATATGGACAGTTTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAA                                     <  1:346304/146‑1                                                                                                                                                     GCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT    <  2:770842/149‑1                                                                                                                                                      CCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATC   <  2:780308/149‑1                                                                                                                                                       CCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCG  >  1:43949/1‑149                                                                                                                                                                                                                                                                                                             CGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATATTTTTATGCTTGAGAAAAAATACGTAACTTATTTATGATATGGACAGTTTGGCCCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/4542837‑4542683 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCG  >  NC_000913_3_bme_pgi/4542997‑4543140

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