New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =4 (0.070)55 (1.350) 27/102 0.1 87.3% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 13 (0.320)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209552‑1209683
                                                                   |                                                                
ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                  >  1:1057000/1‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                  >  1:3355872‑M2/68‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                  >  1:1037512‑M2/68‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                  >  1:1046344‑M2/68‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                <  1:919591‑M2/3‑1 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                >  1:3980267‑M2/66‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                >  1:3839197‑M2/66‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                >  1:2591264‑M2/66‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                >  1:2149948‑M2/66‑68 (MQ=255)
       aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                           >  1:3731590‑M2/61‑68 (MQ=255)
        accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                          >  1:3127921‑M2/60‑68 (MQ=255)
        accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                          >  1:1996742‑M2/60‑68 (MQ=255)
           tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaac                                                       >  1:3826772/1‑67 (MQ=255)
             tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                     >  1:3802502‑M2/55‑68 (MQ=255)
             tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                     >  1:2557016‑M2/55‑68 (MQ=255)
              gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                    >  1:144280‑M2/54‑68 (MQ=255)
                gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:1530075/1‑68 (MQ=255)
                gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:2883571/1‑68 (MQ=255)
                gCCACGAAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:2503181/1‑68 (MQ=255)
                     aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             <  1:1915019‑M2/22‑1 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:1995675‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:1801060‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:320743‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:3359702‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:1658983‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:3683489‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:1066096‑M2/45‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:95339‑M2/44‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:2271468‑M2/44‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:2513766‑M2/44‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:135801‑M2/44‑68 (MQ=255)
                         tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:2971706‑M2/43‑68 (MQ=255)
                          agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        >  1:3655048‑M2/42‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCTCTTATGAAAATGAAGAAATAACAAGCa                                      >  1:2628535‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      >  1:267459‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttcgcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:127001‑M2/29‑1 (MQ=255)
                             cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                     >  1:2839621‑M2/39‑68 (MQ=255)
                              gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    <  1:2745267‑M2/31‑1 (MQ=255)
                                 ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                 >  1:340342/1‑68 (MQ=255)
                                  ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:1466313‑M2/35‑1 (MQ=255)
                                   tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               <  1:1133027‑M2/36‑1 (MQ=255)
                                   tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               <  1:565405‑M2/36‑1 (MQ=255)
                                      attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            >  1:3914044‑M2/30‑68 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:693096‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:4000669‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1843225‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3681593‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1717633‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:3372880‑M2/28‑68 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3400735‑M2/41‑1 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcatc                          >  1:2712951/1‑67 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:473182/1‑68 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1723591/68‑1 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1915025/68‑1 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:3996998‑M2/25‑68 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:1062916‑M2/25‑68 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCtgt                     >  1:1617805‑M2/23‑67 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:1687272‑M2/23‑68 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2044012‑M2/46‑1 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:185343‑M2/46‑1 (MQ=255)
                                             gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:926406/1‑68 (MQ=255)
                                             gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:123645/1‑68 (MQ=255)
                                              cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:2038204/1‑68 (MQ=255)
                                                gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:3708749‑M2/20‑68 (MQ=255)
                                                 ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                 <  1:1155907‑M2/50‑1 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:340896‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:1484754‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:710335‑M2/18‑68 (MQ=255)
                                                    cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc              <  1:1253096‑M2/53‑1 (MQ=255)
                                                    cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc              >  1:3741702‑M2/16‑68 (MQ=255)
                                                    cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc              >  1:3489868‑M2/16‑68 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:282641‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:3522062‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAAGGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:605279‑M2/59‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:958502‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATCGCATCATTCCTGCTTTTACCAGGGGGa       <  1:333687‑M2/60‑1 (MQ=255)
                                                                caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:826600‑M2/4‑68 (MQ=255)
                                                                   |                                                                
GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209552‑1209683

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