New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780518 (0.330)27 (0.670) 14/96 0.5 57.1% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 27 (0.670)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666
                                                               |                                                                
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACTGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1197075‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1337238‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1692472‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGCAATGCAGAGc                                                                 >  1:483429‑M2/1‑64 (MQ=255)
   ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                              <  1:3255200‑M2/65‑5 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:861022‑M2/1‑58 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:3060207‑M2/1‑58 (MQ=255)
                  aGTGGCTGAAGGCCGCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:3769665‑M2/1‑46 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:453609‑M2/1‑46 (MQ=255)
                    tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGttt                                             >  1:2787429/1‑65 (MQ=255)
                          aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                       >  1:336133/1‑65 (MQ=255)
                            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:2556741‑M2/65‑30 (MQ=255)
                            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                     >  1:3718900/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3034893‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:1198848‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagcc                                    <  1:132234‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTg                                    >  1:1212541/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:323732/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:316534/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:487369/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:3507150/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:3534202/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3036614‑M2/65‑31 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2388066‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:782205‑M2/1‑34 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:3682529‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                              >  1:623747/1‑65 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2255235‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2177791‑M2/1‑24 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:2861344‑M2/65‑44 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                       >  1:446843/1‑65 (MQ=255)
                                           acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                      >  1:3668653/1‑65 (MQ=255)
                                            cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcca                     >  1:3974376‑M2/1‑20 (MQ=255)
                                            cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCa                     <  1:1392881/65‑1 (MQ=255)
                                              ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                   <  1:2491624/65‑1 (MQ=255)
                                                 aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                >  1:2151598/1‑65 (MQ=255)
                                                  cAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgat               >  1:638969‑M2/1‑14 (MQ=255)
                                                   aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca              >  1:2749580/1‑65 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:4041077‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:896074‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:3074896‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:2930274/1‑65 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:2836903/1‑65 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:1860399/1‑65 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:2314801/1‑65 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3698746‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1642258‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3322075‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3972182‑M2/65‑60 (MQ=255)
                                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc       >  1:852266/1‑65 (MQ=255)
                                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc       >  1:781514/1‑65 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:2784376‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:794669‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:4293814‑M2/65‑61 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:3638213‑M2/65‑62 (MQ=255)
                                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg     >  1:217212/1‑65 (MQ=255)
                                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:369957/65‑1 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:2179200‑M2/65‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  <  1:2583804/65‑1 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:1334259/1‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCtt  >  1:2923112/1‑65 (MQ=255)
                                                               |                                                                
CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666

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