New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =15 (0.260)38 (0.870) 22/96 0.2 74.3% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 15 (0.340)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683
                                                                |                                                                
atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                  <  1:63086‑M2/1‑1 (MQ=255)
 tgaacccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                                 <  1:3149421‑M2/2‑1 (MQ=255)
  gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                <  1:4317157‑M2/3‑1 (MQ=255)
  aaaccgggtgcgtcttaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                <  1:3812322‑M2/3‑1 (MQ=255)
    aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                              >  1:2289157‑M2/61‑65 (MQ=255)
    aaccggtgcgtcctaatgtagcccttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                              >  1:822092‑M2/61‑65 (MQ=255)
     accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                             >  1:2700702‑M2/60‑65 (MQ=255)
           gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                       >  1:4552352‑M2/54‑65 (MQ=255)
            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                      >  1:3938653/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:3496122/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:4455770/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:1243479/1‑65 (MQ=255)
               cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                   >  1:932641‑M2/50‑65 (MQ=255)
                ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                  >  1:1827523‑M2/49‑65 (MQ=255)
                 taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                 <  1:944128‑M2/18‑1 (MQ=255)
                        gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          <  1:3878255‑M2/25‑1 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATGACaa                                         >  1:592010‑M2/40‑65 (MQ=255)
                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3728884‑M2/40‑65 (MQ=255)
                          cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        >  1:1255207‑M2/39‑65 (MQ=255)
                           gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                       <  1:4543343‑M2/28‑1 (MQ=255)
                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                  <  1:3022170‑M2/33‑1 (MQ=255)
                                 gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                 <  1:406551‑M2/34‑1 (MQ=255)
                                  cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:2528736/65‑1 (MQ=255)
                                   attgagttttgctctttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:3190204‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:638176‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:2456439‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:3961823‑M2/30‑65 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:75627‑M2/30‑65 (MQ=255)
                                   aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:2029359/1‑65 (MQ=255)
                                   aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:1996604/1‑65 (MQ=255)
                                     tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             <  1:3136354‑M2/38‑1 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:4379717/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:4429324/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:349730/1‑65 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1081690‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3695115‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3595832‑M2/41‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:1053479‑M2/25‑65 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:1928940‑M2/25‑65 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2612377‑M2/41‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:4517093/65‑1 (MQ=255)
                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                         >  1:1891486/1‑65 (MQ=255)
                                          tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:2439621‑M2/23‑65 (MQ=255)
                                          tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:897255‑M2/23‑65 (MQ=255)
                                          tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:3186072‑M2/23‑65 (MQ=255)
                                           cAGAGCCTTCACCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCt                       >  1:4209862/1‑65 (MQ=255)
                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:864755/1‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAGGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2564752‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:3927966‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:116325‑M2/46‑1 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:352538‑M2/46‑1 (MQ=255)
                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1596411‑M2/18‑65 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1169357/1‑65 (MQ=255)
                                                ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:573428‑M2/17‑65 (MQ=255)
                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           >  1:1181244‑M2/10‑65 (MQ=255)
                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           >  1:4422079‑M2/10‑65 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:4265778‑M2/59‑1 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:4358646‑M2/59‑1 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:3127118‑M2/59‑1 (MQ=255)
                                                           accagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:441227‑M2/60‑1 (MQ=255)
                                                           accaaCTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGgaga       <  1:4124257‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:825399‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1863085‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1435742‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1585559‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       >  1:500144‑M2/6‑65 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1708021‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1854589‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       >  1:2848345‑M2/6‑65 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:320656‑M2/60‑1 (MQ=255)
                                                            ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt      >  1:316602‑M2/5‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:1579562‑M2/4‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:59150‑M2/4‑65 (MQ=255)
                                                             caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     >  1:3608305‑M2/4‑65 (MQ=255)
                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:899843/1‑65 (MQ=255)
                                                                |                                                                
GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683

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