New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4546905 =	217 (1.370)	6 (0.040)	6/252	NT	2.6%	coding (1810/2637 nt)	fimD	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
	?	NC_000913	4546951 =	237 (1.580)					coding (1856/2637 nt)	fimD	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Rejected: Frequency below/above cutoff threshold.

CGTAATAGAGCTGCTTAATATCATCGCTATGGCTGTAACCGATATTGGCATTGCCGTAACCACCGCGATAATTCAGCGTGGCGTAGCCTGTACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4547181‑4546905 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTA  >  NC_000913/4546951‑4547219                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    CGTAATAGAGCTGCTTAATATCATCGCTATGGCTGTAACCGATATTGGCATTGCCGTAACCACCGCGATAATTCAGCGTGGCGTAGCCTGTACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGCGATTTACTGTCAGAACGCAGCCAG                                                                                                                                                                                                                                                                                 >  1:147181/1‑275                                                                CGCGATAATTCAGCGTGGCGTAGCCTGTACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                >  4:307990/1‑213                                                                                      GCCTGTACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAG                                                                                                                                                                                                                                                                                 <  2:561931/190‑1                                                                                           TACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAG                                                                                                                                                                                                                                                                                     >  4:97857/1‑181                                                                                               TCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAG                                                                                                                                                                                                                                                                                 <  2:97858/181‑1                                                                                                                            AGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAG                                                                                                                                                                                                                                                                                     >  7:645619/1‑148                                                                                                                                 GTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                <  6:855277/148‑1                                                                                                                                             TAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG                                                                                                                                                                                                                                                                               >  6:85924/1‑137                                                                                                                                             TAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG                                                                                                                                                                                                                                                                               <  5:85924/137‑1                                                                                                                                                       TGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                >  7:210195/1‑126                                                                                                                                                       TGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                <  8:210195/126‑1                                                                                                                                                                       CGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                <  7:550919/110‑1                                                                                                                                                                       CGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                                                                                                                                                >  8:550919/1‑110                                                                                                                                                                                                                                           TGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTC                                                                                                                                                                                         >  7:319846/1‑129                                                                                                                                                                                                                                           TGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTC                                                                                                                                                                                         <  8:319846/129‑1                                                                                                                                                                                                                                            GGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCA                                                                                                                                                                                                                                                                   <  6:721584/54‑1                                                                                                                                                                                                                                            GGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCA                                                                                                                                                                                                                                                                   >  5:721584/1‑54                                                                                                                                                                                                                                                  TCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGGCAGCTACAGCATGTCACACGATCTCAACGGTCGG                                                                                                                                                                                                                                            <  6:231625/71‑1                                                                                                                                                                                                                                                  TCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGGCAGCTACAGCATGTCACACGATCTCAACGGTCGG                                                                                                                                                                                                                                            >  5:231625/1‑71                                                                                                                                                                                                                                                                                CCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGGGGTTACGGCAATGCCAATATCGGGTACAGCCATAGCGATGATATTAAGCGGCTCTATTACGGAGTCAACGGCGGGGTACCGGCTCATGCCAcaggggta  >  4:56631/1‑267                                                                                                                                                                                                                                                                                 CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACA                                                                                                                                <  7:955280/148‑1                                                                                                                                                                                                                                                                                 CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACA                                                                                                                                                                                              <  6:1108184/86‑1                                                                                                                                                                                                                                                                                 CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACA                                                                                                                                                                                              >  5:1108184/1‑86                                                                                                                                                                                                                                                                                  AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAG                                                                                                                               <  5:412957/148‑1                                                                                                                                                                                                                                                                                  AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGA                                                                                                                                                                                                      >  8:693852/1‑77                                                                                                                                                                                                                                                                                  AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGA                                                                                                                                                                                                      <  7:693852/77‑1                                                                                                                                                                                                                                                                                   GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGG                                                                                                                              <  7:361640/148‑1                                                                                                                                                                                                                                                                                   GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATG                                                                                                                                                                >  7:61481/1‑114                                                                                                                                                                                                                                                                                   GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATG                                                                                                                                                                <  8:61481/114‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    CGTAATAGAGCTGCTTAATATCATCGCTATGGCTGTAACCGATATTGGCATTGCCGTAACCACCGCGATAATTCAGCGTGGCGTAGCCTGTACTTCCGCTATTTCCATCGCCTCCCCCGGCATAGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4547181‑4546905 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTA  >  NC_000913/4546951‑4547219

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