New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.160)37 (1.300) 24/96 0.0 84.6% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 9 (0.320)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209555‑1209680
                                                                |                                                             
atcaaaccggtgcgtcctatggtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                               <  1:2826314‑M2/1‑1 (MQ=255)
 gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGCCTTCTCCCAAACCAACg                                                              >  1:2162906/1‑65 (MQ=255)
     accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:473275‑M2/60‑65 (MQ=255)
     accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:168034‑M2/60‑65 (MQ=255)
     accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:1545590‑M2/60‑65 (MQ=255)
         gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                      <  1:2468579‑M2/10‑1 (MQ=255)
          tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                     >  1:1587618‑M2/55‑65 (MQ=255)
           gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                    >  1:2238191‑M2/54‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGGCTTCTCCCAAACCAACGCTTATGAAAATg                                                  >  1:2677307/1‑65 (MQ=255)
              tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                 <  1:199217‑M2/15‑1 (MQ=255)
                ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                               <  1:1678180‑M2/17‑1 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:663451‑M2/45‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:2155956‑M2/45‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:1376383‑M2/45‑65 (MQ=255)
                      tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                         >  1:891708‑M2/43‑65 (MQ=255)
                       agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:419340‑M2/42‑65 (MQ=255)
                        gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                       >  1:2183809‑M2/41‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:410442‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:2411050‑M2/40‑65 (MQ=255)
                              gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:33718‑M2/31‑1 (MQ=255)
                              gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:895542‑M2/31‑1 (MQ=255)
                               ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                >  1:2678206‑M2/34‑65 (MQ=255)
                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                               >  1:1821160‑M2/33‑65 (MQ=255)
                                 gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                              <  1:2419743‑M2/34‑1 (MQ=255)
                                  cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  1:1417506‑M2/35‑1 (MQ=255)
                                     tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          <  1:1125098‑M2/38‑1 (MQ=255)
                                     tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:2194553‑M2/28‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:174533/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:2838931/1‑65 (MQ=255)
                                      gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         <  1:2192532‑M2/39‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:697508‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCACCAtt                       <  1:1519575‑M2/41‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:2037948/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:1351460/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:514073/1‑65 (MQ=255)
                                             tagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:403694/64‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:2645865‑M2/20‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:546857‑M2/20‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:2472217‑M2/20‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:697563‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2110494‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:1488569‑M2/46‑1 (MQ=255)
                                             gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:1551818/65‑1 (MQ=255)
                                              ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 <  1:317371‑M2/47‑1 (MQ=255)
                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:1244387‑M2/18‑65 (MQ=255)
                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:112215‑M2/18‑65 (MQ=255)
                                                  ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             <  1:2221729‑M2/51‑1 (MQ=255)
                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg        <  1:30639‑M2/56‑1 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:365603‑M2/59‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:104324‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1440037‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    >  1:934495‑M2/6‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  >  1:1498139‑M2/4‑65 (MQ=255)
                                                                |                                                             
GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209555‑1209680

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