New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =52 (1.020)5 (0.110) 5/110 2.0 9.8% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 47 (1.080)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATT  >  NC_000913/4542618‑4542743
                                                                |                                                             
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccacaagatgaaaaatt                                                               <  1:33512‑M1/1‑1 (MQ=255)
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                               <  1:1199877‑M1/1‑1 (MQ=255)
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                               <  1:2377085‑M1/1‑1 (MQ=255)
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                               <  1:371806‑M1/1‑1 (MQ=255)
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                               <  1:3479670‑M1/1‑1 (MQ=255)
aagagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                               <  1:707754‑M1/65‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                              <  1:1139980‑M1/2‑1 (MQ=255)
 agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                              >  1:3643072‑M1/1‑65 (MQ=255)
   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagctgaaacatTTgg                                                            <  1:730721‑M1/4‑1 (MQ=255)
   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                            <  1:113398‑M1/4‑1 (MQ=255)
    tgttttgtaatcgagtttgactgcccacactgcacagttttcccccaaaagatgaaacatTTggg                                                           >  1:259877‑M1/61‑65 (MQ=255)
    tgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTggg                                                           >  1:3476648‑M1/61‑65 (MQ=255)
    gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGG                                                           <  1:1953480‑M1/65‑1 (MQ=255)
     gttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                          <  1:1355538‑M1/6‑1 (MQ=255)
      agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                         >  1:3215608‑M1/1‑65 (MQ=255)
      agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                         >  1:1636969‑M1/1‑65 (MQ=255)
       tttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGcc                                                        <  1:3014668‑M1/8‑1 (MQ=255)
       tttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGcc                                                        <  1:1550678‑M1/8‑1 (MQ=255)
        ttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCa                                                       >  1:1347121‑M1/57‑65 (MQ=255)
         tgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                      >  1:2367057‑M1/56‑65 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                      >  1:96750/1‑65 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                      >  1:2427962/1‑65 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                      <  1:3758968/65‑1 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                      <  1:1475696/65‑1 (MQ=255)
           tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttt                                                    >  1:2541938‑M1/1‑62 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                   <  1:2970056/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                   <  1:3858812/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                   <  1:1755796/65‑1 (MQ=255)
              gATTTACCTTATTGATAATAAAGTTAAAAAAACAAATAAAAACAAGACAATTGGGGCCAAACTGt                                                 <  1:1954403/65‑1 (MQ=255)
              gATTTACCTTATTGATAATAAAGTTAAAAAAACAAATAAAAAAAAGCCAATTGGGGCCAAACTGt                                                 <  1:1566552/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                 <  1:1949664/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                 <  1:3483624/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                 <  1:1237463/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                 <  1:558212/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATGGGGGCCAAACTGt                                                 <  1:2750918/65‑1 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAACCTGTcc                                               >  1:2296765/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                               >  1:3319662/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                               >  1:1363059/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                               >  1:3752379/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                               >  1:2163222/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                          >  1:2681186/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                          >  1:28721/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                          <  1:890012/65‑1 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                          <  1:1856865/65‑1 (MQ=255)
                      ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCa                                         >  1:1350862/1‑65 (MQ=255)
                       tATTGATAATAAAGTTAAAAAAACAAATAAATACGAGACAATTGGGGCCAAACTGTCCATATCat                                        <  1:618863/65‑1 (MQ=255)
                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagagg                                       <  1:2250630‑M1/65‑17 (MQ=255)
                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                       <  1:3428048/65‑1 (MQ=255)
                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaa                                     >  1:1016866‑M1/1‑47 (MQ=255)
                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                     >  1:212464/1‑65 (MQ=255)
                            ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                   <  1:1238773/65‑1 (MQ=255)
                            ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                   >  1:2736354/1‑65 (MQ=255)
                             taataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                  <  1:3507723/65‑1 (MQ=255)
                              aataaAGTTAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                <  1:1910215/65‑1 (MQ=255)
                               ataaAGTTAAAAAAACAAATAAATGCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                <  1:3818818/65‑1 (MQ=255)
                                actgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGtt                               <  1:3399464‑M1/33‑1 (MQ=255)
                                 ctgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTa                              <  1:2052137‑M1/34‑1 (MQ=255)
                                   aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACg                            >  1:1763964/1‑65 (MQ=255)
                                    gTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                           >  1:1253999/1‑65 (MQ=255)
                                     acagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                          <  1:2719883‑M1/38‑1 (MQ=255)
                                      tAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAt                         <  1:1203030/65‑1 (MQ=255)
                                           aaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaact                    <  1:2657675‑M1/65‑36 (MQ=255)
                                                 ataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAg              <  1:1474271/65‑1 (MQ=255)
                                                  taaataCAAGACAATTGGGGCCAAACTGTCCATACCATAAATAAGTTACGTATTTTTTCTCAAGc             <  1:3539173/65‑1 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa          >  1:2458113/1‑65 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa          >  1:3430705/1‑65 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa          >  1:128693/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaat        >  1:3644289/1‑64 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa        >  1:3175387/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa        >  1:3510959/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa        >  1:2928993/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa        >  1:1327023/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCTATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:2564772/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:2813939/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:51743/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:1710641/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:3431668/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:813804/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat   >  1:1416009/1‑65 (MQ=255)
                                                             catTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt  <  1:3179719‑M1/62‑1 (MQ=255)
                                                                |                                                             
AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATT  >  NC_000913/4542618‑4542743

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