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

GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/380062‑380154
                                                 |                                           
gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                             >  1:2244016/1‑50 (MQ=255)
    gAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgag                                         <  1:358511‑M1/50‑5 (MQ=255)
    gAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                            >  1:2207321‑M1/1‑46 (MQ=255)
     aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgaga                                        <  1:6043258‑M1/50‑6 (MQ=255)
     aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                            >  1:2214405‑M1/1‑45 (MQ=255)
      aTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagat                                       <  1:5024977‑M1/50‑7 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga                                    <  1:1266971‑M1/50‑10 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgg                                     <  1:3590025‑M1/49‑9 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgg                                     <  1:5112150‑M1/49‑9 (MQ=255)
            aTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacgg                                 >  1:6380618‑M1/1‑38 (MQ=255)
             tttCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAcgagatggacggt                                >  1:3826483‑M1/1‑37 (MQ=255)
             tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggt                                <  1:1979390‑M1/50‑14 (MQ=255)
               tCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtc                               <  1:1093888‑M1/49‑15 (MQ=255)
               tCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtc                               <  1:2513981‑M1/49‑15 (MQ=255)
                ccAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcc                              <  1:4308974‑M1/49‑16 (MQ=255)
                ccAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcc                              <  1:3771797‑M1/49‑16 (MQ=255)
                 cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc                            <  1:4325121‑M1/50‑18 (MQ=255)
                 cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc                            <  1:381948‑M1/50‑18 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                           >  1:2497159‑M1/1‑32 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                           <  1:2492481‑M1/50‑19 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                           >  1:4730446‑M1/1‑32 (MQ=255)
                   ggTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctc                          <  1:1336313‑M1/50‑20 (MQ=255)
                     tCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgc                        >  1:5615000‑M1/1‑29 (MQ=255)
                      cAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcc                       >  1:3300477‑M1/1‑28 (MQ=255)
                       aGACGGGCTGCAAGTTGCAGACCGTTAcgcgctggacggtcccctcgccc                      >  1:6318880‑M1/1‑27 (MQ=255)
                       aGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccc                      >  1:4184341‑M1/1‑27 (MQ=255)
                       aGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccc                      >  1:2344429‑M1/1‑27 (MQ=255)
                        gACGGGCTGCAAGTTGCAGACCGTTAcgagatgtacggtcccctcgcccc                     >  1:5674970‑M1/1‑26 (MQ=255)
                        gACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccc                     >  1:2650843‑M1/1‑26 (MQ=255)
                         aCGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccc                    >  1:3653509‑M1/1‑25 (MQ=255)
                         aCGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccc                    >  1:2776085‑M1/1‑25 (MQ=255)
                         aCGGGCTGCAAGTTGCAGACAGTTAcgagatggacggtcccctcgccccc                    <  1:164577‑M1/50‑26 (MQ=255)
                          cGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccct                   >  1:5062140‑M1/1‑24 (MQ=255)
                            ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttg                 >  1:1319998‑M1/1‑22 (MQ=255)
                             gCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgg                >  1:3514537‑M1/1‑21 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg               >  1:3445342‑M1/1‑17 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg               >  1:258996‑M1/1‑17 (MQ=255)
                                  aaGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggga             >  1:3711977‑M1/1‑16 (MQ=255)
                                   aGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagag          >  1:1957106‑M1/1‑15 (MQ=255)
                                    gTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagg         >  1:6473029‑M1/1‑14 (MQ=255)
                                     ttGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagaggg        <  1:1193157‑M1/50‑38 (MQ=255)
                                      tGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt       >  1:2736283‑M1/1‑12 (MQ=255)
                                        cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtt      <  1:6581203‑M1/49‑40 (MQ=255)
                                         aGACCGTTAcgagatggacggtcccctcgcccccttggggagagggttag    <  1:265911‑M1/50‑42 (MQ=255)
                                          gACCGTTAcgagatggacggtcccctcgcccccttggggagagggttagg   <  1:6096446‑M1/50‑43 (MQ=255)
                                           aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttaggg  <  1:4497640‑M1/50‑44 (MQ=255)
                                           aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttaggg  >  1:4328052‑M1/1‑7 (MQ=255)
                                           aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttaggg  <  1:3650776‑M1/50‑44 (MQ=255)
                                                 |                                           
GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/380062‑380154

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