New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =48 (0.810)12 (0.230) 10/110 1.3 17.7% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 70 (1.370)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:1961464/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:3701236/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:1443050/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:524567/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:2363762/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                 >  1:157930/1‑65 (MQ=255)
   atatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGt                                                              <  1:4715673‑M2/4‑1 (MQ=255)
   atatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGt                                                              <  1:1438236‑M2/4‑1 (MQ=255)
   aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGt                                                              <  1:2664382/65‑1 (MQ=255)
    tatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGtt                                                             <  1:202140‑M2/5‑1 (MQ=255)
        tttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCAt                                                         >  1:4721313‑M2/57‑65 (MQ=255)
         tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATc                                                        <  1:59097/65‑1 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:4290156/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:2439616/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:2124277/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:2667047/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:3554632/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:287271/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                     >  1:813591/1‑65 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:957350/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:2233188/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:2294258/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:3845914/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:638550/65‑1 (MQ=255)
                aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                 <  1:3754595/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAGATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:201711/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:3061963/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                            <  1:894088/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:2902888/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:293952/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:3655/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:1386180/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:155440/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:3961582/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:1663409/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  1:4728640/65‑1 (MQ=255)
                         gTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtg                                        <  1:4296561/65‑1 (MQ=255)
                         atacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtg                                        <  1:3117237‑M2/26‑1 (MQ=255)
                          tacgtaacttatttatgatacggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3010308‑M2/27‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3523708/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3375337/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3577156/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:4144630/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:101604/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:3193747/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:1514781/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:2201851/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                       <  1:2547799/65‑1 (MQ=255)
                            cgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGc                                     <  1:936359‑M2/29‑1 (MQ=255)
                            cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                     <  1:608386/65‑1 (MQ=255)
                              taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:3409052‑M2/31‑1 (MQ=255)
                              taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:3711189‑M2/31‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:3220639/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2245928/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2495895/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:3715467/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:2815189/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:1664235/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:2723517/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:3765464/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:2935379/65‑1 (MQ=255)
                               aacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                  <  1:4426902‑M2/32‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:2503539/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:349548/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:388043/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:2304032/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                 <  1:960486/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                <  1:1014717/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                <  1:658874/65‑1 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                          >  1:3669473/1‑65 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                          <  1:4788920/65‑1 (MQ=255)
                                        ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                         <  1:1148481/65‑1 (MQ=255)
                                            aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa                     >  1:1668036/1‑65 (MQ=255)
                                               aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc                  >  1:3639742/1‑65 (MQ=255)
                                                gacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGGGCAGTGTTGGCAGTCAAACTCg                 <  1:4513285‑M2/49‑1 (MQ=255)
                                                gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg                 <  1:721850/65‑1 (MQ=255)
                                                gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg                 <  1:1252051/65‑1 (MQ=255)
                                                 tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                >  1:3988311/1‑65 (MQ=255)
                                                 acagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                >  1:1924920‑M2/16‑65 (MQ=255)
                                                     aaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTac            <  1:1980779/65‑1 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca           >  1:1726960/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca           >  1:1702281/1‑65 (MQ=255)
                                                       tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaa          >  1:4010303‑M2/10‑65 (MQ=255)
                                                         gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa        >  1:276089‑M2/8‑65 (MQ=255)
                                                          ccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac       <  1:4648657‑M2/59‑1 (MQ=255)
                                                           cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca      >  1:3560681‑M2/6‑65 (MQ=255)
                                                            cccaATGTTTCATCTTTTGGGTGAAATCTGTGCAGTGTTGTCAGTCAAACTCGTTTacaaaacaa     <  1:2190918‑M2/61‑1 (MQ=255)
                                                            ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa     <  1:3990466‑M2/61‑1 (MQ=255)
                                                            ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa     <  1:997602‑M2/61‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  >  1:4477922‑M2/2‑65 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:816725‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt  <  1:3583142‑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 < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ 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.