New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)35 (1.070)
+CGAGAT
20/86 0.2 69.0% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 4070444 = 36 (0.960)intergenic (‑43/+71) yihR/yihS putative aldose 1‑epimerase YihR/sulfoquinovose isomerase

ACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070395‑4070493
                                                 |                                                 
tcagaatttccaggtcagacgggctgcaagttgcagaccgttacgagatg                                                   <  1:2431160‑M2/1‑1 (MQ=255)
tcagaatttccaggtcagacgggctgcaagttgcagaccgttacgagatg                                                   <  1:1728895‑M2/1‑1 (MQ=255)
 cGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                  <  1:1043001/50‑1 (MQ=255)
  agaatttccaggtcagacgggctgcaagttgcagaccgttacgagatgg                                                  <  1:1756896‑M2/2‑1 (MQ=255)
   gaatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGAc                                                <  1:2313026‑M2/4‑1 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               >  1:914173‑M2/46‑50 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               <  1:1030856‑M2/5‑1 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               >  1:334042‑M2/46‑50 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               <  1:1655749‑M2/5‑1 (MQ=255)
     aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:2703913/1‑50 (MQ=255)
     aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:2245037/1‑50 (MQ=255)
      tttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGt                                             <  1:2083007‑M2/7‑1 (MQ=255)
      tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGt                                             >  1:90609/1‑50 (MQ=255)
        tccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTcc                                           >  1:2517627‑M2/42‑50 (MQ=255)
        tccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTcc                                           >  1:1965446‑M2/42‑50 (MQ=255)
        cAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                           >  1:323134/1‑50 (MQ=255)
        cAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTc                                            <  1:2042700/49‑1 (MQ=255)
         ccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTccc                                          >  1:760965‑M2/41‑50 (MQ=255)
         ccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTccc                                          >  1:2665731‑M2/41‑50 (MQ=255)
         aGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccc                                          <  1:2483171/50‑1 (MQ=255)
          gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                         >  1:2201574/1‑50 (MQ=255)
          caggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTccc                                          <  1:1139357‑M2/10‑1 (MQ=255)
           aggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCt                                        >  1:1772552‑M2/39‑50 (MQ=255)
           aggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCt                                        <  1:1651762‑M2/12‑1 (MQ=255)
            ggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTc                                       >  1:1429710‑M2/38‑50 (MQ=255)
            cTCCCGTGAGGCTGACGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                       <  1:2108307/50‑1 (MQ=255)
             tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                      >  1:423434/1‑50 (MQ=255)
             tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                      >  1:1166122/1‑50 (MQ=255)
             gtcagacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCg                                      >  1:1805783‑M2/37‑50 (MQ=255)
               ccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                    <  1:1968475/50‑1 (MQ=255)
                cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                   <  1:2542098/50‑1 (MQ=255)
                cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                   >  1:499576/1‑50 (MQ=255)
                agacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGccc                                   >  1:2169048‑M2/34‑50 (MQ=255)
                agacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGccc                                   >  1:1219794‑M2/34‑50 (MQ=255)
                 gTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcccc                                  >  1:1893520/1‑50 (MQ=255)
                  tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                 <  1:1081881/50‑1 (MQ=255)
                   cgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTcgccccct                                >  1:365417‑M2/31‑48 (MQ=255)
                   cgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTcgccccct                                >  1:1102969‑M2/31‑48 (MQ=255)
                   cgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTcgccccct                                <  1:1309526‑M2/20‑3 (MQ=255)
                    gggctgcaagttgcagaccgttacgagctGGACGGTCCCCTCGCCCCctt                               >  1:1717677‑M2/30‑50 (MQ=255)
                    gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                               >  1:2588659‑M2/30‑50 (MQ=255)
                    aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                               <  1:1969785/50‑1 (MQ=255)
                    aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                <  1:1912872/49‑1 (MQ=255)
                     ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                              >  1:2220047/1‑50 (MQ=255)
                     ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                              <  1:1622713/50‑1 (MQ=255)
                      gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg                             <  1:2723901‑M2/23‑1 (MQ=255)
                      gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                             >  1:404894/1‑50 (MQ=255)
                        tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1769519/1‑49 (MQ=255)
                          caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGga                          <  1:1595326‑M2/26‑1 (MQ=255)
                          aTGTCAAGCCCAAAAGCTGCACCGGACGGTCCCCTCGCCCCTTTgggg                           <  1:2454039/48‑1 (MQ=25)
                          aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                           <  1:23039/48‑1 (MQ=255)
                          aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                           <  1:2154253/48‑1 (MQ=255)
                             gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                      >  1:2232897‑M2/21‑50 (MQ=255)
                              cAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                     >  1:1335693/1‑50 (MQ=255)
                               tgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:2435449‑M2/19‑50 (MQ=255)
                               tgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:2426103‑M2/19‑50 (MQ=255)
                                gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:2663696‑M2/18‑49 (MQ=255)
                                    cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg               <  1:1786986/50‑1 (MQ=255)
                                    cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg               <  1:1760822/50‑1 (MQ=255)
                                     ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGTGATTGGGTTagggt              >  1:2495778‑M2/13‑49 (MQ=255)
                                     ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt              <  1:1960815‑M2/38‑2 (MQ=255)
                                      cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt              >  1:1063670‑M2/12‑48 (MQ=255)
                                      cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt              >  1:329954‑M2/12‑48 (MQ=255)
                                        aaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag           <  1:770706/50‑1 (MQ=255)
                                        aaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag           <  1:1366854/50‑1 (MQ=255)
                                         aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg          <  1:898096/50‑1 (MQ=255)
                                           agCACCGGTCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg        <  1:376651/49‑1 (MQ=255)
                                             cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa       <  1:626867/49‑1 (MQ=255)
                                             cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa       >  1:79683/1‑49 (MQ=255)
                                              aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa     <  1:2407417/50‑1 (MQ=255)
                                                cGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc   <  1:265932/50‑1 (MQ=255)
                                                 ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc  >  1:274110/1‑50 (MQ=255)
                                                 |                                                 
ACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070395‑4070493

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