New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.120)47 (1.270) 22/96 0.1 84.3% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 13 (0.350)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209557‑1209681
                                                              |                                                              
gaacccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                              <  1:1566452‑M2/3‑1 (MQ=255)
aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                              <  1:3590711/65‑1 (MQ=255)
  aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                            >  1:604417‑M2/61‑65 (MQ=255)
   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                           >  1:1433109‑M2/60‑65 (MQ=255)
      tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                        >  1:2219987/1‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:3404790‑M2/55‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:3773168‑M2/55‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:145435‑M2/55‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:370269‑M2/55‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:421‑M2/55‑65 (MQ=255)
          ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                    >  1:207621/1‑65 (MQ=255)
           gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                   >  1:1855262‑M2/52‑65 (MQ=255)
           gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                   >  1:208635/1‑65 (MQ=255)
               taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                               <  1:692598‑M2/18‑1 (MQ=255)
                aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                              >  1:1944596‑M2/47‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:2661662‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:329436‑M2/45‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:1659628‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:513043‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:2521819‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:261462‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:2972491‑M2/44‑65 (MQ=255)
                     agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  1:2101534‑M2/42‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:684724‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:1573495‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:2330866‑M2/40‑65 (MQ=255)
                          ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                    <  1:930625‑M2/29‑1 (MQ=255)
                          ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                    >  1:3582286‑M2/37‑65 (MQ=255)
                           acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                   >  1:302887/1‑65 (MQ=255)
                               gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                               <  1:2270630‑M2/34‑1 (MQ=255)
                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                             <  1:3291256‑M2/36‑1 (MQ=255)
                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                             >  1:2959584‑M2/30‑65 (MQ=255)
                                   tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:1319201‑M2/28‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:25741/1‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:90704/1‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGgcatcatg                        <  1:3433034‑M2/41‑2 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        >  1:643745‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        >  1:818007‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:3549420‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:1978597‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        >  1:1262511‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:3846953‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:1256368‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:378395‑M2/41‑1 (MQ=255)
                                       tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                       >  1:3229312/1‑65 (MQ=255)
                                        tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:2608543‑M2/23‑65 (MQ=255)
                                        gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:752907/1‑65 (MQ=255)
                                        gCAGAGCCTTCTCCCAAACCAACGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:1653229/1‑65 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:2385334‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:1953321‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:434750‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:2426296‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:1450871‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:3143318‑M2/46‑1 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:3757566‑M2/18‑65 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:2092492‑M2/18‑65 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:1810168‑M2/18‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAGTGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:2395992/1‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:2149448/1‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:1642897/1‑65 (MQ=255)
                                               cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt               >  1:1019925‑M2/16‑65 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  1:711678‑M2/59‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGGAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:3127919‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:469151‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:1483559‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:2839221‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:90902‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:1222456‑M2/60‑1 (MQ=255)
                                                          ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt    >  1:161299‑M2/5‑65 (MQ=255)
                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt   >  1:2435012‑M2/4‑65 (MQ=255)
                                                            aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:2299487‑M2/63‑1 (MQ=255)
                                                              |                                                              
AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209557‑1209681

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