New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454572348 (0.930)3 (0.060) 3/246 NT 6.1% coding (628/2637 nt) fimD type I fimbriae usher protein
?NC_000913 = 4545737 47 (0.970)coding (642/2637 nt) fimD type I fimbriae usher protein

CCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4545575‑4545723
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atatttaaCCCACTCTGTAGGTTTAAATATGCATAATGGCTGTTACCCCCAATCCGATTCTGTACACTATTTCCGCTGAAATTATAATTGAGCAATCCGGCATTAATACCGGGATCCCATAACTCAGGAGGAATATAACCACGCGCGC  <  NC_000913/4545737‑4545598
                                                                                                                                                                                                                                                                                                 
ctCTCAGGTATTTATGAGTAATCGCGCGAGGGGGTATATTCCTCCTGAGTTATGGGATGCCGCTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCAT                                                                                                                                                    <  1:512658/141‑1
CCCTCAGTGATTTATGAGTAATCGCGTGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCAT                                                                                                                                                    <  2:959529/143‑1
CCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCAT                                                                                                                                                    <  1:284862/143‑1
  CTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGACGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGCGGTAACAGACATTATGCATAT                                                                                                                                                  >  2:149216/1‑143
    CAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGGATTAATGCAGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAAAATAGGATTGGGGTTAAAAGCCATTAAGCATATTc                                                                                                                                                >  2:384386/1‑142
     AGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGAGGTAACCGCCAATATGCATATATA                                                                                                                                               >  2:272643/1‑143
      GGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAA                                                                                                                                              <  2:411250/143‑1
                         GCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAACCCACTCTGTAGGTTTGAA                                                                                                                           <  2:971866/143‑1
                                                                 TTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAACCCACTCTGTAGGTTTAAATAT                                                                                                                        <  1:875076/106‑1
                                                                 TTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAACCCACTCTGTAGGTTTAAATAT                                                                                                                        >  2:875076/1‑106
                                                                                                                                                TTTCACCCCCTCTGTATCTTTAAATATGCATAATGTCTGTTACCCCCAATCCTATTCTCTACACTATTTCCGCTGAAATTATAATTGAGCAATCCGACATTAATACCGGGATCCCATAACTCAGGAGGAATATAACCACGCGC    <  2:546349/143‑1
                                                                                                                                                  TAACCCACTCTGTAGGTTTAAATATGCATAATGGCTGTTACCCCCAATCCGATTCTGTACACTATTTCCGCTGAAATTATAATTGAGCAATCCGGCATTAATACCGGGATCCCATAACTCAGGAGGAATATAAACACGCGCGC  >  2:150741/1‑143
                                                                                                                                                                                                                                                                                                 
CCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4545575‑4545723
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atatttaaCCCACTCTGTAGGTTTAAATATGCATAATGGCTGTTACCCCCAATCCGATTCTGTACACTATTTCCGCTGAAATTATAATTGAGCAATCCGGCATTAATACCGGGATCCCATAACTCAGGAGGAATATAACCACGCGCGC  <  NC_000913/4545737‑4545598

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