New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078055 (0.110)32 (0.960) 16/96 0.3 81.3% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 11 (0.330)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666
                                                                |                                                                
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGAagag                                                                  >  1:3192560/1‑65 (MQ=255)
 ccgccagcctgTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:648465‑M2/3‑64 (MQ=255)
 cagcgagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:143703‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:286272‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:2758801‑M2/1‑64 (MQ=255)
               gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                   <  1:856187‑M2/65‑16 (MQ=255)
               gCGGAGTAGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                   <  1:1104505‑M2/65‑16 (MQ=255)
                   aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:418055‑M2/1‑46 (MQ=255)
                   aGTGGCTGAAGGCCACAGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:210063‑M2/1‑46 (MQ=255)
                         tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                         >  1:1362310/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2414741‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3519030‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3579020‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggatttcgagagcc                                    <  1:1745368‑M2/65‑31 (MQ=255)
                               cgccggaCCGGGGACACCAACAGGTAATGCAGAGcctgctcacaaaccaaggttttcgagagccg                                   >  1:1242820‑M2/4‑34 (MQ=255)
                               ccACGGAGCCGTGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3272654‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:832715‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2535028‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1915285‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:660898‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2049592‑M2/1‑34 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              <  1:2939561‑M2/65‑37 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2031275‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                         <  1:2468724/65‑1 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                       <  1:1621202/65‑1 (MQ=255)
                                            acacCAACAGGTAATGCAGAGccttctcccaaaccacggttttcgagagccgtttgcacagtgcc                      <  1:2780928‑M2/65‑45 (MQ=255)
                                            acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcc                      >  1:2099221‑M2/1‑21 (MQ=255)
                                             cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcca                     <  1:3611886‑M2/65‑46 (MQ=255)
                                                   cAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgat               >  1:1940445‑M2/1‑14 (MQ=255)
                                                    aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca              >  1:1416628/1‑65 (MQ=255)
                                                     ggTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt             <  1:135654‑M2/65‑54 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:950218‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:466930‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:2983000‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:2517012‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:2383946/1‑65 (MQ=255)
                                                          tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat        <  1:554662‑M2/65‑59 (MQ=255)
                                                           gCAGAGcctcctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:2442749‑M2/65‑60 (MQ=255)
                                                           gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc       >  1:1708182/1‑65 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:3069018‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:3101626‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:287723‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:1041839‑M2/65‑61 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:1653760/1‑65 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:2833638/1‑65 (MQ=255)
                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:2841976‑M2/65‑62 (MQ=255)
                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:1106683/65‑1 (MQ=255)
                                                              cagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:3473996/64‑1 (MQ=255)
                                                                gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgagttgatatcgccaa  <  1:131875‑M2/65‑65 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:2945726/1‑65 (MQ=255)
                                                                |                                                                
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666

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