New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =4 (0.090)43 (1.290) 22/96 0.1 86.0% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 11 (0.330)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209559‑1209681
                                                            |                                                              
aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                            >  1:3510371‑M2/61‑65 (MQ=255)
    tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                        >  1:1362310/1‑65 (MQ=255)
      tgcgtcctaatgtggccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:2592431‑M2/55‑65 (MQ=255)
      tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:1814956‑M2/55‑65 (MQ=255)
      tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:2357068‑M2/55‑65 (MQ=255)
          tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                  <  1:1929063‑M2/15‑1 (MQ=255)
           cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                 >  1:651698‑M2/50‑65 (MQ=255)
           cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                 >  1:19047‑M2/50‑65 (MQ=255)
               atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                             <  1:988343‑M2/20‑1 (MQ=255)
               atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                             >  1:600754‑M2/46‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:3393328‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:2960775‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:2558111‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:1696992‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:1246171‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:1784178‑M2/45‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           >  1:1352212‑M2/44‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                           <  1:1344175‑M2/22‑1 (MQ=255)
                   agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  1:1925865‑M2/42‑65 (MQ=255)
                    gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                        <  1:2468724/65‑1 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:1879983‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:2087370‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       <  1:560560‑M2/26‑1 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:2207503‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:3629836‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                       >  1:1307728‑M2/40‑65 (MQ=255)
                      gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                      <  1:1621202/65‑1 (MQ=255)
                         tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                   <  1:2918030‑M2/30‑1 (MQ=255)
                            tgcattgagttttgctccttctcccaaaccaaCGTTTATCAGAATGAATAATTGACAAGCAAATg                                >  1:3490846‑M2/33‑65 (MQ=255)
                              cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                              <  1:1888127‑M2/35‑1 (MQ=255)
                               attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                             >  1:869625‑M2/30‑65 (MQ=255)
                               aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                             >  1:1416628/1‑65 (MQ=255)
                                 gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:2383946/1‑65 (MQ=255)
                                  gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                          <  1:2620391‑M2/39‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2357787‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2471580‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2499416‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2571085‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2885394‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        >  1:1727451‑M2/25‑65 (MQ=255)
                                      gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:1708182/1‑65 (MQ=255)
                                       cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  1:1653760/1‑65 (MQ=255)
                                       cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  1:2833638/1‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:627198‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:1123817‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   >  1:2614985‑M2/20‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:2174026‑M2/46‑1 (MQ=255)
                                         gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:1106683/65‑1 (MQ=255)
                                         cagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:3473996/64‑1 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:1647917‑M2/18‑65 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:1768674‑M2/18‑65 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:3275543‑M2/18‑65 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:3110099‑M2/18‑65 (MQ=255)
                                           gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:2945726/1‑65 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                >  1:2422972‑M2/17‑65 (MQ=255)
                                                ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc            <  1:1473413‑M2/53‑1 (MQ=255)
                                                  cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg          >  1:438368‑M2/11‑65 (MQ=255)
                                                   ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg         <  1:2581695‑M2/56‑1 (MQ=255)
                                                   ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg         >  1:2195318‑M2/10‑65 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg        <  1:2661750‑M2/57‑1 (MQ=255)
                                                      aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  1:915685‑M2/59‑1 (MQ=255)
                                                      aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  1:11253‑M2/59‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:3147893‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:2496268‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:45339‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:54126‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:286531‑M2/60‑1 (MQ=255)
                                                        ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt    >  1:1810807‑M2/5‑65 (MQ=255)
                                                        ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt    >  1:2213597‑M2/5‑65 (MQ=255)
                                                          aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:1307056‑M2/63‑1 (MQ=255)
                                                            |                                                              
TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209559‑1209681

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