New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)38 (0.600)
+CGAGAT
25/86 0.6 58.7% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 4070444 = 61 (0.840)intergenic (‑43/+71) yihR/yihS putative aldose 1‑epimerase YihR/sulfoquinovose isomerase

GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC  >  NC_000913/380062‑380156
                                                 |                                             
gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                               <  1:511334/50‑1 (MQ=255)
gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                               <  1:2924284/50‑1 (MQ=255)
  aGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                              >  1:914907‑M1/1‑48 (MQ=255)
  aGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                              >  1:2306419‑M1/1‑48 (MQ=255)
    gAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                              >  1:3710784‑M1/1‑46 (MQ=255)
    gAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                              >  1:3084726‑M1/1‑46 (MQ=255)
     aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgaga                                          <  1:2858639‑M1/50‑6 (MQ=255)
     aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                              >  1:1131895‑M1/1‑45 (MQ=255)
       tCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg                                        <  1:742434‑M1/50‑8 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga                                      <  1:2278166‑M1/50‑10 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga                                      <  1:2498214‑M1/50‑10 (MQ=255)
         aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgg                                       <  1:5237534‑M1/49‑9 (MQ=255)
          gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggac                                     >  1:1724603‑M1/1‑40 (MQ=255)
          gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggac                                     <  1:4571882‑M1/50‑11 (MQ=255)
           aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg                                    >  1:5037066‑M1/1‑39 (MQ=255)
             tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggt                                  >  1:1509447‑M1/1‑37 (MQ=255)
              ttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtc                                 >  1:5264778‑M1/1‑36 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                             >  1:2070707‑M1/1‑32 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                             >  1:1953615‑M1/1‑32 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                             >  1:2973251‑M1/1‑32 (MQ=255)
                  aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccct                             >  1:4837674‑M1/1‑32 (MQ=255)
                     tCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggccggtcccctcgc                          >  1:2976270‑M1/1‑29 (MQ=255)
                     tCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgc                          >  1:344132‑M1/1‑29 (MQ=255)
                           gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt                    <  1:5322702‑M1/50‑28 (MQ=255)
                           gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt                    <  1:25771‑M1/50‑28 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggg                <  1:4307872‑M1/48‑32 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg                 >  1:3993109‑M1/1‑17 (MQ=255)
                                 cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg                 >  1:5505688‑M1/1‑17 (MQ=255)
                                  aaGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggga               >  1:4662005‑M1/1‑16 (MQ=255)
                                   aGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagag            >  1:3046267‑M1/1‑15 (MQ=255)
                                   aGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagag            <  1:5475860‑M1/50‑36 (MQ=255)
                                    gTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagg           >  1:3172422‑M1/1‑14 (MQ=255)
                                     ttGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagaggg          <  1:1061763‑M1/50‑38 (MQ=255)
                                     ttGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagaggg          <  1:4603521‑M1/50‑38 (MQ=255)
                                     ttGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagaggg          <  1:4931723‑M1/50‑38 (MQ=255)
                                      tGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt         >  1:829751‑M1/1‑12 (MQ=255)
                                       gCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt         >  1:3905649‑M1/1‑11 (MQ=255)
                                        cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtta       <  1:1767362‑M1/50‑41 (MQ=255)
                                        cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtt        <  1:178277‑M1/49‑40 (MQ=255)
                                        cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt         >  1:927468‑M1/1‑10 (MQ=255)
                                         aGACCGTTAcgagatggacggtcccctcgcccccttggggagagggttag      <  1:4161290‑M1/50‑42 (MQ=255)
                                         aGACCGTTAcgagatggacggtcccctcgcccccttggggagagggttag      >  1:3693795‑M1/1‑9 (MQ=255)
                                          gACCGTTAcgagatggacggtcccctcgcccccttggggagagggttagg     >  1:1319448‑M1/1‑8 (MQ=255)
                                            ccGTTAcgagatggacggtcccctcgcccccttggggagagggttagggt   <  1:462770‑M1/50‑45 (MQ=255)
                                             cGTTAcgagatggacggtcccctcgcccccttggggagagggttagggtg  <  1:3974927‑M1/50‑46 (MQ=255)
                                             cGTTAcgagatggacggtcccctcgcccccttggggagagggttagggtg  <  1:1729925‑M1/50‑46 (MQ=255)
                                                 |                                             
GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC  >  NC_000913/380062‑380156

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