46
|
1 <tool id="CompareOverlappingAdapt" name="compare overlapping">
|
|
2 <description>Provide the queries that overlap with a reference.</description>
|
|
3 <requirements>
|
|
4 <requirement type="set_environment">PYTHONPATH</requirement>
|
|
5 </requirements>
|
|
6 <command interpreter="python">
|
60
|
7 ../Java/Python/CompareOverlappingAdapt.py -i $formatType.inputFileName1 -f $formatType.FormatInputFileName1 -j $formatType2.inputFileName2 -g $formatType2.FormatInputFileName2 -o $outputFileGff $InvertMatch $NotOverlapping -m $NbOverlap
|
46
|
8 #if $OptionDistance.Dist == 'Yes':
|
|
9 -d $OptionDistance.distance
|
|
10 #end if
|
60
|
11 #if $OptionCA == 'Collinear':
|
46
|
12 -c
|
60
|
13 #elif $OptionCA == 'AntiSense':
|
46
|
14 -a
|
|
15 #end if
|
|
16 </command>
|
|
17
|
|
18 <inputs>
|
|
19 <conditional name="formatType">
|
|
20 <param name="FormatInputFileName1" type="select" label="Input Query File Format">
|
|
21 <option value="bed">bed</option>
|
|
22 <option value="gff">gff</option>
|
|
23 <option value="gff2">gff2</option>
|
|
24 <option value="gff3">gff3</option>
|
|
25 <option value="sam">sam</option>
|
|
26 <option value="gtf">gtf</option>
|
|
27 </param>
|
|
28 <when value="bed">
|
|
29 <param name="inputFileName1" format="bed" type="data" label="Input File 1"/>
|
|
30 </when>
|
|
31 <when value="gff">
|
|
32 <param name="inputFileName1" format="gff" type="data" label="Input File 1"/>
|
|
33 </when>
|
|
34 <when value="gff2">
|
|
35 <param name="inputFileName1" format="gff2" type="data" label="Input File 1"/>
|
|
36 </when>
|
|
37 <when value="gff3">
|
|
38 <param name="inputFileName1" format="gff3" type="data" label="Input File 1"/>
|
|
39 </when>
|
|
40 <when value="sam">
|
|
41 <param name="inputFileName1" format="sam" type="data" label="Input File 1"/>
|
|
42 </when>
|
|
43 <when value="gtf">
|
|
44 <param name="inputFileName1" format="gtf" type="data" label="Input File 1"/>
|
|
45 </when>
|
|
46 </conditional>
|
|
47
|
|
48 <conditional name="formatType2">
|
|
49 <param name="FormatInputFileName2" type="select" label="Input Reference File Format">
|
|
50 <option value="bed">bed</option>
|
|
51 <option value="gff">gff</option>
|
|
52 <option value="gff2">gff2</option>
|
|
53 <option value="gff3">gff3</option>
|
|
54 <option value="sam">sam</option>
|
|
55 <option value="gtf">gtf</option>
|
|
56 </param>
|
|
57 <when value="bed">
|
|
58 <param name="inputFileName2" format="bed" type="data" label="Input File 2"/>
|
|
59 </when>
|
|
60 <when value="gff">
|
|
61 <param name="inputFileName2" format="gff" type="data" label="Input File 2"/>
|
|
62 </when>
|
|
63 <when value="gff2">
|
|
64 <param name="inputFileName2" format="gff2" type="data" label="Input File 2"/>
|
|
65 </when>
|
|
66 <when value="gff3">
|
|
67 <param name="inputFileName2" format="gff3" type="data" label="Input File 2"/>
|
|
68 </when>
|
|
69 <when value="sam">
|
|
70 <param name="inputFileName2" format="sam" type="data" label="Input File 2"/>
|
|
71 </when>
|
|
72 <when value="gtf">
|
|
73 <param name="inputFileName2" format="gtf" type="data" label="Input File 2"/>
|
|
74 </when>
|
|
75 </conditional>
|
|
76 <conditional name="OptionDistance">
|
|
77 <param name="Dist" type="select" label="Maximum Distance between two reads">
|
|
78 <option value="Yes">Yes</option>
|
|
79 <option value="No" selected="true">No</option>
|
|
80 </param>
|
|
81 <when value="Yes">
|
|
82 <param name="distance" type="integer" value="0"/>
|
|
83 </when>
|
|
84 <when value="No">
|
|
85 </when>
|
|
86 </conditional>
|
60
|
87 <param name="OptionCA" type="select" label="Collinear or anti-sense features only">
|
|
88 <option value="Collinear">Collinear</option>
|
|
89 <option value="AntiSense">AntiSense</option>
|
|
90 <option value="All" selected="true">All</option>
|
|
91 </param>
|
46
|
92 <param name="InvertMatch" type="boolean" truevalue="-x" falsevalue="" checked="false" label="Invert match: the output file will contain all query elements which do NOT overlap"/>
|
|
93 <param name="NotOverlapping" type="boolean" truevalue="-O" falsevalue="" checked="false" label="Also report the query data which do not overlap, with the nbOverlaps tag set to 0."/>
|
60
|
94 <param name="NbOverlap" type="integer" value="1" label="Min. # of overlapping nt. to declare an overlap."/>
|
46
|
95 </inputs>
|
|
96
|
|
97 <outputs>
|
|
98 <data name="outputFileGff" format="gff3"/>
|
|
99 </outputs>
|
|
100
|
|
101 <help>
|
|
102 This script may be the most important one. It basically compares two sets of transcripts and keeps those from the first set which overlap with the second one. The first set is considered as the query set (basically, your data) and the second one is the reference set (RefSeq data, for example).
|
|
103
|
|
104 It is vital to understand that it will output the elements of the first file which overlap with the elements of the second one.
|
|
105
|
|
106 Various modifiers are also available:
|
|
107
|
|
108 -Invert selection (report those which do not overlap).
|
|
109
|
|
110 -Restrict to collinear / anti-sense overlapping data.
|
|
111
|
|
112 -Keep the query data even if they do not strictly overlap with the reference data, but are located not further away than *n* nucleotide from some reference data.
|
|
113
|
|
114 Some option reverses the selection. Put in other words, it performs the comparison as usual, and outputs all those query data which do not overlap.
|
|
115 </help>
|
|
116 </tool>
|