New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078056 (0.140)35 (1.120) 15/96 0.3 80.9% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 12 (0.380)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666
                                                                |                                                                
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                  <  1:2819973/65‑1 (MQ=255)
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGAGACACCCACAGGTAATGCagag                                                                  >  1:2716078/1‑65 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1205486‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:2116858‑M2/1‑64 (MQ=255)
   gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                               <  1:2405312‑M2/65‑4 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:1630780‑M2/1‑58 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:2322109‑M2/1‑58 (MQ=255)
            ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                      >  1:2643754‑M2/1‑53 (MQ=255)
                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCCGAGCCTTCTCCCAAACCAACg                                                >  1:2009316/1‑65 (MQ=255)
                              gCCACGGAAGCGGGACACCAACAGGCAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:1308973‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagaggc                                    <  1:2675056‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:218399‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:172297/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:209132/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCCTCTCCCAAACAAACGTTTATGAAAATg                                    >  1:2427509/1‑65 (MQ=255)
                               ccACGGAGCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1384814‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1763534‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:299740‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3151289‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3353428‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2715681‑M2/1‑34 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:3060905‑M2/1‑29 (MQ=255)
                                         gggACCCCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2012703‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTATTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2948664‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:789121‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:668579‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3153836‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:593114‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2402967‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1655879‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2163958‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:17537‑M2/1‑24 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       >  1:2148688‑M2/1‑22 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:480826‑M2/65‑44 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                       >  1:566019/1‑65 (MQ=255)
                                             cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcca                     >  1:643148‑M2/1‑20 (MQ=255)
                                              acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                    <  1:2003842‑M2/65‑47 (MQ=255)
                                               ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                   >  1:1245601/1‑65 (MQ=255)
                                                 aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                 <  1:3124948/65‑1 (MQ=255)
                                                  aCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                >  1:1732273‑M2/1‑15 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:831233‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:2227512‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattcg            <  1:2039101‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:227719/1‑65 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:2543502/65‑1 (MQ=255)
                                                          tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat        <  1:1528576‑M2/65‑59 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatctgatttgatatc       <  1:1044077‑M2/65‑60 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3365264‑M2/65‑60 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:2062880‑M2/65‑61 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:2811523/1‑65 (MQ=255)
                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:3133141‑M2/65‑62 (MQ=255)
                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:3236245‑M2/65‑62 (MQ=255)
                                                                gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:3208306‑M2/65‑65 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCCAATGGCATCATTCCTGCtt  >  1:1689408/1‑65 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:1835006/1‑65 (MQ=255)
                                                                |                                                                
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666

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