New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078050 (0.000)58 (1.400) 23/96 0.1 99.2% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 1 (0.020)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666
                                                               |                                                                
cggccaggctgTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1485622‑M2/3‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:2829444‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:3939969‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:103636‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:47092‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:2955713‑M2/1‑64 (MQ=255)
 agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGCAATGCagagac                                                                >  1:4169821/1‑63 (MQ=255)
  gccagccTGTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                               <  1:659523‑M2/65‑4 (MQ=255)
   ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                              <  1:3669205‑M2/65‑5 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:1710393‑M2/1‑58 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:365229‑M2/1‑58 (MQ=255)
       ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          <  1:2178950‑M2/65‑9 (MQ=255)
             ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                    <  1:1594335‑M2/65‑15 (MQ=255)
              gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTTATGCAGAGccttctcccaaacca                                                   <  1:3423275‑M2/65‑16 (MQ=255)
                ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                 <  1:1972256‑M2/65‑18 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:226274‑M2/1‑46 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:3961304‑M2/1‑46 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:882579‑M2/1‑46 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:3278828‑M2/1‑46 (MQ=255)
                    tGGCTGAAGGCCACGGAACCGGGACACCAACAGGCAATGCAGAGccttctcccaaaccaagggtt                                             <  1:1653416‑M2/65‑22 (MQ=255)
                        tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcga                                         <  1:2525073‑M2/65‑26 (MQ=255)
                            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:474590‑M2/65‑30 (MQ=255)
                             gCCACGNAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2588452‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGCCACCAACAGGTAATGCAGAGccttctcccaaaacaaggttttcgagagcc                                    <  1:3503448‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:245720‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3501839‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:892207‑M2/65‑31 (MQ=255)
                              tcaCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3082969‑M2/2‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGtcttctcccaaaccaaggttttcgagagccg                                   >  1:2641024‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2568056‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1049904‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:245331‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2182982‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:14526‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:4251963‑M2/1‑34 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:3265660‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              <  1:1885874‑M2/65‑37 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1255036‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgcttgc                              <  1:1315833‑M2/65‑37 (MQ=255)
                                     aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                            <  1:2869355‑M2/65‑39 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:4265645‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1376847‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1879259‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1880513‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcgagagacgtttgcacagt                         >  1:2675805‑M2/1‑24 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:4127271‑M2/65‑44 (MQ=255)
                                             acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                    <  1:2151819‑M2/65‑47 (MQ=255)
                                                aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg                 <  1:2722182‑M2/65‑50 (MQ=255)
                                                aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagcgccatccg                 <  1:3149512‑M2/65‑50 (MQ=255)
                                                 cagggtaaTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                >  1:1206782‑M2/4‑15 (MQ=255)
                                                  cAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgat               >  1:1163777‑M2/1‑14 (MQ=255)
                                                   aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt              <  1:2507085‑M2/65‑53 (MQ=255)
                                                     gTAATGCAGAGccttctcccaacccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:4310542‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:2031740‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:1979994‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:1682316‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:1483608‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:4198085‑M2/65‑55 (MQ=255)
                                                      tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga           >  1:2933922‑M2/1‑10 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:2838067‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3811001‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagcagtttgcacagtgccatccgatttgatatc       <  1:4168306‑M2/65‑60 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:1782067‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:4271677‑M2/65‑61 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:90927/1‑65 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:4137338‑M2/65‑62 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:1993470‑M2/65‑62 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:1739511‑M2/65‑62 (MQ=255)
                                                              agcctcntcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca   <  1:2590920‑M2/65‑64 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgcgatccgatttgatatcgcgaa  <  1:1942018‑M2/65‑65 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccac  <  1:3096409‑M2/65‑65 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  >  1:1878657‑M2/1‑1 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:1390788‑M2/65‑65 (MQ=255)
                                                               |                                                                
CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666

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