New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078052 (0.040)49 (1.290) 18/96 0.3 88.7% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 11 (0.290)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666
                                                               |                                                                
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:106533‑M2/1‑64 (MQ=255)
cacccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 <  1:1439780/62‑1 (MQ=255)
   ccagccTGTTGGCGGAGTGGCTGCAGGCCACGGGACCGGCACACCAACAGGTAATGCAGAGcctt                                                              <  1:3834126‑M2/65‑5 (MQ=255)
   ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                              <  1:1625058‑M2/65‑5 (MQ=255)
     tgcctGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                            >  1:978806‑M2/2‑59 (MQ=255)
     agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                            >  1:1893670‑M2/1‑59 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:2554034‑M2/1‑58 (MQ=255)
       ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          >  1:3010476‑M2/1‑57 (MQ=255)
               cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaagccaa                                                  <  1:469499‑M2/65‑17 (MQ=255)
                          aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                       >  1:3589668/1‑65 (MQ=255)
                            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:192060‑M2/65‑30 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:3074273‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:4033809‑M2/65‑31 (MQ=255)
                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:3889196/1‑65 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:693691‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1205333‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2905210‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1023395‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3250648‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1274955‑M2/1‑34 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1324700‑M2/1‑34 (MQ=255)
                                   gAAGCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1993208‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:345078‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:3784707‑M2/1‑29 (MQ=255)
                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcgagagccgtttgc                              <  1:2101534‑M2/65‑37 (MQ=255)
                                     aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                            <  1:2977278‑M2/65‑39 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:4040852‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:469123‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:387817‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3162100‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1358153‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:3657099‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:634585‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1552576‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:10938‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1544718‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:888723‑M2/1‑24 (MQ=255)
                                        gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1808951‑M2/1‑24 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:381341‑M2/65‑44 (MQ=255)
                                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                       >  1:364228/1‑65 (MQ=255)
                                              ccAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacagtgcgatc                   <  1:3881467‑M2/65‑48 (MQ=255)
                                                aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                 >  1:3350054/1‑65 (MQ=255)
                                                   aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt              <  1:48841‑M2/65‑53 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgctttg            <  1:2485802‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:570834‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:791522‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:3396915‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgcttgcacagtgccatccgatttg            <  1:1363273‑M2/65‑55 (MQ=255)
                                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:1885884/1‑65 (MQ=255)
                                                      tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga           >  1:197357‑M2/1‑10 (MQ=255)
                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:980932/65‑1 (MQ=255)
                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:3252761/65‑1 (MQ=255)
                                                         tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat        <  1:2783387‑M2/65‑59 (MQ=255)
                                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc        >  1:1278946/1‑65 (MQ=255)
                                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc        >  1:2362191/1‑65 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:4002441‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1804158‑M2/65‑60 (MQ=255)
                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:1933358‑M2/65‑60 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      >  1:3742830‑M2/1‑5 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:919355‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      >  1:1552142‑M2/1‑5 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:811623‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:2612806‑M2/65‑61 (MQ=255)
                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttaatatcg      <  1:1982214‑M2/65‑61 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:3435171‑M2/65‑62 (MQ=255)
                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:1643173‑M2/65‑62 (MQ=255)
                                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:2972821/65‑1 (MQ=255)
                                                               gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:3466735‑M2/65‑65 (MQ=255)
                                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:479985/1‑65 (MQ=255)
                                                               |                                                                
CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209539‑1209666

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