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.

GCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGT  >  NC_000913/4542932‑4543059
                                                                |                                                               
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:901244/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:197425/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:1758731/1‑65 (MQ=255)
  caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                               <  1:2886321/65‑1 (MQ=255)
  caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                               <  1:2412119/65‑1 (MQ=255)
  caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                               >  1:34033/1‑65 (MQ=255)
  caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                               >  1:608663/1‑65 (MQ=255)
  aatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATg                                                               <  1:3179719‑M2/3‑1 (MQ=255)
   aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                              <  1:2149341/65‑1 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:3050439/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:1162547/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:1230496/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:2361677/1‑65 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCTAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:1138081/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:3087882/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:1667454/65‑1 (MQ=255)
                    aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                             >  1:3246310/1‑65 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:3075266/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:79535/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:2998603/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaacc                                          <  1:738729/65‑3 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:2310905/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:1052183/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTGTCATCTTTTGGGGGAAAAc                                          <  1:33186/65‑1 (MQ=255)
                          tacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:2719883‑M2/27‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3168512/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:507346/65‑1 (MQ=255)
                           ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                      >  1:395624/1‑65 (MQ=255)
                              taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2052137‑M2/31‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:1505927/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2802314/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:1237726/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2452677/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:1949590/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:2099852/65‑1 (MQ=255)
                               aacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:3399464‑M2/32‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:291687/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:1974685/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:1586541/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                <  1:3699720/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                <  1:230891/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                <  1:118713/65‑1 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                          <  1:720566/65‑1 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                          >  1:704741/1‑65 (MQ=255)
                                        ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                         <  1:2747133/65‑1 (MQ=255)
                                        ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                         <  1:275806/65‑1 (MQ=255)
                                              gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                   >  1:3938716/1‑65 (MQ=255)
                                               aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc                  >  1:2831080/1‑65 (MQ=255)
                                                 tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGCTGGCAGTCAAACTcct                >  1:3022297/1‑63 (MQ=255)
                                                      ttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca           >  1:2367057‑M2/11‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca           >  1:2406500/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca           >  1:2377745/1‑65 (MQ=255)
                                                       tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaa          >  1:1347121‑M2/10‑65 (MQ=255)
                                                        ggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa         <  1:3014668‑M2/57‑1 (MQ=255)
                                                        ggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa         <  1:1550678‑M2/57‑1 (MQ=255)
                                                          ccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac       <  1:1355538‑M2/59‑1 (MQ=255)
                                                           cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca      >  1:3476648‑M2/6‑65 (MQ=255)
                                                           cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTGGGCAGTCAAACTCGATTacaaaaca      >  1:259877‑M2/6‑65 (MQ=255)
                                                            ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa     <  1:113398‑M2/61‑1 (MQ=255)
                                                            ccaaATGTTTCAGCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa     <  1:730721‑M2/61‑1 (MQ=255)
                                                              aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg   <  1:1139980‑M2/63‑1 (MQ=255)
                                                               aATTTTTCATCTTGTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:33512‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:2377085‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:371806‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:3479670‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:1199877‑M2/64‑1 (MQ=255)
                                                                |                                                               
GCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGT  >  NC_000913/4542932‑4543059

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