New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780515 (0.260)28 (0.640) 15/96 0.5 68.1% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 15 (0.340)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666
                                                                |                                                                
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCagag                                                                  >  1:3737148/1‑65 (MQ=255)
 gcgccatcctgTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:471275‑M2/3‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:51508‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1490226‑M2/1‑64 (MQ=255)
 cacccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 <  1:2196168/62‑1 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:305019‑M2/1‑58 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:2591677‑M2/1‑58 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:1796771‑M2/1‑58 (MQ=255)
        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          >  1:831104‑M2/1‑57 (MQ=255)
        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          <  1:1860071‑M2/65‑9 (MQ=255)
            ttGGAGGAGTGGCTGAAGGCCACGGAACCGGGAAGCCAACAGGTAATGCAGAGccttctcccaaa                                                      >  1:752354‑M2/1‑53 (MQ=255)
              ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                    >  1:1117037‑M2/1‑51 (MQ=255)
                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctcccaaacccagg                                                <  1:3003716‑M2/65‑19 (MQ=255)
                   aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                               >  1:1335409‑M2/1‑46 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                     >  1:3938653/1‑65 (MQ=255)
                             ggCCACGGAACCGGGACAACAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:4485095‑M2/65‑30 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:851236‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2804109‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2364804‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2433801‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:4455770/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:1243479/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:3496122/1‑65 (MQ=255)
                               ccACGGACCCGGGACACCAACAGGTAATGCAGAGcctgctcccaaaccaaggttttcgagagccg                                   >  1:721014‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2869352‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:611825‑M2/1‑34 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1606313‑M2/1‑29 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1463363‑M2/1‑29 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              <  1:1230985‑M2/65‑37 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:1458809‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:335475‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2040593‑M2/1‑24 (MQ=255)
                                                   cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc               <  1:2528736/65‑1 (MQ=255)
                                                    aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt              <  1:1785227‑M2/65‑53 (MQ=255)
                                                    aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca              >  1:2029359/1‑65 (MQ=255)
                                                    aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca              >  1:1996604/1‑65 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacaggggcatccgatttg            <  1:2774009‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:349730/1‑65 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:4429324/1‑65 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:4379717/1‑65 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:4517093/65‑1 (MQ=255)
                                                          tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat        <  1:1439225‑M2/65‑59 (MQ=255)
                                                          tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc        >  1:1891486/1‑65 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:1145926‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:3944579‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:1529848‑M2/65‑61 (MQ=255)
                                                            cAGAGCCTTCACCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCt      >  1:4209862/1‑65 (MQ=255)
                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:4509849‑M2/65‑62 (MQ=255)
                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc     <  1:1684438‑M2/65‑62 (MQ=255)
                                                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg     >  1:864755/1‑65 (MQ=255)
                                                               agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca   <  1:4451817‑M2/65‑64 (MQ=255)
                                                                gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  >  1:4379891‑M2/1‑1 (MQ=255)
                                                                gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa  <  1:1816337‑M2/65‑65 (MQ=255)
                                                                gccttctcccaaaccaaggttttcgagagccgtttgcacagagccatccgatttgatatcgccaa  <  1:779825‑M2/65‑65 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:1169357/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 < 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.