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

ACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070395‑4070493
                                                 |                                                 
tcagaatttccaggtcagacgggctgcaagttgcagaccgttacgagatg                                                   >  1:2364513‑M2/50‑50 (MQ=255)
 cagaatttccaggtcagacgggctgcaagttgcagaccgttacgagatg                                                   >  1:1697561‑M2/49‑49 (MQ=255)
 cGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                  <  1:1494020/50‑1 (MQ=255)
 cGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCg                                                   >  1:635788/1‑49 (MQ=255)
 cGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCg                                                   >  1:1608265/1‑49 (MQ=255)
  agaatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGa                                                 >  1:159338‑M2/48‑50 (MQ=255)
   aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                <  1:1334595/50‑1 (MQ=255)
   aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                >  1:348369/1‑50 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               >  1:752070‑M2/46‑50 (MQ=255)
    aatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACg                                               <  1:1481071‑M2/5‑1 (MQ=255)
     aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:875317/1‑50 (MQ=255)
      tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGt                                             <  1:1000739/50‑1 (MQ=255)
      tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGt                                             <  1:2698392/50‑1 (MQ=255)
         aGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccc                                          >  1:130297/1‑50 (MQ=255)
          caggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTcccc                                         <  1:2201215‑M2/11‑1 (MQ=255)
           ccTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCt                                        <  1:2687188/50‑1 (MQ=255)
            cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                       <  1:578305/50‑1 (MQ=255)
             tCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                      >  1:1621904/1‑50 (MQ=255)
                agacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGccc                                   >  1:2078592‑M2/34‑50 (MQ=255)
                 gacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGcccc                                  <  1:501216‑M2/18‑1 (MQ=255)
                 gTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcccc                                  <  1:2327815/50‑1 (MQ=255)
                   cgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTcgccccct                                >  1:1476150‑M2/31‑48 (MQ=255)
                    gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                               >  1:1264525‑M2/30‑50 (MQ=255)
                    aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                               >  1:1991838/1‑50 (MQ=255)
                     ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg                              <  1:2642827‑M2/22‑1 (MQ=255)
                     ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                              <  1:2037300/50‑1 (MQ=255)
                        tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1889791/1‑49 (MQ=255)
                        tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1918929/1‑49 (MQ=255)
                          caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgggg                           <  1:1094846‑M2/25‑1 (MQ=255)
                          caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                            >  1:2187642‑M2/24‑47 (MQ=255)
                          caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGgag                         <  1:1143684‑M2/27‑1 (MQ=255)
                             tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                       <  1:2132906/49‑1 (MQ=255)
                             gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                      >  1:1577121‑M2/21‑50 (MQ=255)
                               tgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:3813‑M2/19‑50 (MQ=255)
                                gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:11692‑M2/18‑49 (MQ=255)
                                 cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTa                  <  1:175924‑M2/34‑1 (MQ=255)
                                 aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                    >  1:2401358/1‑48 (MQ=255)
                                   gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg                <  1:1379725‑M2/36‑1 (MQ=255)
                                    accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg               <  1:2089348‑M2/37‑1 (MQ=255)
                                     aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                >  1:1842193/1‑48 (MQ=255)
                                      cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTtagggtg             <  1:339808‑M2/39‑3 (MQ=255)
                                       gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt              >  1:978776‑M2/11‑47 (MQ=255)
                                         aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg          >  1:292514/1‑50 (MQ=255)
                                         aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg          >  1:1952222/1‑50 (MQ=255)
                                           tGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg         >  1:2615148/1‑49 (MQ=255)
                                            gCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg        >  1:185724/1‑49 (MQ=255)
                                             cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg        <  1:2350952/48‑1 (MQ=255)
                                             cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg        <  1:990530/48‑1 (MQ=255)
                                              cccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa     <  1:525608/49‑1 (MQ=255)
                                              aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa     <  1:1874673/50‑1 (MQ=255)
                                                 ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc  >  1:2356171/1‑50 (MQ=255)
                                                 |                                                 
ACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070395‑4070493

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