Mercurial > repos > devteam > join
comparison gops_join.py @ 3:860e7e4899b1 draft
planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
author | devteam |
---|---|
date | Tue, 13 Oct 2015 12:50:54 -0400 |
parents | cb9ed67cb329 |
children | ddf5484243d2 |
comparison
equal
deleted
inserted
replaced
2:83c18f310619 | 3:860e7e4899b1 |
---|---|
6 -1, --cols1=N,N,N,N: Columns for start, end, strand in first file | 6 -1, --cols1=N,N,N,N: Columns for start, end, strand in first file |
7 -2, --cols2=N,N,N,N: Columns for start, end, strand in second file | 7 -2, --cols2=N,N,N,N: Columns for start, end, strand in second file |
8 -m, --mincols=N: Require this much overlap (default 1bp) | 8 -m, --mincols=N: Require this much overlap (default 1bp) |
9 -f, --fill=N: none, right, left, both | 9 -f, --fill=N: none, right, left, both |
10 """ | 10 """ |
11 import sys, traceback, fileinput | 11 import fileinput |
12 from warnings import warn | 12 import sys |
13 from bx.intervals import * | 13 from bx.intervals.io import NiceReaderWrapper |
14 from bx.intervals.io import * | 14 from bx.intervals.operations.join import join |
15 from bx.intervals.operations.join import * | |
16 from bx.cookbook import doc_optparse | 15 from bx.cookbook import doc_optparse |
17 from galaxy.tools.util.galaxyops import * | 16 from bx.tabular.io import ParseError |
17 from galaxy.tools.util.galaxyops import fail, parse_cols_arg, skipped | |
18 | 18 |
19 assert sys.version_info[:2] >= ( 2, 4 ) | 19 assert sys.version_info[:2] >= ( 2, 4 ) |
20 | 20 |
21 | |
21 def main(): | 22 def main(): |
22 mincols = 1 | 23 mincols = 1 |
23 upstream_pad = 0 | |
24 downstream_pad = 0 | |
25 leftfill = False | 24 leftfill = False |
26 rightfill = False | 25 rightfill = False |
27 | 26 |
28 options, args = doc_optparse.parse( __doc__ ) | 27 options, args = doc_optparse.parse( __doc__ ) |
29 try: | 28 try: |
30 chr_col_1, start_col_1, end_col_1, strand_col_1 = parse_cols_arg( options.cols1 ) | 29 chr_col_1, start_col_1, end_col_1, strand_col_1 = parse_cols_arg( options.cols1 ) |
31 chr_col_2, start_col_2, end_col_2, strand_col_2 = parse_cols_arg( options.cols2 ) | 30 chr_col_2, start_col_2, end_col_2, strand_col_2 = parse_cols_arg( options.cols2 ) |
32 if options.mincols: mincols = int( options.mincols ) | 31 if options.mincols: |
32 mincols = int( options.mincols ) | |
33 if options.fill: | 33 if options.fill: |
34 if options.fill == "both": | 34 if options.fill == "both": |
35 rightfill = leftfill = True | 35 rightfill = leftfill = True |
36 else: | 36 else: |
37 rightfill = options.fill == "right" | 37 rightfill = options.fill == "right" |