New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 1209367 =5 (0.100)40 (1.180) 7/78 0.3 85.5% coding (290/630 nt) ycfK hypothetical protein
?W3110S.gb 1211196 = 10 (0.300)pseudogene (37/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment

TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCA  >  W3110S.gb/1211136‑1211249
                                                            |                                                     
aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                      >  1:12648‑M2/61‑62 (MQ=255)
aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                      >  1:520929‑M2/61‑62 (MQ=255)
aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                      >  1:3145097‑M2/61‑62 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaacgc                                                     <  1:3783684‑M2/3‑2 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:3474283‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:809804‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:801709‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:767749‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:3766984‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:3751924‑M2/3‑1 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:2125650‑M2/3‑1 (MQ=255)
     gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                 <  1:703760‑M2/7‑1 (MQ=255)
     gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                 <  1:429189‑M2/7‑1 (MQ=255)
     gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                 <  1:566907‑M2/7‑1 (MQ=255)
     gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                 <  1:1594804‑M2/7‑1 (MQ=255)
      aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                 >  1:3623850/1‑61 (MQ=255)
      aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                 >  1:2542252/1‑61 (MQ=255)
      aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                 >  1:2547027/1‑61 (MQ=255)
      aaGGCCACGGAACCAGGACACCAACAAGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                 >  1:472135/1‑61 (MQ=255)
        ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaa                                              >  1:1922880/1‑62 (MQ=255)
            aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                          <  1:160245/62‑1 (MQ=255)
            aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                          <  1:4197231/62‑1 (MQ=255)
            aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                          <  1:1576896/62‑1 (MQ=255)
            aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                          <  1:4257517/62‑1 (MQ=255)
            aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                          <  1:437666/62‑1 (MQ=255)
                   agccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                     <  1:4624623‑M2/3‑1 (MQ=255)
                    gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                  <  1:847664‑M2/22‑1 (MQ=255)
                    gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                  <  1:4099506‑M2/22‑1 (MQ=255)
                    gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                  <  1:61620‑M2/22‑1 (MQ=255)
                             gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                          <  1:73202‑M2/30‑1 (MQ=255)
                             gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                          <  1:4554597‑M2/30‑1 (MQ=255)
                             gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                          <  1:3136378‑M2/30‑1 (MQ=255)
                             gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                          <  1:2961523‑M2/30‑1 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                      >  1:809677‑M2/28‑61 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:2504581‑M2/28‑62 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:2499185‑M2/28‑62 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:54723‑M2/28‑62 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:1947947‑M2/28‑62 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:36672‑M2/28‑62 (MQ=255)
                                 tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                     >  1:1691430‑M2/28‑62 (MQ=255)
                                          gctccttctcccaaacaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:3770587‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:4305442‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:4134739‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:1727434‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:3641114‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:2149019‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:2434653‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:3328945‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:2613833‑M2/43‑1 (MQ=255)
                                          ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc             <  1:2727863‑M2/43‑1 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg           <  1:2569823‑M2/45‑1 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc          <  1:3813306‑M2/46‑1 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc          <  1:1414661‑M2/46‑1 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc          <  1:1431879‑M2/46‑1 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc          <  1:2944504‑M2/46‑1 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:2043956‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:2474488‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:2813604‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:4576510‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:3171073‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:3475523‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:3565999‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                           >  1:392344‑M2/15‑43 (MQ=255)
                                              ttctcccaaaccaaCCTATATGAAAATGAAGAAATa                                  >  1:3365287‑M2/15‑36 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                   >  1:3823186‑M2/9‑45 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa  >  1:2493132‑M2/9‑62 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa  >  1:1243612‑M2/9‑62 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa  >  1:1990154‑M2/9‑62 (MQ=255)
                                                    caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTACTGCTTTTa     >  1:2059251‑M2/9‑59 (MQ=255)
                                                            |                                                     
TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCA  >  W3110S.gb/1211136‑1211249

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