New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =43 (0.970)13 (0.390) 10/102 0.9 24.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 46 (1.370)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209553‑1209680
                                                                  |                                                             
gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                              >  1:2985275/1‑68 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             >  1:94800‑M2/66‑68 (MQ=255)
  ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGtt                                                            >  1:1355101/1‑68 (MQ=255)
       accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                       >  1:245510‑M2/60‑68 (MQ=255)
       accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                       >  1:2389503‑M2/60‑68 (MQ=255)
             gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                 >  1:1235250‑M2/54‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:3047447/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:1988910/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:21787/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:2326464/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:2346354/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:2918330/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:3185163/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:844714/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  1:1259198/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTACCAAACCAACGTTTATGAAAATgaag                                               >  1:2722472/1‑68 (MQ=255)
                      tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:3094736‑M2/45‑68 (MQ=255)
                      tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:435688‑M2/45‑68 (MQ=255)
                       gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                       >  1:1143148‑M2/44‑68 (MQ=255)
                         cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                     <  1:1175577/68‑1 (MQ=255)
                         cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                     <  1:1871655/68‑1 (MQ=255)
                          gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                    <  1:111442/68‑1 (MQ=255)
                             acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACTTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:1596232/68‑1 (MQ=255)
                              cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                >  1:727499/1‑68 (MQ=255)
                                ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                              <  1:499627/68‑1 (MQ=255)
                                 cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  1:1869813/68‑1 (MQ=255)
                                 cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  1:571587/68‑1 (MQ=255)
                                     gggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAANAACAAGCAAATGGcatca                         <  1:2573481/67‑1 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACCAGCAAATGGcatca                         >  1:3469355/1‑68 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         >  1:3115811/1‑68 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         >  1:2352352/1‑68 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         <  1:1363783/68‑1 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         >  1:1397607/1‑68 (MQ=255)
                                      ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                        <  1:1540187‑M2/40‑1 (MQ=255)
                                      ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                        <  1:2733493/68‑1 (MQ=255)
                                      ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                        <  1:3209554/68‑1 (MQ=255)
                                       gTAGTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:2673516/1‑68 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:681815/68‑1 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:2507234/1‑68 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:3072067/1‑68 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:2397173/1‑68 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:199725/68‑1 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:510596/1‑68 (MQ=255)
                                          gttttgctccttctcccaaacccaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:2347881‑M2/25‑68 (MQ=255)
                                            tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:3353408‑M2/46‑1 (MQ=255)
                                            gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2409869/68‑1 (MQ=255)
                                            gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:662202/68‑1 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:3304633/1‑68 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:1331016/1‑68 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:3001020/1‑68 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:257992/1‑68 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:1869021/1‑68 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  1:1671618/1‑68 (MQ=255)
                                              agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:1339293/1‑68 (MQ=255)
                                                 tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             >  1:347703‑M2/18‑68 (MQ=255)
                                                 tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             <  1:1935099‑M2/51‑1 (MQ=255)
                                                 tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             >  1:863044‑M2/18‑68 (MQ=255)
                                                 gCCTTCTCCCAAACCAACGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             >  1:2506638/1‑68 (MQ=255)
                                                 gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             >  1:1148028/1‑68 (MQ=255)
                                                  ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc            >  1:1133897‑M2/17‑68 (MQ=255)
                                                  ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc            >  1:2236121‑M2/17‑68 (MQ=255)
                                                  ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc            >  1:2865145‑M2/17‑68 (MQ=255)
                                                   cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc           >  1:2502921‑M2/16‑68 (MQ=255)
                                                      ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTttaccctg        >  1:2806551‑M2/13‑65 (MQ=255)
                                                        cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg      >  1:345076‑M2/11‑68 (MQ=255)
                                                        cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg      >  1:2314026‑M2/11‑68 (MQ=255)
                                                        cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg      >  1:2024715‑M2/11‑68 (MQ=255)
                                                         ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:2463740‑M2/59‑1 (MQ=255)
                                                         ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:1645438‑M2/59‑1 (MQ=255)
                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1048010‑M2/60‑1 (MQ=255)
                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:170345‑M2/60‑1 (MQ=255)
                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1411816‑M2/60‑1 (MQ=255)
                                                            aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  >  1:545035‑M2/7‑68 (MQ=255)
                                                                  |                                                             
GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209553‑1209680

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.