New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)15 (0.710)
+CGAGATGGACG
14/74 0.2 59.2% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 624576 = 27 (0.990)intergenic (‑66/‑309) fepB/entC ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC

AATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAG  >  NC_000913/624530‑624623
                                              |                                               
tccaggtcagacgggctgcaagttgcagaccgttacgagatggacgGt                                                <  1:1946184‑M2/2‑1 (MQ=255)
aataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGt                                                <  1:785955/48‑1 (MQ=255)
 ccaggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcc                                              >  1:615496‑M2/46‑49 (MQ=255)
 ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTcc                                              >  1:1976968/1‑49 (MQ=255)
   aggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcccc                                            >  1:1859147‑M2/44‑49 (MQ=255)
   aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTcccc                                            >  1:963090/1‑49 (MQ=255)
    aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCt                                           >  1:269063/1‑49 (MQ=255)
     gtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTc                                          <  1:563903‑M2/8‑1 (MQ=255)
       ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGc                                        >  1:1950690/1‑49 (MQ=255)
       cagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGc                                        >  1:1296422‑M2/40‑49 (MQ=255)
         aaGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGccc                                      >  1:1829598/1‑49 (MQ=255)
           gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCt                                    >  1:13756/1‑49 (MQ=255)
            gggctgcaagttgcagaccgttacgagatggacgGTCCCCTcgccccct                                   >  1:1699956‑M2/35‑47 (MQ=255)
            ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                   >  1:868930/1‑49 (MQ=255)
              cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg                                 >  1:630667/1‑49 (MQ=255)
              cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg                                 <  1:193513/49‑1 (MQ=255)
               ctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTgg                                <  1:1997745‑M2/18‑1 (MQ=255)
                 aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                               >  1:1261613/1‑48 (MQ=255)
                  aTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                             <  1:760376/49‑1 (MQ=255)
                   tGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                               >  1:181298/1‑46 (MQ=255)
                    gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                             >  1:437839/1‑47 (MQ=255)
                    agttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga                             >  1:2038162‑M2/27‑47 (MQ=255)
                     aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                           <  1:1926703/48‑1 (MQ=255)
                     aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                          >  1:627239/1‑49 (MQ=255)
                      cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg                         <  1:658320/49‑1 (MQ=255)
                       aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAggg                        >  1:2090153/1‑49 (MQ=255)
                        gcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAggg                        <  1:1556872‑M2/26‑1 (MQ=255)
                          ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                       >  1:1648309/1‑47 (MQ=255)
                          ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                     <  1:867276/49‑1 (MQ=255)
                           gaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg                    >  1:2070648‑M2/20‑49 (MQ=255)
                             gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                  >  1:1895895/1‑49 (MQ=255)
                             gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                  <  1:1016127/49‑1 (MQ=255)
                             gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                   <  1:1581616/48‑1 (MQ=255)
                              cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                 <  1:2105891‑M2/33‑1 (MQ=255)
                               gttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                 >  1:1890420‑M2/16‑48 (MQ=255)
                                ttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGa               <  1:1687698‑M2/35‑1 (MQ=255)
                                aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                 >  1:1952532/1‑47 (MQ=255)
                                 tacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg              >  1:687980‑M2/14‑49 (MQ=255)
                                  acgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg              <  1:1744505‑M2/36‑1 (MQ=255)
                                  aGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg             <  1:753285/49‑1 (MQ=255)
                                    gagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgggg           >  1:509892‑M2/11‑49 (MQ=255)
                                    gagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgggg           >  1:764998‑M2/11‑49 (MQ=255)
                                        caATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAca       <  1:1781682/49‑1 (MQ=255)
                                             cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag  >  1:2103254/1‑49 (MQ=255)
                                              |                                               
AATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAG  >  NC_000913/624530‑624623

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.