New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078058 (0.180)39 (1.170) 18/96 0.2 77.2% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 17 (0.510)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666
                                                               |                                                                
cagccggcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1212139‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:405717‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:3075927‑M2/1‑64 (MQ=255)
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1546685‑M2/1‑64 (MQ=255)
                 gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacccagg                                                <  1:1572529‑M2/65‑19 (MQ=255)
                 gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                <  1:1570364‑M2/65‑19 (MQ=255)
                 gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                <  1:3459737/65‑1 (MQ=255)
                  aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:1384369‑M2/1‑46 (MQ=255)
                    tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt                                             <  1:547342‑M2/65‑22 (MQ=255)
                       cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcg                                          <  1:414889‑M2/65‑25 (MQ=255)
                            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:2551208‑M2/65‑30 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgcgagcc                                    <  1:669565‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagtc                                    <  1:1851227‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:1746904‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:1359272‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:1046760/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:2228342/1‑65 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:1566971/1‑65 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1761165‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1993368‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:613997‑M2/1‑34 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:2315556‑M2/1‑29 (MQ=255)
                                        gggACGCCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2394565‑M2/1‑24 (MQ=255)
                                        gggACACCAACATGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1820062‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:856865‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:351703‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3440993‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:308656‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1131680‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                         <  1:2530935/65‑1 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:2823930‑M2/65‑44 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:2912869‑M2/65‑44 (MQ=255)
                                             acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                    >  1:769384/1‑65 (MQ=255)
                                              ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                   <  1:1080884‑M2/65‑48 (MQ=255)
                                              ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                   <  1:3012429/65‑1 (MQ=255)
                                                  cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc               <  1:1514973/65‑1 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:2877891‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:1386968‑M2/65‑55 (MQ=255)
                                                      tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca           <  1:2319318/65‑1 (MQ=255)
                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:1079167/65‑1 (MQ=255)
                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:1490592/65‑1 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1798950‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:3209133‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggatttcgagagccgtatgacaagtgccatccgatctgatatc       <  1:2497386‑M2/65‑60 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:2964914‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:2641301‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      >  1:251826‑M2/1‑5 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:868677/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:949114/1‑65 (MQ=255)
                                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:3364987/1‑65 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     >  1:2226315‑M2/1‑4 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:2181416‑M2/65‑62 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:1072673‑M2/65‑62 (MQ=255)
                                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:427901/65‑1 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:170442‑M2/65‑65 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:1689925‑M2/65‑65 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:1377238‑M2/65‑65 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:869473‑M2/65‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:602001/1‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:1673560/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 < 26 ≤ 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.