New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =7 (0.120)38 (0.850) 21/96 0.2 79.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 14 (0.310)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209560‑1209680
                                                           |                                                             
accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:468692‑M2/60‑65 (MQ=255)
   tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                       >  1:2732647/1‑65 (MQ=255)
      gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                    >  1:2426618‑M2/54‑65 (MQ=255)
       ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAAGCAACGTTTATGAAAAt                                                   >  1:4184806/1‑65 (MQ=255)
        gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:4584871/1‑65 (MQ=255)
        gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:2934148/1‑65 (MQ=255)
          cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                >  1:4015946‑M2/50‑65 (MQ=255)
            taatttagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  1:1421183‑M2/18‑1 (MQ=255)
            taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  1:3817347‑M2/18‑1 (MQ=255)
               tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:1176891‑M2/45‑65 (MQ=255)
               tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:2282899‑M2/45‑65 (MQ=255)
                  agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAGATGAAGAAATAAc                                        >  1:1328108‑M2/42‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1498930‑M2/40‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:3404123‑M2/40‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1342368‑M2/40‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1406778‑M2/40‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1882351‑M2/40‑65 (MQ=255)
                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:159416‑M2/40‑65 (MQ=255)
                     cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                     <  1:758701‑M2/27‑1 (MQ=255)
                     cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                     >  1:1239515‑M2/39‑65 (MQ=255)
                        tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                  <  1:3502880‑M2/30‑1 (MQ=255)
                         gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:542288‑M2/31‑1 (MQ=255)
                              attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                            >  1:1812582‑M2/30‑65 (MQ=255)
                              attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                            <  1:1193397‑M2/36‑1 (MQ=255)
                               ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                           >  1:643721/1‑65 (MQ=255)
                               ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                           >  1:1109080/1‑65 (MQ=255)
                                gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:3140235/1‑65 (MQ=255)
                                 gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         >  1:1682889‑M2/27‑65 (MQ=255)
                                   gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:2712635‑M2/41‑1 (MQ=255)
                                   gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:794883‑M2/25‑65 (MQ=255)
                                   gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:4516008‑M2/25‑65 (MQ=255)
                                   gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:2612621‑M2/25‑65 (MQ=255)
                                   gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:3804536‑M2/25‑65 (MQ=255)
                                   aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:1422043/65‑1 (MQ=255)
                                     tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:21060‑M2/23‑65 (MQ=255)
                                     tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:4541768‑M2/23‑65 (MQ=255)
                                      cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:1503892/1‑65 (MQ=255)
                                      cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:3235719/1‑65 (MQ=255)
                                      cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:839308/1‑65 (MQ=255)
                                        gctccttctcccaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:3665880‑M2/20‑65 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:4597249‑M2/46‑1 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:4619485‑M2/46‑1 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:4730221‑M2/46‑1 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:865567‑M2/20‑65 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:3511110‑M2/46‑1 (MQ=255)
                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:3149623‑M2/46‑1 (MQ=255)
                                        gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:4761709/65‑1 (MQ=255)
                                         ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 <  1:4153130‑M2/47‑1 (MQ=255)
                                          tccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:823346‑M2/18‑65 (MQ=255)
                                          tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:4404449‑M2/18‑65 (MQ=255)
                                          gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:3516660/1‑65 (MQ=255)
                                          gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCtt                >  1:4192390/1‑65 (MQ=255)
                                           tcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt               >  1:3120384‑M2/17‑65 (MQ=255)
                                            cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt              <  1:2485882‑M2/50‑1 (MQ=255)
                                                     aaccaaCGTTTATGAAAATGAAGGAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:2735843‑M2/59‑1 (MQ=255)
                                                     aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:1367317‑M2/59‑1 (MQ=255)
                                                     aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:1582362‑M2/59‑1 (MQ=255)
                                                     aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:789506‑M2/59‑1 (MQ=255)
                                                     aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:4720088‑M2/59‑1 (MQ=255)
                                                      accgaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:169375‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:4577509‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:41050‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1981333‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1766883‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:740150‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:4030900‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:4027642‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:3517342‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:138479‑M2/60‑1 (MQ=255)
                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:3824883‑M2/60‑1 (MQ=255)
                                                      aCCAGCGTTTATCAAAATGAAGGAATAACCAGCAAATGGCATCAGTCCTGCTTTGACCAGGGGGa    <  1:2260433/65‑1 (MQ=255)
                                                        caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  <  1:3447985‑M2/62‑1 (MQ=255)
                                                           |                                                             
GGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209560‑1209680

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