New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4540250 =71 (0.450)4 (0.030) 4/132 14.2 5.7% intergenic (+49/‑433) fimE/fimA FimE/major type 1 subunit fimbrin
?NC_000913 4540564 = 71 (0.510)intergenic (+363/‑119) fimE/fimA FimE/major type 1 subunit fimbrin
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAC  >  NC_000913/4540489‑4540639
                                                                           |                                                                           
tAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                             >  2:1916007/1‑76 (MQ=255)
 aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAt                                                                            <  1:2286906/76‑1 (MQ=255)
  gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                                           <  1:992884/76‑1 (MQ=255)
  gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                                           <  1:2606703/76‑1 (MQ=255)
   aaTATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                                          >  1:673787/1‑76 (MQ=255)
    atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtt                                                                         <  2:1939873/76‑1 (MQ=255)
    atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtt                                                                         >  2:370323/1‑76 (MQ=255)
      atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtctt                                                                       >  2:3109076/1‑73 (MQ=255)
      atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                       <  2:4843118/76‑1 (MQ=255)
      atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGACAAAATGGCCTCGAATGTTTc                                                                       >  1:2985453/1‑76 (MQ=255)
       ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCa                                                                      >  1:4211721/1‑76 (MQ=255)
       ttAAGCCAACAAATAAACTGAAAAAGATTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCa                                                                      >  2:1635902/1‑76 (MQ=255)
         aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATc                                                                    <  1:1125582/76‑1 (MQ=255)
           gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtt                                                                  <  2:1306617/76‑1 (MQ=255)
           gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtt                                                                  >  2:1861533/1‑76 (MQ=255)
             caccaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                <  2:2114302/73‑1 (MQ=255)
                 aaaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgggg                                                            >  2:4342209/1‑76 (MQ=255)
                 aaaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgggg                                                            >  2:2968297/1‑76 (MQ=255)
                  aaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                           >  1:4711115/1‑76 (MQ=255)
                   aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGa                                                          >  2:2703166/1‑76 (MQ=255)
                   aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGa                                                          <  2:4036537/76‑1 (MQ=255)
                    tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                             <  2:3855020/72‑1 (MQ=255)
                    tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                             >  1:3855021/1‑72 (MQ=255)
                    tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                                         <  1:1868295/76‑1 (MQ=255)
                     aaaCTGAAAAAGTTTGTCCGCGATGCTGTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                        <  2:5305715/76‑1 (MQ=255)
                      aaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                       >  2:2599779/1‑76 (MQ=255)
                      aaCTGAAAAAGTTTGTCCGCGATGCTTCCCTCTATGAGTCAAAATGGCCCCACATGTTTAATCTTTTGCGGGacaa                                                       >  2:2610617/1‑76 (MQ=255)
                        cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACt                                                     >  2:1309700/1‑76 (MQ=255)
                             aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGGAAACTGCGCa                                                >  1:1511765/1‑76 (MQ=255)
                               aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                              >  1:4954202/1‑76 (MQ=255)
                               aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                              <  1:1916007/76‑1 (MQ=255)
                                gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGACAACTGTGCAgtg                                             >  1:4067194/1‑76 (MQ=255)
                                gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                             >  2:735709/1‑76 (MQ=255)
                                gTTTGTCCGCGATGCGTGCCTCGATGAGTCAAAATGGCCGCAAATGTTTCATCTTTTGGGGGACAAGTGTGCAgtg                                             >  1:4092307/1‑76 (MQ=255)
                                 tttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                            >  2:159300/1‑76 (MQ=255)
                                 tttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                            >  2:763891/1‑76 (MQ=255)
                                  ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                           >  1:4616995/1‑76 (MQ=255)
                                   tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                                          >  1:1035933/1‑76 (MQ=255)
                                    gTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                         <  1:1777194/76‑1 (MQ=255)
                                    gTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                         >  1:1428897/1‑76 (MQ=255)
                                     tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGc                                        >  2:5403207/1‑76 (MQ=255)
                                     tCCGCGATGATTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgga                                        >  2:583059/1‑75 (MQ=255)
                                      ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGNGAAAACTGTGCAGTGTTGGCa                                       <  1:3292479/76‑1 (MQ=255)
                                      ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                       <  2:1199078/76‑1 (MQ=255)
                                      ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                       >  2:5834657/1‑76 (MQ=255)
                                      ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                       >  1:3302695/1‑76 (MQ=255)
                                       cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGTCGGAAAACTGTGCAGTGTTGGCAg                                      >  1:59334/1‑76 (MQ=255)
                                       cgcgAGGATTTCCTCTATGAGTCAAAATGGCCAAAAATGGTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                                      >  1:49972/1‑76 (MQ=255)
                                        gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                     >  2:1761509/1‑76 (MQ=255)
                                        gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                     >  1:3660680/1‑76 (MQ=255)
                                        gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                     >  1:3725833/1‑76 (MQ=255)
                                        gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                     <  1:615188/76‑1 (MQ=255)
                                         cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                                    <  1:120214/76‑1 (MQ=255)
                                          gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                                   >  2:596928/1‑76 (MQ=255)
                                          gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                                   >  1:4931815/1‑76 (MQ=255)
                                          gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                                   >  2:511983/1‑76 (MQ=255)
                                              cTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                               >  1:523377/1‑76 (MQ=255)
                                               tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc                              <  2:4616994/76‑1 (MQ=255)
                                               tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc                              >  2:5663914/1‑76 (MQ=255)
                                                 tCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                            >  1:2644294/1‑76 (MQ=255)
                                                 tCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGGGGTAGCAGTCAAACTCGt                            >  1:4326933/1‑76 (MQ=255)
                                                      tATGAGTCAAAATAGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCATTGTTGGCCGTCAAACTCGTTTaca                       <  1:1635902/76‑1 (MQ=255)
                                                        tGAGTAAAAATGGCCCCAAATGTTTCATCGTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTagaaa                     >  1:1393683/1‑76 (MQ=255)
                                                         gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                    <  2:5051560/76‑1 (MQ=255)
                                                          aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac                   >  2:294047/1‑76 (MQ=255)
                                                           gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca                  >  2:4029659/1‑76 (MQ=255)
                                                           gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca                  <  2:2269539/76‑1 (MQ=255)
                                                             cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGGGTTGGCAGTCAAACTCGTGTaccaaacaaa                >  2:1716889/1‑76 (MQ=255)
                                                              aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACCAAg               >  2:408224/1‑76 (MQ=255)
                                                              aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg               <  1:4342210/76‑1 (MQ=255)
                                                                 aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGCCAAACTCGTTTACAAAACAAAgtga            >  2:90236/1‑75 (MQ=255)
                                                                   ggCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc          >  2:3957760/1‑76 (MQ=255)
                                                                   ggCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc          >  1:1903340/1‑76 (MQ=255)
                                                                    gCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACa         >  1:59890/1‑76 (MQ=255)
                                                                      cccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGa       <  1:3234613/76‑1 (MQ=255)
                                                                       ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa      >  2:4924128/1‑76 (MQ=255)
                                                                       ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa      >  1:4632775/1‑76 (MQ=255)
                                                                        cAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc     >  2:4526535/1‑76 (MQ=255)
                                                                        cAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc     <  1:2968297/76‑1 (MQ=255)
                                                                          aaTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa   <  2:3660680/76‑1 (MQ=255)
                                                                           aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc  >  1:5066169/1‑76 (MQ=255)
                                                                           |                                                                           
TAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAC  >  NC_000913/4540489‑4540639

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