New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 566000 =NA (NA)68 (1.370) 32/98 NT 100% intergenic (+90/‑56) ybcD/insE DLP12 prophage region; ECK0532:JW0527:b4508; predicted replication protein fragment/IS3 element protein InsE
?W3110S.gb = 576913 0 (0.000)coding (78/498 nt) ybcS predicted lysozyme

ATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGG  >  W3110S.gb/565932‑566056
                                                                    |                                                        
aTTACCTGGGGACTCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTg                                                         <  1:713371/70‑1 (MQ=255)
     cTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATATGATCTTa                                                   >  1:1642528/1‑71 (MQ=255)
        gggACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTAc                                                  <  1:1377016/69‑1 (MQ=255)
             gCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCaa                                           <  1:2693016/71‑1 (MQ=255)
              cATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAAt                                          <  1:1123361/71‑1 (MQ=255)
               aTAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGa                                                        <  1:4713752/56‑1 (MQ=255)
               aTAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGa                                                        <  1:3118948/56‑1 (MQ=255)
                tAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTg                                                         <  1:4411118/54‑1 (MQ=255)
                tAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAg                                        >  1:887788/1‑71 (MQ=255)
                tAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAg                                        <  1:4456650/71‑1 (MQ=255)
                  agcgCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:1551638/4‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:3376605/1‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:3880666/1‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:3909331/1‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:4213396/1‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:4578348/1‑71 (MQ=255)
                  agagCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      >  1:4754797/1‑71 (MQ=255)
                     gCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTAccc                                                >  1:127315/1‑58 (MQ=255)
                     gCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAAt                                          <  1:309106/64‑1 (MQ=255)
                     gCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGa                                    >  1:3141362/1‑70 (MQ=255)
                        tCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAg                                        >  1:4415761/1‑63 (MQ=255)
                        tCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACAcg                                >  1:2025316/1‑71 (MQ=255)
                         cTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAAt                                          <  1:2365887/60‑1 (MQ=255)
                         cTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACAcgc                               >  1:3142237/1‑71 (MQ=255)
                          tgtgCACAGGACGTTGCAGCGGTAACAGGCGCAACCGTAATCTGAt                                                       <  1:1853937/46‑1 (MQ=255)
                            tgCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTAcc                                                 <  1:2443916/50‑1 (MQ=255)
                            tgCACAGGACGTTGCCGCGCTAACAGGCGCAACCGTAATCTATCTTACCCAGCAATAGTGGACACGCGGCt                           <  1:4704026/71‑1 (MQ=37)
                              cacaGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCt                           >  1:4364095/1‑70 (MQ=255)
                              cacaGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCt                           <  1:178364/70‑1 (MQ=255)
                               acaGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTaa                         <  1:3208710/71‑1 (MQ=255)
                                 aGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTg                                                         <  1:993691/37‑1 (MQ=255)
                                 aGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTg                                                         <  1:220959/37‑1 (MQ=255)
                                 aGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                       <  1:1197304/71‑1 (MQ=255)
                                   gACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                       >  1:753632/1‑69 (MQ=255)
                                   gACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGa                     <  1:3528274/71‑1 (MQ=255)
                                     cGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCaa                                           >  1:1720054/1‑47 (MQ=255)
                                     cGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATa                                         <  1:3432370/49‑1 (MQ=255)
                                     cGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATa                                         <  1:2685655/49‑1 (MQ=255)
                                     cGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATa                                         <  1:4382992/49‑1 (MQ=255)
                                      gTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCaa                                           <  1:3726404/46‑1 (MQ=255)
                                       ttGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTAcc                                                 <  1:1420160/39‑1 (MQ=255)
                                       ttGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTg                                      <  1:3258368/50‑1 (MQ=255)
                                        tGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATa                                         >  1:194383/1‑46 (MQ=255)
                                        tGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaaa                >  1:2819817/1‑71 (MQ=255)
                                         gCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCa                                               >  1:328318/1‑39 (MQ=255)
                                         gCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCa                                               >  1:479147/1‑39 (MQ=255)
                                           cgcgGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGT‑‑ACACGCGGCTAAGTGAGTa                  <  1:74525/64‑1 (MQ=21)
                                              ggTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAActct            >  1:3962876/1‑69 (MQ=255)
                                                tAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAg                                        >  1:2916535/1‑39 (MQ=255)
                                                tAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt        >  1:713170/1‑71 (MQ=255)
                                                  aCAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTa                          >  1:2564735/1‑51 (MQ=255)
                                                   cAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtcag     <  1:105905/71‑1 (MQ=255)
                                                   cAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt        >  1:1678193/1‑68 (MQ=255)
                                                   cAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt        >  1:3610648/1‑68 (MQ=255)
                                                      gcgcAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAgg  >  1:4471287/1‑71 (MQ=255)
                                                         cAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtcag     <  1:3220695/65‑1 (MQ=39)
                                                           aCCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtca      >  1:44388/1‑62 (MQ=35)
                                                            ccGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                       <  1:2467777/44‑1 (MQ=35)
                                                                    |                                                        
ATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGG  >  W3110S.gb/565932‑566056

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