New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.160)48 (1.260) 23/102 0.1 82.7% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 14 (0.370)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NC_000913/1209552‑1209686
                                                                   |                                                                   
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                     >  1:3948964‑M2/68‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                     >  1:1997057‑M2/68‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                     >  1:1459516‑M2/68‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                   >  1:1724162‑M2/66‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                   >  1:2341657‑M2/66‑68 (MQ=255)
       aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                              >  1:2383213‑M2/61‑68 (MQ=255)
        accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                             >  1:1653729‑M2/60‑68 (MQ=255)
           tGAAGGCCACGGAACCGGGACACCAAAAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                          >  1:4020387/1‑68 (MQ=255)
             tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                        >  1:3031641‑M2/55‑68 (MQ=255)
             tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                        >  1:6456‑M2/55‑68 (MQ=255)
                gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                     >  1:853225/1‑68 (MQ=255)
                gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                     >  1:3878709/1‑68 (MQ=255)
                   ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                  >  1:2613483‑M2/49‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                              >  1:232343‑M2/45‑68 (MQ=255)
                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                              >  1:2596843‑M2/45‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                             >  1:3083555‑M2/44‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                             >  1:2404372‑M2/44‑68 (MQ=255)
                        aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                             >  1:59235/1‑68 (MQ=255)
                          agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                           >  1:1448938‑M2/42‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGACGAAATAACAAGCa                                         >  1:711820‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGCAATAACAAGCa                                         >  1:3445653‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  1:4005592‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  1:3431461‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  1:134439‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  1:3230183‑M2/40‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  1:1904977‑M2/40‑68 (MQ=255)
                             gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                        >  1:1105303/1‑68 (MQ=255)
                             cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                        >  1:1417357‑M2/39‑68 (MQ=255)
                              gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                       <  1:2864260‑M2/31‑1 (MQ=255)
                              acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                       <  1:1120779/68‑1 (MQ=255)
                               ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                      >  1:828810‑M2/37‑68 (MQ=255)
                               cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                      >  1:2552484/1‑68 (MQ=255)
                                acCAACAGGTAATGCAGAGCCTTCTGCCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAACTg                                     >  1:463261/1‑68 (MQ=255)
                                 gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                    <  1:739369‑M2/34‑1 (MQ=255)
                                   tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGACGAAATAACAAGCAAATGGca                                  >  1:2201962‑M2/33‑68 (MQ=255)
                                   tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                  >  1:2240705‑M2/33‑68 (MQ=255)
                                   aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                  >  1:2864027/1‑68 (MQ=255)
                                      attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                               >  1:831744‑M2/30‑68 (MQ=255)
                                      attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                               <  1:1213507‑M2/39‑1 (MQ=255)
                                       ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                              <  1:3323647‑M2/40‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             <  1:3644121‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             <  1:486725‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             <  1:1043500‑M2/41‑1 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             >  1:1046500‑M2/28‑68 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             >  1:2355073‑M2/28‑68 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             >  1:224564‑M2/28‑68 (MQ=255)
                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             <  1:2101065‑M2/41‑1 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             >  1:1670929/1‑68 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                          >  1:370180‑M2/25‑68 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                          >  1:1499215‑M2/25‑68 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                          >  1:1946452‑M2/25‑68 (MQ=255)
                                            tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTg                         >  1:2707404/1‑68 (MQ=255)
                                            tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                         >  1:1952180/1‑68 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                        <  1:3306408‑M2/46‑1 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                        <  1:303609‑M2/46‑1 (MQ=255)
                                               agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                      >  1:7095/1‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   >  1:3974041‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   >  1:801913‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   >  1:1443675‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   >  1:3807621‑M2/18‑68 (MQ=255)
                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   >  1:2047892‑M2/18‑68 (MQ=255)
                                                  gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                   <  1:3715096/68‑1 (MQ=255)
                                                         cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg            <  1:3062473‑M2/58‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:3469702‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:171581‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:2225694‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:2350933‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:4020732‑M2/59‑1 (MQ=255)
                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg           <  1:2419894‑M2/59‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:1026214‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:903096‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:587318‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:3802334‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:3271658‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:2989239‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:2238900‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa          <  1:1788414‑M2/60‑1 (MQ=255)
                                                            aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt         <  1:2641646‑M2/61‑1 (MQ=255)
                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt       >  1:1471917‑M2/6‑68 (MQ=255)
                                                                caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa     >  1:2677175‑M2/4‑68 (MQ=255)
                                                                   cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt  >  1:45792/1‑68 (MQ=255)
                                                                   |                                                                   
GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NC_000913/1209552‑1209686

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