New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 120938213 (0.260)17 (0.500) 16/70 NT 68.1% coding (305/630 nt) ycfK hypothetical protein
?AP009048 = 1211179 7 (0.200)pseudogene (54/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAC  >  AP009048/1211115‑1211247
                                                                |                                                                    
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc                                                                 >  1:1135538‑M2/1‑65 (MQ=255)
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc                                                                 >  1:1140971‑M2/1‑65 (MQ=255)
   gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                              <  1:668617‑M2/70‑9 (MQ=255)
       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                         <  1:1745426‑M2/71‑14 (MQ=255)
         cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                        <  1:1383313‑M2/70‑15 (MQ=255)
             tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                   >  1:915811/1‑71 (MQ=255)
              ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                  <  1:4001002‑M2/71‑21 (MQ=255)
                 ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttttt                                               >  1:3925584/1‑69 (MQ=255)
                         tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                       <  1:2598576‑M2/71‑32 (MQ=255)
                         tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccca                                                            >  1:2090866‑M2/1‑40 (MQ=255)
                          gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagatccgt                                      <  1:3598600‑M2/71‑33 (MQ=255)
                          gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagatccgt                                      <  1:3593122‑M2/71‑33 (MQ=255)
                          gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                          <  1:3872‑M2/51‑13 (MQ=255)
                           aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                     <  1:3371156/71‑1 (MQ=255)
                           aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                       <  1:2483918/69‑1 (MQ=255)
                            aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                    <  1:1074324/71‑1 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                  <  1:503482/71‑1 (MQ=255)
                                  cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                              <  1:3006181‑M2/71‑41 (MQ=255)
                                     aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                           <  1:2212546‑M2/71‑44 (MQ=255)
                                        cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                        <  1:17067/71‑1 (MQ=255)
                                            acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                    >  1:2274415‑M2/1‑21 (MQ=255)
                                             cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgat                   <  1:2424059‑M2/71‑52 (MQ=255)
                                             cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                    >  1:2976591/1‑70 (MQ=255)
                                              acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                            <  1:646697‑M2/61‑43 (MQ=255)
                                                cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                             >  1:3523856‑M2/1‑17 (MQ=255)
                                                 aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                             <  1:3846395‑M2/41‑26 (MQ=255)
                                                  aCAGGTAATGCAGCGccttctcccaaaccaaggttttcgag                                            <  1:3374355‑M2/41‑27 (MQ=255)
                                                   cAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata             <  1:3297396‑M2/71‑58 (MQ=255)
                                                   cAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata             >  1:194466‑M2/1‑14 (MQ=255)
                                                   cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                      >  1:2157128/1‑46 (MQ=255)
                                                    aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat            <  1:3778663‑M2/71‑59 (MQ=255)
                                                       tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc                      >  1:2055616‑M2/1‑10 (MQ=255)
                                                         aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc        <  1:655554‑M2/70‑63 (MQ=255)
                                                           gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa     <  1:825630‑M2/71‑66 (MQ=255)
                                                            cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac    <  1:1892484‑M2/71‑67 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc            >  1:1564835/1‑63 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt     <  1:3382553/70‑1 (MQ=255)
                                                              gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg  >  1:682437‑M2/1‑3 (MQ=255)
                                                              gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg  >  1:3864788‑M2/1‑3 (MQ=255)
                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCa                         <  1:4304446/48‑1 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc  >  1:1861870/1‑69 (MQ=255)
                                                                |                                                                    
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAC  >  AP009048/1211115‑1211247

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.