New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269052 (1.020)4 (0.090) 4/110 2.3 8.0% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 47 (1.080)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

ATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACA  >  NC_000913/4542925‑4543050
                                                             |                                                                
atatTAAGCCGACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:2141121‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGCGTCAAAAtgg                                                               <  1:209133‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:3393916‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:3707667‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:3579907‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:658507‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               >  1:726930‑M2/1‑62 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:1466873‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:138185‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:1294578‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:2434233‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:3220743‑M2/65‑4 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                               <  1:3117478‑M2/65‑4 (MQ=255)
  atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                             <  1:2282263‑M2/65‑6 (MQ=255)
  atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                             <  1:3210042‑M2/65‑6 (MQ=255)
  atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                             <  1:3483101‑M2/65‑6 (MQ=255)
  atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                             <  1:3660579‑M2/65‑6 (MQ=255)
    tAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccac                                                           <  1:3351030‑M2/65‑8 (MQ=255)
       gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                        >  1:1758731/1‑65 (MQ=255)
       gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                        >  1:197425/1‑65 (MQ=255)
       gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                        >  1:901244/1‑65 (MQ=255)
         caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                      >  1:608663/1‑65 (MQ=255)
         caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                      <  1:2412119/65‑1 (MQ=255)
         caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                      <  1:2886321/65‑1 (MQ=255)
         caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                      >  1:34033/1‑65 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                     <  1:2149341/65‑1 (MQ=255)
                   aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                            >  1:2361677/1‑65 (MQ=255)
                   aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                            >  1:1230496/1‑65 (MQ=255)
                   aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                            >  1:1162547/1‑65 (MQ=255)
                   aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                            >  1:3050439/1‑65 (MQ=255)
                       aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCTAAATGGCCCCAAATGTTTCATCTTTTggg                                        <  1:1138081/65‑1 (MQ=255)
                       aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                        <  1:1667454/65‑1 (MQ=255)
                       aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                        <  1:3087882/65‑1 (MQ=255)
                           aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgttt                                    <  1:2657675‑M2/65‑31 (MQ=255)
                           aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                    >  1:3246310/1‑65 (MQ=255)
                            gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                   <  1:79535/65‑1 (MQ=255)
                            gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                   <  1:3075266/65‑1 (MQ=255)
                            gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                   <  1:2998603/65‑1 (MQ=255)
                              ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaacc                                 <  1:738729/65‑3 (MQ=255)
                              ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                 <  1:1052183/65‑1 (MQ=255)
                              ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                 <  1:2310905/65‑1 (MQ=255)
                              ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTGTCATCTTTTGGGGGAAAAc                                 <  1:33186/65‑1 (MQ=255)
                                 tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                              <  1:3168512/65‑1 (MQ=255)
                                 tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                              <  1:507346/65‑1 (MQ=255)
                                  ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                             >  1:395624/1‑65 (MQ=255)
                                     cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                          <  1:2802314/65‑1 (MQ=255)
                                     cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                          <  1:1237726/65‑1 (MQ=255)
                                     cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                          <  1:1505927/65‑1 (MQ=255)
                                     cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                          <  1:2452677/65‑1 (MQ=255)
                                      gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                         <  1:2099852/65‑1 (MQ=255)
                                      gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                         <  1:1949590/65‑1 (MQ=255)
                                       aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                        <  1:1974685/65‑1 (MQ=255)
                                       aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                        <  1:1586541/65‑1 (MQ=255)
                                       aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                        <  1:291687/65‑1 (MQ=255)
                                        tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                       <  1:230891/65‑1 (MQ=255)
                                        tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                       <  1:118713/65‑1 (MQ=255)
                                        tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                       <  1:3699720/65‑1 (MQ=255)
                                            ttCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatca                   >  1:1016866‑M2/1‑18 (MQ=255)
                                              ccTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaat                 <  1:2250630‑M2/65‑50 (MQ=255)
                                              ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                 >  1:704741/1‑65 (MQ=255)
                                              ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                 <  1:720566/65‑1 (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)
                                                           aaatggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaa    >  1:2541938‑M2/1‑3 (MQ=255)
                                                             aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca  >  1:2377745/1‑65 (MQ=255)
                                                             aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca  >  1:2406500/1‑65 (MQ=255)
                                                             |                                                                
ATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACA  >  NC_000913/4542925‑4543050

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