New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.160)44 (1.170) 25/96 0.1 88.0% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 6 (0.160)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683
                                                                |                                                                
atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                  <  1:2548878‑M2/1‑1 (MQ=255)
   aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                               >  1:1357956‑M2/62‑65 (MQ=255)
    aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                              >  1:3474651‑M2/61‑65 (MQ=255)
    aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                              >  1:3519632‑M2/61‑65 (MQ=255)
          tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                        >  1:1863258‑M2/55‑65 (MQ=255)
               cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                   >  1:1228914‑M2/50‑65 (MQ=255)
                   atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                               >  1:2430862‑M2/46‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:2978093‑M2/45‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:2398273‑M2/45‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:622543‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:3014908‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:3195327‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:2520011‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             <  1:802517‑M2/22‑1 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:1427497‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:1737227‑M2/44‑65 (MQ=255)
                        gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:1814269‑M2/41‑65 (MQ=255)
                        gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          <  1:3706066‑M2/25‑1 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3660341‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:2101758‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:2108392‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:234547‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3767933‑M2/40‑65 (MQ=255)
                          cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        <  1:3506263‑M2/27‑1 (MQ=255)
                           gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                       <  1:1819916‑M2/28‑1 (MQ=255)
                              gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:2503226‑M2/35‑65 (MQ=255)
                              gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:387757‑M2/35‑65 (MQ=255)
                                  cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:3278052‑M2/35‑1 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:3368995‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:3162390‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:3981888‑M2/30‑65 (MQ=255)
                                     tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:3946743‑M2/28‑65 (MQ=255)
                                     tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             <  1:1861526‑M2/38‑1 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:1180460/1‑65 (MQ=255)
                                      gagttttgctccttctcccaaaccaaCGTTTATAAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:1467269‑M2/39‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:263873‑M2/25‑65 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3196983‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3484677‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2234905‑M2/41‑1 (MQ=255)
                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGCAAATGAACAAATAACAAGCAACTGGCATCATTc                         >  1:814240/1‑65 (MQ=255)
                                          tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:413971‑M2/23‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:160374/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:910171/1‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2897224‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:3404296‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:260359‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:709963‑M2/46‑1 (MQ=255)
                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:3360816‑M2/18‑65 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1989933/1‑65 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1813966/1‑65 (MQ=255)
                                                ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:3327141‑M2/17‑65 (MQ=255)
                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                 >  1:2992367‑M2/16‑65 (MQ=255)
                                                  ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                <  1:2636229‑M2/51‑1 (MQ=255)
                                                  ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                <  1:455439‑M2/51‑1 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:2050970‑M2/59‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2401868‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:3694472‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1157387‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1475611‑M2/60‑1 (MQ=255)
                                                            ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt      >  1:2752109‑M2/5‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:2946874‑M2/4‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:1990521‑M2/4‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:1762605‑M2/4‑65 (MQ=255)
                                                              aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt    <  1:918968‑M2/63‑1 (MQ=255)
                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:2170220/1‑65 (MQ=255)
                                                                |                                                                
GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683

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