New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078056 (0.160)28 (0.980) 12/96 0.4 80.6% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 9 (0.320)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGC  >  NC_000913/1209538‑1209664
                                                                |                                                              
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACCCCAACAGGTAATGCagag                                                                >  1:2298623/1‑65 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                               >  1:1568280‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                               >  1:2764314‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                               >  1:2135479‑M2/1‑64 (MQ=255)
  agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                              >  1:589787‑M2/1‑63 (MQ=255)
      agccTGTTGGCGGAGCGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                          >  1:67478‑M2/1‑59 (MQ=255)
            ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                    >  1:869898‑M2/1‑53 (MQ=255)
                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGCCTTCTCCCAAACCAACg                                              >  1:2162906/1‑65 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                   <  1:1720748‑M2/65‑30 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                  <  1:1032705‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                  <  1:2860451‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                  <  1:1834810‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                  <  1:263858‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagcc                                  <  1:2457430‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGGCTTCTCCCAAACCAACGCTTATGAAAATg                                  >  1:2677307/1‑65 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                 >  1:1569720‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                 >  1:2326056‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                 >  1:2452913‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                 >  1:89624‑M2/1‑34 (MQ=255)
                                 gcgGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                               >  1:2105779‑M2/2‑32 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacgaaggttttcgagagccgtttgc                            <  1:2163897‑M2/65‑37 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                            <  1:2434223‑M2/65‑37 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                            >  1:530097‑M2/1‑29 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                            >  1:1169592‑M2/1‑29 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                       >  1:2950133‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                       >  1:2567694‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                       >  1:2146306‑M2/1‑24 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                     <  1:2351825‑M2/65‑44 (MQ=255)
                                            acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcc                    <  1:272186‑M2/65‑45 (MQ=255)
                                            acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcc                    <  1:1345267‑M2/65‑45 (MQ=255)
                                                cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc                <  1:1188179‑M2/65‑49 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc          >  1:174533/1‑65 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc          >  1:2838931/1‑65 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt       <  1:2037948/65‑1 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt       <  1:1351460/65‑1 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt       >  1:514073/1‑65 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc     <  1:554067‑M2/65‑60 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc     <  1:1742993‑M2/65‑60 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg    <  1:82452‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg    <  1:2545022‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtcccatccgatttgatatcg    <  1:404457‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaagggtttcgagagccgtttgcacagtgccatccgatttgatatcg    <  1:1690273‑M2/65‑61 (MQ=255)
                                                              tagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc  <  1:403694/64‑1 (MQ=255)
                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc  <  1:1551818/65‑1 (MQ=255)
                                                                |                                                              
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGC  >  NC_000913/1209538‑1209664

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