Mercurial > repos > saskia-hiltemann > cgatools_v17
comparison tools/cgatools17/junctiondiff_v17.xml @ 1:3a2e0f376f26 draft
Minor change to tv2vcf.xml to allow for workflow automation
author | dgdekoning |
---|---|
date | Wed, 21 Oct 2015 10:09:15 -0400 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
0:751b62d30ae1 | 1:3a2e0f376f26 |
---|---|
1 <tool id="cg_junctiondiff" name="JunctionDiff" version="1.7.1"> | |
2 | |
3 <description>Reports difference between junction calls</description> <!--adds description in toolbar--> | |
4 | |
5 <requirements> | |
6 <requirement type="package" version="1">cgatools17</requirement> | |
7 </requirements> | |
8 | |
9 <command > <!--run executable--> | |
10 cgatools | head -1; | |
11 cgatools junctiondiff --beta | |
12 --reference ${crr.fields.crr_path} | |
13 --junctionsA $inputA | |
14 --junctionsB $inputB | |
15 --scoreThresholdA $scoreA | |
16 --scoreThresholdB $scoreB | |
17 --distance $distance | |
18 --minlength $minlength | |
19 $stat | |
20 --output-prefix cg_ | |
21 ; | |
22 mv cg_diff-* cg_diff.tsv | |
23 </command> | |
24 | |
25 <inputs> | |
26 <!-- reference --> | |
27 <param name="crr" type="select" label="Reference Build"> | |
28 <options from_data_table="cg_anno_files" /> | |
29 </param> | |
30 | |
31 <!--input files--> | |
32 <param name="inputA" type="data" format="tabluar" label="Junction file A"/> | |
33 <param name="inputB" type="data" format="tabluar" label="Junction file B"/> | |
34 | |
35 <!--parameters--> | |
36 <param name="stat" type="select" label="Print input file stats"> | |
37 <option value="">no</option> | |
38 <option value="--statout">yes</option> | |
39 </param> | |
40 | |
41 <param name="scoreA" type="integer" label="Score threshold value for input file A (default 10)" value="10"/> | |
42 <param name="scoreB" type="integer" label="Score threshold value for input file B (default 0)" value="0"/> | |
43 <param name="distance" type="integer" label="Max distance between coordinates of potentially compatible junctions (default 200)" value="200"/> | |
44 <param name="minlength" type="integer" label="Minimum deletion junction length to be included into the difference file (default 500)" value="500"/> | |
45 | |
46 <!-- prefix for output file so you dont have to manually rename history items --> | |
47 <param name="fname" type="text" value="" label="Prefix for your output file" help="Optional"/> | |
48 </inputs> | |
49 | |
50 <outputs> | |
51 <data format="tabular" name="corrected_junctions" from_work_dir="cg_diff.tsv" label="$fname ${tool.name} on ${on_string}: diff"/> | |
52 <data format="tabular" name="report" from_work_dir="cg_report.tsv" label="$fname ${tool.name} on ${on_string}: report"> | |
53 <filter>(stat == '--statout')</filter> | |
54 </data> | |
55 </outputs> | |
56 | |
57 | |
58 <help> | |
59 **What it does** | |
60 | |
61 This tool reports difference between junction calls of Complete Genomics junctions files | |
62 | |
63 **cgatools 1.7.1 Documentation** | |
64 | |
65 Userguide: http://cgatools.sourceforge.net/docs/1.7.1/cgatools-user-guide.pdf | |
66 | |
67 Release notes: http://cgatools.sourceforge.net/docs/1.7.1/cgatools-release-notes.pdf | |
68 | |
69 **Command line reference**:: | |
70 | |
71 COMMAND NAME | |
72 junctiondiff - Reports difference between junction calls of Complete Genomics junctions files. | |
73 | |
74 DESCRIPTION | |
75 junctiondiff takes two junction files A and B as input and produces the | |
76 following output: | |
77 - "diff-inputFileName" - the junctions from an input file A that are not | |
78 present in input file B. | |
79 - "report.txt" - a brief summary report (if --statout is used) | |
80 | |
81 Two junctions are considered equivalent if: | |
82 - they come from different files | |
83 - left and right positions of one junction are not more than "--distance" | |
84 bases apart from the corresponding positions of another junction | |
85 - the junction scores are equal or above the scoreThreshold | |
86 - they are on the same strands | |
87 | |
88 OPTIONS | |
89 -h [ --help ] | |
90 Print this help message. | |
91 | |
92 --beta | |
93 This is a beta command. To run this command, you must pass the --beta | |
94 flag. | |
95 | |
96 -s [ --reference ] arg | |
97 Reference file. | |
98 | |
99 -a [ --junctionsA ] arg | |
100 input junction file A. | |
101 | |
102 -b [ --junctionsB ] arg | |
103 input junction file B. | |
104 | |
105 -A [ --scoreThresholdA ] arg (=10) | |
106 score threshold value for the input file A. | |
107 | |
108 -B [ --scoreThresholdB ] arg (=0) | |
109 score threshold value for the input file B. | |
110 | |
111 -d [ --distance ] arg (=200) | |
112 Max distance between coordinates of potentially compatible junctions. | |
113 | |
114 -l [ --minlength ] arg (=500) | |
115 Minimum deletion junction length to be included into the difference | |
116 file. | |
117 | |
118 -o [ --output-prefix ] arg | |
119 The path prefix for all the output reports. | |
120 | |
121 -S [ --statout ] | |
122 (Debug) Report various input file statistics. Experimental feature. | |
123 | |
124 SUPPORTED FORMAT_VERSION | |
125 1.5 or later | |
126 </help> | |
127 </tool> |