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

GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC  >  NC_000913/380062‑380156
                                                 |                                             
gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                               <  1:1172584/50‑1 (MQ=255)
gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                               >  1:2338467/1‑50 (MQ=255)
       tCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg                                        >  1:2364513‑M1/1‑43 (MQ=255)
        cAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg                                        >  1:1697561‑M1/1‑42 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga                                      >  1:159338‑M1/1‑41 (MQ=255)
           aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg                                    >  1:752070‑M1/1‑39 (MQ=255)
           aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg                                    <  1:1481071‑M1/50‑12 (MQ=255)
                 cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc                              <  1:2201215‑M1/50‑18 (MQ=255)
                       aGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccc                        >  1:2078592‑M1/1‑27 (MQ=255)
                        gACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccc                       <  1:501216‑M1/50‑25 (MQ=255)
                          cGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccct                     >  1:1476150‑M1/1‑24 (MQ=255)
                           gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt                    >  1:1264525‑M1/1‑23 (MQ=255)
                            ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttg                   <  1:2642827‑M1/50‑29 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggag              <  1:1143684‑M1/50‑34 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggg                <  1:1094846‑M1/48‑32 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg                 >  1:2187642‑M1/1‑17 (MQ=255)
                                    gTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagg           >  1:1577121‑M1/1‑14 (MQ=255)
                                      tGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt         >  1:3813‑M1/1‑12 (MQ=255)
                                       gCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt         >  1:11692‑M1/1‑11 (MQ=255)
                                        cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtta       <  1:175924‑M1/50‑41 (MQ=255)
                                          gACCGTTAcgagatggacggtcccctcgcccccttggggagagggttagg     <  1:1379725‑M1/50‑43 (MQ=255)
                                           aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttaggg    <  1:2089348‑M1/50‑44 (MQ=255)
                                             cGTTAcgagatggacggtcccctcgcccccttggggagagggttagggtg  <  1:339808‑M1/50‑46 (MQ=255)
                                              gTTAcgagatggacggtcccctcgcccccttggggagagggttagggt   >  1:978776‑M1/1‑4 (MQ=255)
                                                 |                                             
GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC  >  NC_000913/380062‑380156

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.