New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)20 (0.950)
+CGAGATGGACG
17/76 0.1 63.4% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 624576 = 30 (1.090)intergenic (‑66/‑309) fepB/entC ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC

ATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAG  >  NC_000913/624527‑624623
                                                 |                                               
atttccaggtcagacgggctgcaagttgcagaccgttacgagatggacgg                                                 >  1:1512139‑M2/50‑50 (MQ=255)
 ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGt                                                <  1:551335/50‑1 (MQ=255)
 ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGt                                                <  1:76587/50‑1 (MQ=255)
  taataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTc                                               >  1:18831/1‑50 (MQ=255)
  taataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTc                                               >  1:618504/1‑50 (MQ=255)
   tccaggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcc                                              >  1:1713311‑M2/47‑50 (MQ=255)
   aataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTcc                                              >  1:1271951/1‑50 (MQ=255)
    ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTccc                                             >  1:88351/1‑50 (MQ=255)
     caggtcagacgggctgcaagttgcagaccgttacgagatggacgGTcccc                                            <  1:2329682‑M2/6‑1 (MQ=255)
     caggtcagacgggctgcaagttgcagaccgttacgagatggacgGTccc                                             <  1:1334832‑M2/5‑1 (MQ=255)
      aggtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCt                                           >  1:269085‑M2/44‑50 (MQ=255)
      aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCt                                           >  1:95329/1‑50 (MQ=255)
       ggtcagacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTc                                          <  1:603963‑M2/8‑1 (MQ=255)
       aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCt                                           <  1:2042621/49‑1 (MQ=255)
       aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTc                                          >  1:1904341/1‑50 (MQ=255)
           agacgggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGccc                                      <  1:533001‑M2/12‑1 (MQ=255)
            aaGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGcccc                                     >  1:2245391/1‑50 (MQ=255)
             aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGcccc                                     <  1:723372/49‑1 (MQ=255)
             aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCt                                    >  1:2181376/1‑50 (MQ=255)
             aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCt                                    <  1:209326/50‑1 (MQ=255)
              gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                   >  1:1041964/1‑50 (MQ=255)
              gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                   >  1:1313026/1‑50 (MQ=255)
               gggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCctt                                  <  1:774621‑M2/16‑1 (MQ=255)
               gggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCctt                                  >  1:396766‑M2/35‑50 (MQ=255)
               gggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCctt                                  <  1:1797286‑M2/16‑1 (MQ=255)
               ggCGTAATGACAAATTCGACAAAGCGCACAATCTGTCCCCTCGCCCCttt                                  <  1:1557896/50‑1 (MQ=255)
                ggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTg                                 >  1:1194805‑M2/34‑50 (MQ=255)
                   tgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg                               >  1:619021‑M2/31‑49 (MQ=255)
                     caagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg                               >  1:2338924‑M2/29‑47 (MQ=255)
                     aTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                            <  1:611387/50‑1 (MQ=255)
                      tGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                             >  1:1978066/1‑48 (MQ=255)
                        aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg                         >  1:172915/1‑50 (MQ=255)
                          tgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGt                       >  1:174048‑M2/24‑50 (MQ=255)
                          aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                       >  1:988089/1‑50 (MQ=255)
                          aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                       >  1:1015833/1‑50 (MQ=255)
                            aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                       >  1:604917/1‑48 (MQ=255)
                            aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                       >  1:1126601/1‑48 (MQ=255)
                             agaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg                    >  1:222185‑M2/21‑50 (MQ=255)
                               accgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                  <  1:1806552‑M2/32‑1 (MQ=255)
                                gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                 <  1:139991/50‑1 (MQ=255)
                                ccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                 <  1:89375‑M2/33‑1 (MQ=255)
                                 cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt                 >  1:2200908‑M2/17‑49 (MQ=255)
                                 cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTg                <  1:2059308‑M2/34‑1 (MQ=255)
                                   aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAg              >  1:345902/1‑50 (MQ=255)
                                    tacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg             >  1:165572‑M2/14‑50 (MQ=255)
                                          acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa         >  1:1666028/1‑48 (MQ=255)
                                               ccGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag  >  1:2150411/1‑50 (MQ=255)
                                               ccGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag  >  1:644796/1‑50 (MQ=255)
                                                cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag  >  1:684250/1‑49 (MQ=255)
                                                cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag  >  1:909137/1‑49 (MQ=255)
                                                 |                                               
ATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAG  >  NC_000913/624527‑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.