New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =0 (0.000)53 (1.280) 23/96 0.1 99.1% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 1 (0.020)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209559‑1209683
                                                            |                                                                
aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccctctcccaaaccaaCGttt                                                              >  1:3968586‑M2/61‑65 (MQ=255)
      tgcgtcataatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                        >  1:2013678‑M2/55‑65 (MQ=255)
       gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                       >  1:1735755‑M2/54‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:978648‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:667471‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:2855837‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:44914‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:4076005‑M2/45‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:1152895‑M2/44‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:3733162‑M2/44‑65 (MQ=255)
                  tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                            >  1:2953066‑M2/43‑65 (MQ=255)
                   agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:4201210‑M2/42‑65 (MQ=255)
                    gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          <  1:1293111‑M2/25‑1 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3199832‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:264376‑M2/40‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:2178503‑M2/40‑65 (MQ=255)
                      cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        <  1:3925530‑M2/27‑1 (MQ=255)
                      cgttgctgcattgagttgtgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        <  1:3730016‑M2/27‑1 (MQ=255)
                       gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                       >  1:2071421‑M2/38‑65 (MQ=255)
                          gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:2648538‑M2/35‑65 (MQ=255)
                           ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                   >  1:2127286‑M2/34‑65 (MQ=255)
                            tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                  >  1:129778‑M2/33‑65 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             <  1:1630288‑M2/38‑1 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:4139125‑M2/28‑65 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:1408275‑M2/28‑65 (MQ=255)
                                  gagttttgctccttctcccaaaccaaCGTCTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:3338708‑M2/39‑1 (MQ=255)
                                   agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:2422020‑M2/40‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:993201‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:633178‑M2/25‑65 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3528345‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:592181‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:3646112‑M2/25‑65 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3680491‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:579728‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:3818178‑M2/25‑65 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:118121‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:122222‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1331020‑M2/41‑1 (MQ=255)
                                      tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:3068570‑M2/23‑65 (MQ=255)
                                      tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:132235‑M2/23‑65 (MQ=255)
                                      tttgctccttctaccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:2666219‑M2/23‑65 (MQ=255)
                                       cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:90927/1‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:296794‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:47586‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:1574605‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:1586178‑M2/20‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:4158206‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:268188‑M2/20‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2265254‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2303047‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:1906091‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:2453753‑M2/20‑65 (MQ=255)
                                         gctccttcgcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:1658041‑M2/46‑1 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcctgtct                   >  1:4306838‑M2/18‑62 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                <  1:3473727‑M2/51‑1 (MQ=255)
                                                   ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           >  1:1453371‑M2/10‑65 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTNCTGCTTTTACCAGGGGGa       <  1:3246377‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2340221‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:755676‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:350344‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:549669‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:4138933‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1737937‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:4044048‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1786677‑M2/60‑1 (MQ=255)
                                                         caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:2089519‑M2/4‑65 (MQ=255)
                                                            cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:1361924/1‑65 (MQ=255)
                                                            |                                                                
TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209559‑1209683

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