New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =62 (1.400)4 (0.100) 4/110 2.1 6.9% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 54 (1.420)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

GCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTG  >  NC_000913/4542932‑4543060
                                                                |                                                                
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                  >  1:813962/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                  >  1:742276/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                  >  1:15937/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                  >  1:488358/1‑65 (MQ=255)
gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                  >  1:3325006/1‑65 (MQ=255)
  caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                                <  1:2836711/65‑1 (MQ=255)
   cacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                               <  1:3038415/64‑1 (MQ=255)
   aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                               <  1:379182/65‑1 (MQ=255)
     caaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                             <  1:1085352/65‑1 (MQ=255)
         tgaaCTGAGAAACTGTTTCCGCTATTCTTTCCTCTATGAGGCAAAATGTCCCCAAATGTTTCATc                                                         <  1:2284864/63‑1 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                      >  1:1539044/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                      >  1:672953/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                      >  1:1303904/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                      >  1:2924667/1‑65 (MQ=255)
            aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                      >  1:1003772/1‑65 (MQ=255)
                    aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                              >  1:3201269/1‑65 (MQ=255)
                    aGTTTATCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                              <  1:176989/65‑1 (MQ=255)
                     gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                             <  1:2654013/65‑1 (MQ=255)
                     ctttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                             <  1:3048541/64‑1 (MQ=255)
                      tttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                            >  1:2136709/1‑65 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCGCAAATGTTTCATCTTTTGGGGGAAAAc                                           <  1:1991234/65‑1 (MQ=255)
                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                           <  1:1723352/65‑1 (MQ=255)
                       aaatacgtaacttagttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAAc                                           <  1:2597021‑M2/24‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                        <  1:3416876/65‑1 (MQ=255)
                          tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                        <  1:1335027/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                    <  1:520493/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                    <  1:3399492/65‑1 (MQ=255)
                              cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                    <  1:1648921/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                   <  1:2228805/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                   <  1:1738620/65‑1 (MQ=255)
                               gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                   <  1:2832364/65‑1 (MQ=255)
                               aacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                   <  1:2150530‑M2/32‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:1569946/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:2702490/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:2078892/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:2394105/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:2337371/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:2212076/65‑1 (MQ=255)
                                aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                  <  1:764788/65‑1 (MQ=255)
                                 tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGTAgtgt                                 <  1:397376/65‑1 (MQ=255)
                                  gCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                >  1:2476837/1‑65 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                           <  1:1373036/65‑1 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                           >  1:3118979/1‑65 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                           <  1:720210/65‑1 (MQ=255)
                                       ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                           >  1:2974625/1‑65 (MQ=255)
                                         tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAgct                         >  1:2630695/1‑63 (MQ=255)
                                              gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                    <  1:968555/65‑1 (MQ=255)
                                              gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                    >  1:2113074/1‑65 (MQ=255)
                                                gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg                  <  1:2496348/65‑1 (MQ=255)
                                                 tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                 >  1:3161354/1‑65 (MQ=255)
                                                 acagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                 >  1:2947063‑M2/16‑65 (MQ=255)
                                                   gaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGttt               <  1:2803262/64‑1 (MQ=255)
                                                   agtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGttt               >  1:3436308‑M2/14‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca            >  1:2704220/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca            >  1:1563158/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca            >  1:229118/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca            >  1:1795957/1‑65 (MQ=255)
                                                      aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTCTGCAGTGTTAGAAATCATACTCGTTTaca            >  1:3274786/1‑65 (MQ=255)
                                                        ggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa          <  1:2345350‑M2/57‑1 (MQ=255)
                                                           cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca       >  1:1362238‑M2/6‑65 (MQ=255)
                                                            ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa      <  1:1757183‑M2/61‑1 (MQ=255)
                                                            ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa      <  1:582861‑M2/61‑1 (MQ=255)
                                                              aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg    <  1:977074‑M2/63‑1 (MQ=255)
                                                              aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg    <  1:876704‑M2/63‑1 (MQ=255)
                                                              aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg    <  1:3192207‑M2/63‑1 (MQ=255)
                                                              aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg    <  1:1578615‑M2/63‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt   <  1:2988647‑M2/64‑1 (MQ=255)
                                                               aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGt   <  1:2716468‑M2/64‑1 (MQ=255)
                                                                aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtg  >  1:133379/1‑65 (MQ=255)
                                                                |                                                                
GCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTG  >  NC_000913/4542932‑4543060

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