New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4541255	118 (0.920)	11 (0.090)	8/168	NT	8.8%	coding (299/603 nt)	fimB	tyrosine recombinase/inversion of on/off regulator of fimA
	?	NC_000913	= 4541272	117 (0.970)					coding (316/603 nt)	fimB	tyrosine recombinase/inversion of on/off regulator of fimA
Rejected: Frequency below cutoff threshold.

TATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4541153‑4541255 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tacccACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  <  NC_000913/4541272‑4541174                                                                                                                                                                                                            TATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTA                                                                                                        >  1:2981883/1‑100   TATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACC                                                                                                      >  1:157597/1‑100   TATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACC                                                                                                      >  1:1681683/1‑100   TATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACC                                                                                                      >  2:3282639/1‑100            CGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTC                                                                                             >  2:1298198/1‑100                          CTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTA                                                                               <  1:1990904/100‑1                           TTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTAC                                                                              <  1:1999439/100‑1                               CAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAG                                                                          <  1:1022933/100‑1                               CAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAG                                                                          <  2:303768/100‑1                                AACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGT                                                                         <  2:352633/100‑1                                AACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGT                                                                         <  2:1055954/100‑1                                AACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGT                                                                         <  1:312550/100‑1                                     CGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGA                                                                    <  2:507551/100‑1                                          CCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTC                                                               <  1:1298191/100‑1                                                                        AAAAACTGGTTGAGTATCCGTACTTCGTACCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTA                                 <  1:444612/100‑1                                                                                                    ACCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAG      <  1:2392147/99‑1                                                                                                     CCCACTCGCTCTCAGCTTGCGGGTACGAAGTACGGATACTCTACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAA       <  1:2880882/97‑1                                                                                                     CCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCC    >  1:2834839/1‑100                                                                                                     CCCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCC    <  2:423452/100‑1                                                                                                      CCACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCT   >  2:2933258/1‑100                                                                                                       CACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  >  1:3129926/1‑100                                                                                                       CACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  >  2:1884979/1‑100                                                                                                       CACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  >  2:2693634/1‑100                                                                                                       CACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  >  1:2079992/1‑100                                                                                                       CACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  >  1:1969714/1‑100                                                                                                                                                                                                            TATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4541153‑4541255 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tacccACTCGCTCTCAGCATGCGGGTACGAAGTACGGATACTCAACCAGTTTTTTAAAGCCTGAACTTCTTTATTCAATAGCGGGTGCGTTGTTGAAAAGCCTT  <  NC_000913/4541272‑4541174

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