|
48
|
1 <tool id="falco-call" name="falco-call" version="0.0.1">
|
|
|
2
|
|
|
3 <requirements>
|
|
|
4
|
|
|
5 <requirement type="package" version="0.1.19">samtools</requirement>
|
|
69
|
6 <requirement type="package" version="3.0.1">R</requirement>
|
|
86
|
7 <requirement type="set_environment">FALCO_PATH</requirement>
|
|
48
|
8
|
|
|
9 </requirements>
|
|
|
10
|
|
|
11 <description>FALCO: Amplicon Analysis Pipeline</description>
|
|
|
12
|
|
|
13 <command interpreter="bash">
|
|
86
|
14 falco-call.sh $falco_cfg \$FALCO_PATH;
|
|
48
|
15 </command>
|
|
|
16
|
|
|
17 <stdio>
|
|
|
18 <!-- Anything higher than 0 means the bash script didnt finish -->
|
|
|
19 <exit_code range="1:" level="fatal" description="Bash script didnt finish correctly, check log" />
|
|
|
20 </stdio>
|
|
|
21
|
|
|
22 <inputs>
|
|
|
23 <!-- ==================== -->
|
|
|
24 <!-- General inputs -->
|
|
|
25 <!-- ==================== -->
|
|
90
|
26 <param name="jobName" type="text" optional="false" label="Analysis/ouput name" help="Supply a name for the outputs to remind you what they contain. The string entered here will also be the sample-name in the VCF" value="TEST">
|
|
|
27 <validator type="empty_field" />
|
|
|
28 <validator type="regex" message="This field may contain only non-whitespace characters">\S+</validator>
|
|
|
29 </param>
|
|
|
30
|
|
48
|
31 <conditional name="genomeSource">
|
|
|
32 <param name="refGenomeSource" type="select" label="You can select a reference genome from your history or use a built-in index">
|
|
|
33 <option value="indexed">Use a built-in index</option>
|
|
|
34 <option value="history">Use one from history</option>
|
|
|
35 </param>
|
|
|
36 <when value="indexed">
|
|
|
37 <param name="indices" type="select" label="Select a reference genome">
|
|
|
38 <options from_data_table="all_fasta">
|
|
|
39 <filter type="sort_by" column="2" />
|
|
|
40 <validator type="no_options" message="No indexes are available in table" />
|
|
|
41 </options>
|
|
|
42 </param>
|
|
|
43 </when>
|
|
|
44 <when value="history">
|
|
|
45 <param name="ownFile" type="data" format="fasta" metadata_name="dbkey" label="Select a reference from history" />
|
|
|
46 </when>
|
|
|
47 </conditional>
|
|
|
48
|
|
|
49 <param name="bam" type="data" multiple="false" optional="false" format="bam" label="Input BAM" help="Select BAM file" />
|
|
|
50
|
|
90
|
51 <param format="tabular" name="filter_file" optional="True" type="data" label="FILTER File" help="When left empty, built-in list of genomic positions is used"/>
|
|
|
52 <param format="txt" name="manifest_file" optional="True" type="data" label="MANIFEST File" help="When left empty, built-in panel is used (TruSeq Amplicon Cancer Panel)"/>
|
|
48
|
53
|
|
|
54 <!-- ==================== -->
|
|
90
|
55 <!-- Optional hisotry output: copynumber TXT file -->
|
|
48
|
56 <!-- ==================== -->
|
|
90
|
57 <param name="txt2history" type="select" label="Also output falco specific TXT files to history" help="These text files are much larger but are required for the falco-report tool, so usually you only want these in a workflow and auto-delete them at the end">
|
|
|
58 <option value="false">no</option>
|
|
|
59 <option value="true">yes</option>
|
|
|
60 </param>
|
|
|
61
|
|
|
62 <!-- ==================== -->
|
|
|
63 <!-- Option to perform a test run with built in data -->
|
|
|
64 <!-- ==================== -->
|
|
|
65 <param name="testRun" type="select" label="Run with test data" help="Use built-in HCT116 data instead of input BAMs. This dataset contains only one amplicon.">
|
|
|
66 <option value="false">no</option>
|
|
|
67 <option value="true">yes</option>
|
|
|
68 </param>
|
|
48
|
69
|
|
|
70 </inputs>
|
|
|
71
|
|
|
72 <!-- ==================== -->
|
|
|
73 <!-- This config is sourced by tool -->
|
|
|
74 <!-- ==================== -->
|
|
|
75 <configfiles>
|
|
|
76 <configfile name="falco_cfg">
|
|
|
77 filter_file=$filter_file
|
|
|
78 manifest_file=$manifest_file
|
|
|
79 bam_file=$bam
|
|
|
80 bam_name=$bam.name
|
|
|
81 html_out=$html_output
|
|
|
82 vcf_out=$vcf_output
|
|
59
|
83 qc_ann_qual_out=$qc_ann_qual_output
|
|
|
84 qc2_ann_txt_out=$qc2_ann_txt_output
|
|
|
85 qc_targets_txt_out=$qc_targets_txt_output
|
|
90
|
86 also_output_txt=$txt2history
|
|
|
87 test_run=$testRun
|
|
59
|
88
|
|
48
|
89 out_path=${html_output.files_path}
|
|
58
|
90 job_name=$jobName
|
|
48
|
91 ## reference source
|
|
|
92 REF_SOURCE="${genomeSource.refGenomeSource}"
|
|
|
93 #if $genomeSource.refGenomeSource == "history":
|
|
|
94 ##build index on the fly
|
|
|
95 REF_FILE="${genomeSource.ownFile}"
|
|
|
96 DB_KEY="${dbkey}"
|
|
|
97 #else:
|
|
|
98 ##use precomputed indexes
|
|
52
|
99 REF_FILE="${genomeSource.indices.fields.path}"
|
|
48
|
100 #end if
|
|
|
101 </configfile>
|
|
|
102 </configfiles>
|
|
|
103
|
|
|
104 <!-- ==================== -->
|
|
|
105 <!-- Main output is an html based report, additional on request -->
|
|
|
106 <!-- ==================== -->
|
|
|
107 <outputs>
|
|
90
|
108 <data format="html" name="html_output" label="FALCO-calling (${jobName})">
|
|
|
109 <filter>(txt2history == 'true')</filter>
|
|
|
110 </data>
|
|
48
|
111 <data format="vcf" name="vcf_output" label="FALCO-calling (${jobName}): VCF" />
|
|
90
|
112 <data format="tabular" name="qc_ann_qual_output" label="FALCO-calling (${jobName}): qc_ann_qual">
|
|
|
113 <filter>(txt2history == 'true')</filter>
|
|
|
114 </data>
|
|
|
115 <data format="tabular" name="qc2_ann_txt_output" label="FALCO-calling (${jobName}): qc2_ann">
|
|
|
116 <filter>(txt2history == 'true')</filter>
|
|
|
117 </data>
|
|
|
118 <data format="tabular" name="qc_targets_txt_output" label="FALCO-calling (${jobName}): qc_targets">
|
|
|
119 <filter>(txt2history == 'true')</filter>
|
|
|
120 </data>
|
|
48
|
121 </outputs>
|
|
|
122
|
|
|
123 <!-- ==================== -->
|
|
|
124 <!-- Tests still to be done -->
|
|
|
125 <!-- ==================== -->
|
|
|
126
|
|
|
127 <!--
|
|
|
128 <tests>
|
|
|
129 <test>
|
|
|
130 <param name="input1" value="input1" />
|
|
|
131 <param name="input2" value="input2" />
|
|
|
132 </test>
|
|
|
133 </tests>
|
|
|
134 -->
|
|
|
135
|
|
|
136 <help>
|
|
|
137 .. class:: infomark
|
|
|
138
|
|
|
139 **Introduction**
|
|
|
140
|
|
|
141 This tool is a wrapper for the variant caller FALCO_ which is part of the Amplicon Analysis Pipeline (AAP)
|
|
|
142
|
|
|
143 .. _FALCO: https://github.com/tgac-vumc/falco/
|
|
|
144
|
|
|
145 Calls and annotates genomic variants for each amplicon in a design.
|
|
|
146
|
|
|
147 **Notes**
|
|
|
148
|
|
90
|
149 Because each amplicon is analysed separately, when two amplicons overlap this can cause the same variant position to be present twice in the output VCF. This is intentional, you can use this to evaluate the quality of the variant call beyond the amplicon reads.
|
|
48
|
150
|
|
|
151 -----
|
|
|
152
|
|
|
153 **Citation**
|
|
|
154
|
|
90
|
155 For the underlying tool please cite: Daoud Sie et al. Performance of amplicon-based next generation DNA sequencing for diagnostic gene mutation profiling in oncopathology (Cell Oncol 2014 Oct;37(5):353-61). Pubmed_
|
|
48
|
156
|
|
|
157 .. _Pubmed: http://www.ncbi.nlm.nih.gov/pubmed/25209392
|
|
|
158
|
|
|
159 </help>
|
|
|
160
|
|
|
161 </tool>
|
