New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =31 (0.780)10 (0.340) 10/96 0.7 24.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 37 (1.240)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209555‑1209681
                                                                |                                                              
ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                <  1:320636/65‑1 (MQ=255)
        tGAAGGCCACGGAACCGGGGCACCAACAGGTAATGCCGAGCCTTCTTCCAAACCAACGTTTATGa                                                        >  1:1014723/1‑65 (MQ=255)
        tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                        >  1:2778477/1‑65 (MQ=255)
        tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                        >  1:471835/1‑65 (MQ=255)
         ctgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                       <  1:1744529‑M2/10‑1 (MQ=255)
          tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                      >  1:692380‑M2/55‑65 (MQ=255)
          aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                      >  1:529154/1‑65 (MQ=255)
            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                    >  1:3185566/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTGTCCCAAACCAACGTTTATGGAAATg                                                   >  1:2190891/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                   >  1:2480178/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                   >  1:259423/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                   >  1:1666387/1‑65 (MQ=255)
               cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                 >  1:2879812‑M2/50‑65 (MQ=255)
                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTGCCAAACCAACGATTATGAAAATGAAGaaa                                             >  1:3080392/1‑65 (MQ=255)
                   gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                             >  1:194601/1‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            >  1:2399470‑M2/45‑65 (MQ=255)
                        gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                        >  1:174142‑M2/41‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCCTTTATGAAAATGAAGCAATAACaa                                       >  1:1300248‑M2/40‑65 (MQ=255)
                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                      >  1:1679737/1‑65 (MQ=255)
                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                      >  1:1129911/1‑65 (MQ=255)
                          cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                      >  1:1550915‑M2/39‑65 (MQ=255)
                              ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                  <  1:3062400/65‑1 (MQ=255)
                                aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                >  1:313205/1‑65 (MQ=255)
                                  cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                              <  1:532102‑M2/35‑1 (MQ=255)
                                  cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                              <  1:1347968/65‑1 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:2108688/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:266032/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:2860803/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:592973/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                           >  1:1189358/1‑65 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:1034928‑M2/41‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:731384/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:2474365/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:839665/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  1:1259511/65‑1 (MQ=255)
                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                       >  1:261966/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:493429/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                      >  1:2688988/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGGAGAGATAACAAAGAAATGGCATCATTCCt                     >  1:3088569/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  1:246684/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  1:2277967/1‑65 (MQ=255)
                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                    >  1:1713316/1‑65 (MQ=255)
                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                    >  1:1520035/1‑65 (MQ=255)
                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:1061894/65‑1 (MQ=255)
                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 <  1:1506381‑M2/48‑1 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:395945/1‑65 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                 >  1:2199766/1‑65 (MQ=255)
                                                ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                >  1:623066‑M2/17‑65 (MQ=255)
                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt               >  1:1176734‑M2/16‑65 (MQ=255)
                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt               >  1:75235‑M2/16‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg          >  1:3096702‑M2/11‑65 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  1:1115390‑M2/59‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:3226316‑M2/60‑1 (MQ=255)
                                                              agcGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:1567197/63‑1 (MQ=255)
                                                              aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:1023854‑M2/63‑1 (MQ=255)
                                                                |                                                              
GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209555‑1209681

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.