New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)26 (0.660)
+CGAGATGGACG
23/76 0.2 56.0% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 624576 = 53 (1.030)intergenic (‑66/‑309) fepB/entC ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC

ATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCA  >  NC_000913/624527‑624625
                                                 |                                                 
atttccaggtcagacgggctgcaagttgcagaccgttacgagatggacgg                                                   <  1:838008‑M2/1‑1 (MQ=255)
atttccaggtcagacgggctgcaagttgcagaccgttacgagatggacgg                                                   <  1:216859‑M2/1‑1 (MQ=255)
 ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGt                                                  >  1:2203390/1‑50 (MQ=255)
 ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGt                                                  >  1:1762275/1‑50 (MQ=255)
  taataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTc                                                 >  1:2443647/1‑50 (MQ=255)
   aataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTcc                                                >  1:1091144/1‑50 (MQ=255)
   aataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTc                                                 <  1:1423934/49‑1 (MQ=255)
    ccaggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcc                                                <  1:2073064‑M2/4‑1 (MQ=255)
    ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTccc                                               >  1:2752049/1‑50 (MQ=255)
    ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTccc                                               >  1:1341552/1‑50 (MQ=255)
    ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTcc                                                <  1:736671/49‑1 (MQ=255)
     caggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcccc                                              <  1:2885137‑M2/6‑1 (MQ=255)
      aggtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCt                                             >  1:925608‑M2/44‑50 (MQ=255)
      aggtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCt                                             >  1:222645‑M2/44‑50 (MQ=255)
      aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCt                                             >  1:3445382/1‑50 (MQ=255)
       ggtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTc                                            <  1:1355543‑M2/8‑1 (MQ=255)
       aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTc                                            >  1:3791688/1‑50 (MQ=255)
        gtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTCg                                           <  1:1939435‑M2/9‑1 (MQ=255)
        ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                           >  1:150546/1‑50 (MQ=255)
        ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                           >  1:2945872/1‑50 (MQ=255)
         gTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                          >  1:1880963/1‑50 (MQ=255)
           agacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGccc                                        >  1:2130769‑M2/39‑50 (MQ=255)
            gacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGcccc                                       >  1:3023207‑M2/38‑50 (MQ=255)
              gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                     <  1:265130/50‑1 (MQ=255)
              gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                     <  1:2019874/50‑1 (MQ=255)
              cgggctgcaagttgcagaccgttacgagatggacgGTCCCCTcgccccct                                     >  1:1788088‑M2/36‑48 (MQ=255)
               ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCttt                                    >  1:3559808/1‑50 (MQ=255)
                gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg                                   <  1:127225/50‑1 (MQ=255)
                  gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:3127635/1‑50 (MQ=255)
                  gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:2454732/1‑50 (MQ=255)
                  ctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg                                 >  1:3177555‑M2/32‑50 (MQ=255)
                   tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:1024357/1‑49 (MQ=255)
                    aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:2148142/1‑48 (MQ=255)
                    aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:18219/1‑48 (MQ=255)
                    aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:1997995/1‑48 (MQ=255)
                     caagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTgggg                                <  1:1725616‑M2/20‑1 (MQ=255)
                     caagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg                                 >  1:3630949‑M2/29‑47 (MQ=255)
                     caagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga                               <  1:3230867‑M2/21‑1 (MQ=255)
                     aTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                 >  1:1619704/1‑47 (MQ=255)
                     aTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                              <  1:866792/50‑1 (MQ=255)
                      aagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga                               >  1:74448‑M2/28‑48 (MQ=255)
                       gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                            >  1:2816291/1‑50 (MQ=255)
                       gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                            >  1:1100528/1‑50 (MQ=255)
                       agttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGagag                            >  1:910787‑M2/27‑50 (MQ=255)
                         ttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAggg                          <  1:3130230‑M2/26‑1 (MQ=255)
                         cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAggg                          <  1:2390503/50‑1 (MQ=255)
                          aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                         >  1:216033/1‑50 (MQ=255)
                           gcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGt                         >  1:641820‑M2/23‑49 (MQ=255)
                            cagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTa                       <  1:3251438‑M2/29‑1 (MQ=255)
                             ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                      <  1:820074/50‑1 (MQ=255)
                             agaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg                      >  1:3150214‑M2/21‑50 (MQ=255)
                              tCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                     >  1:1138332/1‑50 (MQ=255)
                              tCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                     >  1:1144054/1‑50 (MQ=255)
                               cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                    >  1:3908588/1‑50 (MQ=255)
                               accgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                    >  1:1157006‑M2/19‑50 (MQ=255)
                                gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                    <  1:1153448/49‑1 (MQ=255)
                                gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                    <  1:2506161/49‑1 (MQ=255)
                                gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                   >  1:1767979/1‑50 (MQ=255)
                                ccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                   <  1:1929626‑M2/33‑1 (MQ=255)
                                 cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                   >  1:1689742‑M2/17‑49 (MQ=255)
                                 cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTg                  <  1:1831861‑M2/34‑1 (MQ=255)
                                  gttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGa                 <  1:896410‑M2/35‑1 (MQ=255)
                                  cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                   >  1:3387809/1‑48 (MQ=255)
                                  cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                   >  1:3080360/1‑48 (MQ=255)
                                   aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAg                >  1:2463529/1‑50 (MQ=255)
                                    aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg               <  1:498947/50‑1 (MQ=255)
                                    aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg               <  1:841663/50‑1 (MQ=255)
                                      gcgcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgggg             <  1:2169100/50‑1 (MQ=255)
                                      gcgcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg              >  1:1716741/1‑49 (MQ=255)
                                        gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa           >  1:2003338/1‑50 (MQ=255)
                                        gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa           >  1:631933/1‑50 (MQ=255)
                                         cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa           >  1:2897392/1‑49 (MQ=255)
                                         cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAc          <  1:182518/50‑1 (MQ=255)
                                          acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa           >  1:717078/1‑48 (MQ=255)
                                          acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa           >  1:3087546/1‑48 (MQ=255)
                                           caATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcag        <  1:1804044/50‑1 (MQ=255)
                                            aaTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc       >  1:3354024/1‑50 (MQ=255)
                                               ccGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag    >  1:642258/1‑50 (MQ=255)
                                                cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag    >  1:3504614/1‑49 (MQ=255)
                                                 gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCa  >  1:3166568/1‑50 (MQ=255)
                                                 |                                                 
ATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCA  >  NC_000913/624527‑624625

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