New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078057 (0.120)52 (1.160) 15/96 0.5 84.4% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 14 (0.310)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666
                                                                |                                                                
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                  <  1:4692079/65‑1 (MQ=255)
 cagccggcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:3782751‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1531946‑M2/1‑64 (MQ=255)
 cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                 >  1:1606339‑M2/1‑64 (MQ=255)
    ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                              <  1:471163‑M2/65‑5 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                           >  1:1620922‑M2/1‑58 (MQ=255)
        tctGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          <  1:2348251‑M2/64‑9 (MQ=255)
        ccTGTTGGCGGAGTGGCTGAAGGGCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          >  1:22149‑M2/1‑57 (MQ=255)
        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                          >  1:1881070‑M2/1‑57 (MQ=255)
            ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                      >  1:1703625‑M2/1‑53 (MQ=255)
                 ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                 <  1:3715500‑M2/65‑18 (MQ=255)
                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                <  1:1498706‑M2/65‑19 (MQ=255)
                         tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                         >  1:2732647/1‑65 (MQ=255)
                          gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                        >  1:1430512‑M2/1‑39 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     >  1:516962‑M2/1‑36 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:4185271‑M2/65‑30 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     >  1:598995‑M2/1‑36 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:2020365‑M2/65‑30 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                     <  1:683373‑M2/65‑30 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAAGCAACGTTTATGAAAAt                                     >  1:4184806/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:205717‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2873976‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2159262‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                    <  1:2184361‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggatttcgagagcc                                    <  1:879415‑M2/65‑31 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:4584871/1‑65 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:2934148/1‑65 (MQ=255)
                               gcaCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2477006‑M2/2‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCGGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:4448888‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCCGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:164436‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaatgttttcgagagccg                                   >  1:3679537‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:4480276‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:506506‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:569036‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:677258‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3912751‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:75000‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3729451‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3392315‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3616311‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3615734‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3489750‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:3342740‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2103881‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1661672‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:2244083‑M2/1‑34 (MQ=255)
                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                   >  1:1646628‑M2/1‑34 (MQ=255)
                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                              >  1:1148632‑M2/1‑29 (MQ=255)
                                      aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcca                            <  1:372481‑M2/65‑39 (MQ=255)
                                         gggACACCGACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:4615526‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:195448‑M2/1‑24 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         <  1:537391‑M2/65‑42 (MQ=255)
                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                         >  1:2667046‑M2/1‑24 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                       <  1:1170645‑M2/65‑44 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacagtgc                       <  1:15400‑M2/65‑44 (MQ=255)
                                                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat             >  1:643721/1‑65 (MQ=255)
                                                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat             >  1:1109080/1‑65 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:3925085‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:4333984‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg            <  1:3203359‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttcgcacagtgccatccgatttg            <  1:2198425‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggatttcgagagccgtttgcacagtgccatccgatttg            <  1:4515596‑M2/65‑55 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc            >  1:3140235/1‑65 (MQ=255)
                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt         <  1:1422043/65‑1 (MQ=255)
                                                          tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat        <  1:3799520‑M2/65‑59 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:2093321‑M2/65‑60 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc       <  1:259904‑M2/65‑60 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:641342‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg      <  1:1136586‑M2/65‑61 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggtcttcgagagccgtttgtacagtgccatccgatttgatatcg      <  1:3973701‑M2/65‑61 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:1503892/1‑65 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:3235719/1‑65 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt      >  1:839308/1‑65 (MQ=255)
                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc    <  1:4761709/65‑1 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt  >  1:3516660/1‑65 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCtt  >  1:4192390/1‑65 (MQ=255)
                                                                |                                                                
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT  >  NC_000913/1209538‑1209666

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