New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542682 =	16 (0.280)	5 (0.090)	5/472	3.7	17.7%	intergenic (+49/‑433)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
	?	NC_000913	4542996 =	31 (0.550)					intergenic (+363/‑119)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACA  >  NC_000913/4542496‑4542919                                                                                                                                                                                           |                                                                                                                                                                                                                                               cGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa                                                                                                                                                                                 >  1:66846/1‑250 (MQ=255)      aattaaggattagcgagtacagcGAAATATTCGCCATACTGTGCGTTATAACGCCATTAATGCTGCGCGTTATGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAAGTAAAAAGAGAAGAGGTTTGAGATAACTTATTGATAATAAAGTTAAAAAAACAAAAAAATACAAGAAAATTTGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAGTAAAAAAc                                                                                                                                                                            <  2:359745/228‑1 (MQ=255)                 atatcgggaagcaaaaTATTCGCCAGATTGTGGGTTAATCCGCCAGTAATGCTGCTCGTGTTGACGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAAGAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGTCCAAACTTTCAATATCATAAATAAGTTACGTATTTTTTCGCAGGCAGAAAAATATTAATAAAC‑GACAAAAAGCa                                                                                                                                                                <  2:367539/237‑1 (MQ=255)                                              gtgtgTTATACCGGCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAATAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAattattgct                                                                                                                                   <  1:204280/246‑4 (MQ=255)                                                    tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTa                                                                                                                             >  1:214871/1‑250 (MQ=255)                                                     ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaa                                                                                                                            >  2:534659/1‑250 (MQ=255)                                                                      gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                          <  1:292503/251‑1 (MQ=255)                                                                                   gATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTa                                                                                              >  1:147571/1‑250 (MQ=255)                                                                                            aaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGt                                                                                    >  1:516270/1‑251 (MQ=255)                                                                                               gaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTc                                                                                 >  1:4864/1‑251 (MQ=255)                                                                                                                          aaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCAttct                                                      >  1:373314/1‑251 (MQ=255)                                                                                                                                 agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTTTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTcc                                               <  1:351870/251‑1 (MQ=255)                                                                                                                                    agtttgaatttaattaattaataaaaaGTGAAAAAAACAAAAAAATATAAAACAATTGGGGTCAAACTTGCCATAACAAAAATAAGTTACTGATTTTTTATCAAGCATAAAAATATTAAAAAACTGCCAAAAAAC‑TCTTACCGTTTGATATGGAAATTATTTCTATTTTAAATTAACTTCCCATCTCCCTCGCTATATTTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                            <  2:423452/227‑1 (MQ=255)                                                                                                                                                     tGATAATAAAGTTAAAAAAACAAAAAAATACAAGACAATTGGTGCCAAACTGTCCATATAATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑ATCAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTTTAAATTAATGTCACATCCCCTCCGCTATATGTAAATCTAACGTTTCTGCGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCaag                            <  2:357002/250‑2 (MQ=255)                                                                                                                                                                            ataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCCTGCAATATAAACATCTATAAATAAAGATa     >  1:228443/1‑250 (MQ=255)                                                                                                                                                                               aataccagaacaTTGGGGCCACACTGTCCATCTCATAACTAAGTTACGTCTTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACACAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa  >  2:249996/12‑250 (MQ=255)                                                                                                                                                                                           |                                                                                                                                                                                                                                               CGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAC‑GACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACA  >  NC_000913/4542496‑4542919

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