New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 1209367 =12 (0.200)26 (0.650) 21/70 NT 71.0% coding (290/630 nt) ycfK hypothetical protein
?AP009048 1211196 = 13 (0.320)pseudogene (37/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment

TGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  AP009048/1211128‑1211266
                                                                    |                                                                      
agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                      >  1:5287655‑M2/69‑71 (MQ=255)
 gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                     <  1:4819033‑M2/4‑1 (MQ=255)
 gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                     <  1:1475001‑M2/4‑1 (MQ=255)
    ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                  >  1:4590590/1‑71 (MQ=255)
    atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                  >  1:3773606‑M2/65‑71 (MQ=255)
    atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                  >  1:4375770‑M2/65‑71 (MQ=255)
     tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                  >  1:4787840‑M2/64‑70 (MQ=255)
          gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                            >  1:4553394/1‑71 (MQ=255)
           cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                           <  1:2364650‑M2/14‑1 (MQ=255)
              tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                          <  1:3933631‑M2/15‑1 (MQ=255)
              tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                           >  1:1173106‑M2/55‑68 (MQ=255)
              aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                    <  1:2960083/59‑1 (MQ=255)
                cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                        >  1:2324482‑M2/53‑53 (MQ=255)
                  tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                             <  1:3061510‑M2/12‑1 (MQ=255)
                  gtcctaatgtagccgttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                    <  1:130584‑M2/21‑1 (MQ=255)
                  gtcctaatgtagccgttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                    <  1:124347‑M2/21‑1 (MQ=255)
                  ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaa                                                      <  1:5145182/69‑1 (MQ=255)
                    ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                  >  1:918112‑M2/49‑71 (MQ=255)
                     cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                 <  1:3824327/71‑1 (MQ=255)
                      aatgtaggcgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                             >  1:1771719‑M2/47‑58 (MQ=255)
                       atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                  <  1:2477998‑M2/23‑1 (MQ=255)
                        aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                              >  1:4425525/1‑71 (MQ=255)
                            gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATTAAGAAATAACAAGCaaa                                          <  1:2979889‑M2/31‑1 (MQ=255)
                             ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                         <  1:2297183‑M2/32‑1 (MQ=255)
                             ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                         <  1:3804611‑M2/32‑1 (MQ=255)
                              cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                         >  1:1015723‑M2/39‑70 (MQ=255)
                              cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                        <  1:1341765‑M2/33‑1 (MQ=255)
                               gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                       <  1:3964514‑M2/34‑1 (MQ=255)
                               gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                        >  1:361132‑M2/38‑70 (MQ=255)
                                cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                       >  1:4502948/1‑70 (MQ=255)
                                 acCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                     >  1:4824467/1‑71 (MQ=255)
                                  gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                     >  1:2252296‑M2/35‑70 (MQ=255)
                                     gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                           >  1:4585993‑M2/32‑61 (MQ=255)
                                      cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                >  1:1789640‑M2/31‑71 (MQ=255)
                                      cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                >  1:2525466‑M2/31‑71 (MQ=255)
                                        ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                  <  1:3231797‑M2/39‑1 (MQ=255)
                                        ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                               <  1:4610350/70‑1 (MQ=255)
                                          tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                            <  1:3032281/71‑1 (MQ=255)
                                          gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                            >  1:263277‑M2/27‑71 (MQ=255)
                                           aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaa                                                      >  1:2355332/1‑44 (MQ=255)
                                           aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                >  1:5296293/1‑50 (MQ=255)
                                            aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                     >  1:4876985/1‑44 (MQ=255)
                                             ttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                           <  1:5283440‑M2/14‑1 (MQ=255)
                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                       >  1:3262388/1‑40 (MQ=255)
                                               ttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                      >  1:519922‑M2/22‑40 (MQ=255)
                                               ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                   <  1:3522913‑M2/22‑1 (MQ=255)
                                               ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                        <  1:4801843‑M2/49‑1 (MQ=255)
                                               ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                        <  1:4334783‑M2/49‑1 (MQ=255)
                                                agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                      >  1:2746812/1‑71 (MQ=255)
                                                    ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                >  1:437519‑M2/17‑57 (MQ=255)
                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                      <  1:783392‑M2/35‑1 (MQ=255)
                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                           <  1:4356169‑M2/46‑1 (MQ=255)
                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                 >  1:1647094‑M2/16‑71 (MQ=255)
                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                 >  1:518603‑M2/16‑71 (MQ=255)
                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                          <  1:1467524‑M2/31‑1 (MQ=255)
                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                   <  1:1125566‑M2/38‑1 (MQ=255)
                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                         >  1:574951‑M2/15‑62 (MQ=255)
                                                         tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                  >  1:3925528‑M2/12‑50 (MQ=255)
                                                         tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                             >  1:4386790‑M2/12‑55 (MQ=255)
                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCt                             <  1:4698964‑M2/44‑1 (MQ=255)
                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                     <  1:3750735‑M2/52‑1 (MQ=255)
                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg              <  1:4712681‑M2/59‑1 (MQ=255)
                                                             aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt          >  1:4193697‑M2/8‑70 (MQ=255)
                                                             aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt          >  1:330629‑M2/8‑70 (MQ=255)
                                                              aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa        >  1:4734735‑M2/7‑71 (MQ=255)
                                                                ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa      >  1:4199143‑M2/5‑71 (MQ=255)
                                                                ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa      >  1:3989136‑M2/5‑71 (MQ=255)
                                                                ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa      >  1:2297229‑M2/5‑71 (MQ=255)
                                                                ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa      >  1:1853502‑M2/5‑71 (MQ=255)
                                                                ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa      <  1:950293‑M2/67‑1 (MQ=255)
                                                                 caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt     >  1:4434804‑M2/4‑71 (MQ=255)
                                                                   aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                       >  1:1243359‑M2/2‑51 (MQ=255)
                                                                    cGTTTATGAAAATGAAGAAATAACAAGCAAATGGct                                     >  1:2094524/1‑35 (MQ=255)
                                                                    cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt  >  1:2130319/1‑71 (MQ=255)
                                                                    cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt  >  1:469669/1‑71 (MQ=255)
                                                                    |                                                                      
TGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  AP009048/1211128‑1211266

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