New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4546905 =126 (1.350)6 (0.070) 6/186 NT 4.3% coding (1810/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
?NC_000913 4546951 = 150 (1.730)coding (1856/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Rejected: Frequency below/above cutoff threshold.

AGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4547058‑4546905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGG  >  NC_000913/4546951‑4547095
                                                                                                                                                                                                                                                                                                           
AGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAG                                                                                                                                                         >  3:645619/1‑148
     GTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                    <  2:855277/148‑1
                 TAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG                                                                                                                                                   <  1:85924/137‑1
                 TAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG                                                                                                                                                   >  2:85924/1‑137
                           TGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                    <  4:210195/126‑1
                           TGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                    >  3:210195/1‑126
                                           CGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                    >  4:550919/1‑110
                                           CGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGT                                                                                                                                                    <  3:550919/110‑1
                                                                                                               TGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTC                                                             <  4:319846/129‑1
                                                                                                               TGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTC                                                             >  3:319846/1‑129
                                                                                                                GGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCA                                                                                                                                       >  1:721584/1‑54
                                                                                                                GGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGCCAGCTACAGCA                                                                                                                                       <  2:721584/54‑1
                                                                                                                      TCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGGCAGCTACAGCATGTCACACGATCTCAACGGTCGG                                                                                                                >  1:231625/1‑71
                                                                                                                      TCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTGGCAGCTACAGCATGTCACACGATCTCAACGGTCGG                                                                                                                <  2:231625/71‑1
                                                                                                                                                     CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACA    <  3:955280/148‑1
                                                                                                                                                     CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACA                                                                  <  2:1108184/86‑1
                                                                                                                                                     CAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACA                                                                  >  1:1108184/1‑86
                                                                                                                                                      AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAG   <  1:412957/148‑1
                                                                                                                                                      AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGA                                                                          <  3:693852/77‑1
                                                                                                                                                      AGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGA                                                                          >  4:693852/1‑77
                                                                                                                                                       GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGG  <  3:361640/148‑1
                                                                                                                                                       GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATG                                    >  3:61481/1‑114
                                                                                                                                                       GTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATG                                    <  4:61481/114‑1
                                                                                                                                                                                                                                                                                                           
AGCCGGTTTGCACGCTATAGCTGAGGTTGTTGTCTTCCAGCAACGTACCGTATACACCAGCCAGATTGGTCATCCGACCGTTGAGATCGTGTGACATGCTGTAGCTGGCACTGGCATGTCGCCACTGAGATTTACTGTCAGAACGCAGCCAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4547058‑4546905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagtgCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGG  >  NC_000913/4546951‑4547095

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