0
|
1 <tool id="gemini_recessive_and_dominant" name="GEMINI autosomal recessive/dominant" version="@VERSION@.0">
|
|
2 <description>Find variants meeting an autosomal recessive/dominant model</description>
|
|
3 <expand macro="requirements" />
|
|
4 <expand macro="version_command" />
|
|
5 <macros>
|
|
6 <import>gemini_macros.xml</import>
|
|
7 </macros>
|
|
8 <command>
|
|
9 <![CDATA[
|
|
10 gemini
|
|
11
|
2
|
12 #if str($rec_or_dom) == 'recessive':
|
0
|
13 ## start autosomal_recessive
|
|
14 autosomal_recessive
|
|
15 #else:
|
|
16 ## start autosomal_dominant
|
|
17 autosomal_dominant
|
|
18 #end if
|
|
19
|
2
|
20 #if str($report.report_selector) != 'all':
|
0
|
21 --columns "${report.columns}"
|
|
22 #end if
|
|
23
|
4
|
24 @CMDLN_SQL_FILTER_FILTER_OPTION@
|
0
|
25
|
|
26 -d $d
|
|
27 #if int($min_kindreds) > 0:
|
|
28 --min-kindreds $min_kindreds
|
|
29 #end if
|
|
30
|
|
31 "${ infile }"
|
|
32 > "${ outfile }"
|
|
33 ]]>
|
|
34 </command>
|
|
35 <expand macro="stdio" />
|
|
36 <inputs>
|
|
37
|
|
38 <param name="rec_or_dom" type="select" label="Autosomal ..." help="">
|
|
39 <option value="recessive">recessive</option>
|
|
40 <option value="dominant">dominant</option>
|
|
41 </param>
|
|
42
|
2
|
43 <expand macro="infile" />
|
0
|
44 <expand macro="column_filter" />
|
|
45 <expand macro="filter" />
|
|
46 <expand macro="min_sequence_depth" />
|
2
|
47 <param name="min_kindreds" size="4" type="integer" value="-1" label="The min. number of kindreds that must have a candidate variant in a gene"
|
0
|
48 help="-1 means default values (--min-kindreds)" />
|
|
49
|
|
50 </inputs>
|
|
51 <outputs>
|
2
|
52 <data name="outfile" format="tabular" />
|
0
|
53 </outputs>
|
|
54 <tests>
|
|
55 <test>
|
|
56 </test>
|
|
57 </tests>
|
|
58 <help>
|
|
59 **What it does**
|
|
60
|
|
61 Assuming you have defined the familial relationships between samples when loading your VCF into GEMINI, one can leverage a
|
|
62 built-in tool for identifying variants that meet an autosomal recessive or dominant inheritance pattern.
|
|
63 The reported variants will be restricted to those variants having the potential to impact the function of affecting protein coding transcripts.
|
|
64
|
|
65 @CITATION@
|
|
66 </help>
|
|
67 <expand macro="citations"/>
|
|
68 </tool>
|