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

AAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCC  >  NC_000913/380063‑380150
                                                |                                       
aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                        <  1:2053335‑M1/50‑2 (MQ=255)
 aGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcg                                       <  1:336424‑M1/50‑3 (MQ=255)
    aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgaga                                    <  1:446502‑M1/50‑6 (MQ=255)
       cAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg                                  >  1:632363‑M1/1‑42 (MQ=255)
        aGAATCTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgg                                 <  1:1208150‑M1/49‑9 (MQ=255)
         gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggac                               <  1:1261839‑M1/50‑11 (MQ=255)
         gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggac                               >  1:2233346‑M1/1‑40 (MQ=255)
           aTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacgg                             >  1:913277‑M1/1‑38 (MQ=255)
           aTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacgg                             <  1:2227528‑M1/50‑13 (MQ=255)
              tCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtc                           <  1:1261561‑M1/49‑15 (MQ=255)
                cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc                        >  1:1664135‑M1/1‑33 (MQ=255)
                cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc                        >  1:171193‑M1/1‑33 (MQ=255)
                     cAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcc                   >  1:940143‑M1/1‑28 (MQ=255)
                         cGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccct               <  1:348486‑M1/50‑27 (MQ=255)
                          gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt              <  1:1352855‑M1/50‑28 (MQ=255)
                          gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt              <  1:81795‑M1/50‑28 (MQ=255)
                           ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttg             >  1:1398580‑M1/1‑22 (MQ=255)
                            gCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgg            <  1:1025323‑M1/50‑30 (MQ=255)
                                cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggg          <  1:145944‑M1/48‑32 (MQ=255)
                                 aaGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggga         >  1:1528390‑M1/1‑16 (MQ=255)
                                  aGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagag      <  1:1636427‑M1/50‑36 (MQ=255)
                                   gTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagg     <  1:462883‑M1/50‑37 (MQ=255)
                                      gCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt   >  1:1115953‑M1/1‑11 (MQ=255)
                                       cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtt  <  1:636809‑M1/49‑40 (MQ=255)
                                       cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt   >  1:2130177‑M1/1‑10 (MQ=255)
                                                |                                       
AAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCC  >  NC_000913/380063‑380150

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