New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780515 (0.380)26 (0.870) 10/96 0.7 56.4% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 29 (0.970)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666
                                                               |                                                                
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCCCGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1854551‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTATTGCAGAGc                                                                 >  1:2425641‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:483986‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1257592‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGAGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1784459‑M2/1‑64 (MQ=255)
     agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                            >  1:1228782‑M2/1‑59 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:1612230‑M2/1‑58 (MQ=255)
           ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                      >  1:1047067‑M2/1‑53 (MQ=255)
           ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                      >  1:1346357‑M2/1‑53 (MQ=255)
              gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                   <  1:1948916‑M2/65‑16 (MQ=255)
                  acaggctgAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:297922‑M2/4‑46 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                               >  1:558249/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagtgcc                                    <  1:116076‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:805851‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2279338‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:1077805/1‑65 (MQ=255)
                              ccACGGGACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3186706‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1644113‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3265621‑M2/1‑34 (MQ=255)
                                  ggAACCGGGCCACCAACAGGTAATGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATgaagaa                               >  1:347053/1‑65 (MQ=255)
                                  ggAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttg                               <  1:1907652‑M2/65‑36 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1266378‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                              >  1:243900/1‑65 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                              >  1:2342342/1‑65 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                              >  1:2380111/1‑65 (MQ=255)
                                        gggACCCCAACCGGTAATGCAAAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3206595‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAGTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2012939‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3171676‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1308768‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1863204‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAAGGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3269393‑M2/1‑24 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:1437485‑M2/65‑44 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:2530081‑M2/65‑44 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                       >  1:3288686/1‑65 (MQ=255)
                                                aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                 >  1:857570/1‑65 (MQ=255)
                                                 aCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                <  1:2479570‑M2/65‑51 (MQ=255)
                                                  cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc               <  1:3123791/65‑1 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:789946‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:891717‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagcagtttgcacagtgccatccgatttg            <  1:2812204‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:2511208/1‑65 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:935826/1‑65 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:549187/1‑65 (MQ=255)
                                                      tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca           <  1:3057358/65‑1 (MQ=255)
                                                       aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat          <  1:822234/65‑1 (MQ=255)
                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:859501/65‑1 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1047752‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:2936105‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1087228‑M2/65‑60 (MQ=255)
                                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc       >  1:384679/1‑65 (MQ=255)
                                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc       >  1:2927282/1‑65 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:916528‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttaatatcc      <  1:3290920‑M2/65‑61 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:258239/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:2529887/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:654273/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:1683509/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:1584584/1‑65 (MQ=255)
                                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg     >  1:754556/1‑65 (MQ=255)
                                                             gcgCCTTCTCCCAAACCAACGTTTATGAAAGTGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:1345387/63‑1 (MQ=255)
                                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:385796/65‑1 (MQ=255)
                                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:680443/65‑1 (MQ=255)
                                                              agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt   <  1:2139681/65‑1 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:2418380/1‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:1773481/1‑65 (MQ=255)
                                                               |                                                                
CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666

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