New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =15 (0.380)25 (0.840) 17/96 0.2 55.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 29 (0.970)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGA  >  NC_000913/1209557‑1209678
                                                              |                                                           
aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                           >  1:558249/1‑65 (MQ=255)
   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                        >  1:3230902‑M2/60‑65 (MQ=255)
   accggtgcgtcctaatatagccgttgctgcattgagttttgctccttctcccaaaccaaCgtttt                                                        >  1:2002696‑M2/60‑64 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                   >  1:2097070‑M2/55‑65 (MQ=255)
           gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                >  1:1077805/1‑65 (MQ=255)
                ggAACCGGGCCACCAACAGGTAATGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATgaagaa                                           >  1:347053/1‑65 (MQ=255)
                 gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                          >  1:243900/1‑65 (MQ=255)
                 gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                          >  1:2380111/1‑65 (MQ=255)
                 gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                          >  1:2342342/1‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                         >  1:1377319‑M2/45‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                        >  1:2367687‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                        >  1:1191621‑M2/44‑65 (MQ=255)
                      gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                     <  1:805639‑M2/25‑1 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAGAATGAAGAAATAACaa                                    >  1:3171871‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                    >  1:1365462‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                    >  1:236504‑M2/40‑65 (MQ=255)
                        gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                   >  1:3288686/1‑65 (MQ=255)
                        cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                   >  1:1192392‑M2/39‑65 (MQ=255)
                            gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                               <  1:2053806‑M2/31‑1 (MQ=255)
                             ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                              >  1:342763‑M2/34‑65 (MQ=255)
                              aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                             >  1:857570/1‑65 (MQ=255)
                                cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                           <  1:3123791/65‑1 (MQ=255)
                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                          >  1:278468‑M2/30‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                        >  1:549187/1‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                        >  1:2511208/1‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                        >  1:935826/1‑65 (MQ=255)
                                    tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                       <  1:3057358/65‑1 (MQ=255)
                                     aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                      <  1:822234/65‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                     <  1:206429‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                     >  1:1908641‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                     >  1:744316‑M2/25‑65 (MQ=255)
                                      aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                     <  1:859501/65‑1 (MQ=255)
                                        tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                   >  1:1028671‑M2/23‑65 (MQ=255)
                                        gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                   >  1:2927282/1‑65 (MQ=255)
                                        gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                   >  1:384679/1‑65 (MQ=255)
                                         cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                  >  1:654273/1‑65 (MQ=255)
                                         cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                  >  1:258239/1‑65 (MQ=255)
                                         cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                  >  1:2529887/1‑65 (MQ=255)
                                         cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                  >  1:1683509/1‑65 (MQ=255)
                                         cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                  >  1:1584584/1‑65 (MQ=255)
                                          agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                 >  1:754556/1‑65 (MQ=255)
                                           tctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:1701186‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:2708353‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                >  1:737452‑M2/20‑65 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:709043‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:445434‑M2/46‑1 (MQ=255)
                                           gcgCCTTCTCCCAAACCAACGTTTATGAAAGTGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:1345387/63‑1 (MQ=255)
                                           gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:385796/65‑1 (MQ=255)
                                           gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                <  1:680443/65‑1 (MQ=255)
                                            ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt               <  1:2481624‑M2/47‑1 (MQ=255)
                                            agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt               <  1:2139681/65‑1 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt              >  1:36750‑M2/18‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt              >  1:2418380/1‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt              >  1:1773481/1‑65 (MQ=255)
                                                ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa           <  1:424401‑M2/51‑1 (MQ=255)
                                                ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa           <  1:2109020‑M2/51‑1 (MQ=255)
                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg       >  1:1996062‑M2/11‑65 (MQ=255)
                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg      >  1:111741‑M2/10‑65 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg   <  1:990062‑M2/59‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg   <  1:1808922‑M2/59‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg   <  1:2014775‑M2/59‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg   <  1:380161‑M2/59‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa  <  1:624221‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa  <  1:1322701‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa  <  1:2576652‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa  <  1:3080814‑M2/60‑1 (MQ=255)
                                                              |                                                           
AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGA  >  NC_000913/1209557‑1209678

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.