New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =30 (0.560)3 (0.060) 3/116 3.1 7.1% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 53 (1.150)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAC  >  NC_000913/4542930‑4543063
                                                                  |                                                                   
aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTGTATGAGTCAAAATGGCCCCAAAt                                                                    <  1:2070849/68‑1 (MQ=255)
aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAt                                                                    >  1:3164949/1‑68 (MQ=255)
 aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                                   <  1:2095011/68‑1 (MQ=255)
 aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                                   <  1:2013778/68‑1 (MQ=255)
  gCCAACAAATAAACTGAAAAAGTTTGTCCGTGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                                  >  1:1503855/1‑68 (MQ=255)
  gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATgc                                                                  >  1:3795617/1‑67 (MQ=255)
    caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                <  1:419412/68‑1 (MQ=255)
              aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                      >  1:1013318/1‑68 (MQ=255)
              aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                      >  1:1184615/1‑68 (MQ=255)
              aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                      >  1:2091627/1‑68 (MQ=255)
               cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                     <  1:2604424/68‑1 (MQ=255)
               cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                     <  1:2193803/68‑1 (MQ=255)
                    aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGNGGGaaa                                                <  1:2448500/68‑1 (MQ=255)
                    aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                <  1:990074/68‑1 (MQ=255)
                      aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                              <  1:3416627/68‑1 (MQ=255)
                      aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                              <  1:2969789/68‑1 (MQ=255)
                      aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                              >  1:80519/1‑68 (MQ=255)
                             ccacgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:3613708/65‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:1863350/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:114594/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:2683272/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:270924/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:2105007/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:423103/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:2789763/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:3916616/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:1936585/68‑1 (MQ=255)
                             ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                       <  1:2885333/68‑1 (MQ=255)
                              cgtgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                      <  1:3417466/65‑1 (MQ=255)
                              cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                      <  1:641280/68‑1 (MQ=255)
                               gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGATAACTGTGCAgtg                                     <  1:2548461/68‑1 (MQ=255)
                               gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                     <  1:2259694/68‑1 (MQ=255)
                               gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                     <  1:2282586/68‑1 (MQ=255)
                                cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                    <  1:2123553/68‑1 (MQ=255)
                                  aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                                  >  1:3750517/1‑68 (MQ=255)
                                       ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                             <  1:1482351/68‑1 (MQ=255)
                                       ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                             <  1:2616649/68‑1 (MQ=255)
                                         ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                           >  1:3695602/1‑68 (MQ=255)
                                         ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                           >  1:1919395/1‑68 (MQ=255)
                                         ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                           >  1:1927255/1‑68 (MQ=255)
                                         ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                           >  1:141623/1‑68 (MQ=255)
                                         ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                           >  1:4064877/1‑68 (MQ=255)
                                          ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaa                          <  1:3863941/68‑1 (MQ=255)
                                           tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa                         >  1:1582457/1‑68 (MQ=255)
                                           tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa                         <  1:523356/68‑1 (MQ=255)
                                             tATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                       <  1:2878721/68‑1 (MQ=255)
                                              aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc                      >  1:4046640/1‑68 (MQ=255)
                                               tGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg                     <  1:3930531/68‑1 (MQ=255)
                                                gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                    >  1:3037071/1‑68 (MQ=255)
                                                   tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                 >  1:1707886/1‑68 (MQ=255)
                                                        ttGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa            >  1:2218168/2‑68 (MQ=255)
                                                        aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa            >  1:1448789/1‑68 (MQ=255)
                                                        aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa            >  1:533994/1‑68 (MQ=255)
                                                        aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa            >  1:966624/1‑68 (MQ=255)
                                                         tGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac           >  1:1604454/1‑68 (MQ=255)
                                                           gCCTCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa         <  1:2955839/68‑1 (MQ=255)
                                                           gCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa         <  1:3040106/68‑1 (MQ=255)
                                                             cccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg       <  1:1750052/68‑1 (MQ=255)
                                                             cccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg       <  1:207759/68‑1 (MQ=255)
                                                              ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt      <  1:3521516/68‑1 (MQ=255)
                                                              ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt      <  1:1070642/68‑1 (MQ=255)
                                                              ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt      <  1:1865448/68‑1 (MQ=255)
                                                               cAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtg     <  1:1452792/68‑1 (MQ=255)
                                                               cAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtg     <  1:3832742/68‑1 (MQ=255)
                                                                aaaTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt    <  1:1630007/68‑1 (MQ=255)
                                                                aaaTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt    <  1:2274622/68‑1 (MQ=255)
                                                                aaaTGTGTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACNAAgtgt    <  1:2458611/68‑1 (MQ=255)
                                                                  aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc  >  1:2142749/1‑68 (MQ=255)
                                                                  |                                                                   
AAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAC  >  NC_000913/4542930‑4543063

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