|
0
|
1 <tool id="bcbio-nextgen" name="BCBio" version="0.1.0">
|
|
|
2 <command detect_errors="exit_code"><![CDATA[
|
|
3
|
3 PATH="\$BCBIO_DIRECTORY/tool_dir/bin:\$BCBIO_DIRECTORY/anaconda/bin:\$PATH" &&
|
|
0
|
4 mkdir $everything_else.files_path &&
|
|
|
5 cd $everything_else.files_path &&
|
|
|
6 python $__tool_directory__/script.py --workdir `pwd` --output-dir `pwd`/output --template $main --output-conf `pwd`/config.yaml
|
|
|
7 #if str( $region_file.has_region_file_select ) == "yes":
|
|
|
8 --bed ${region_file.bed}
|
|
|
9 #end if
|
|
|
10 #for $i, $sample in enumerate( $samples )
|
|
|
11 --input $sample.forward:$sample.reverse:$sample.phenotype
|
|
|
12 #end for
|
|
4
|
13 && bcbio_nextgen.py \$BCBIO_DIRECTORY/galaxy/bcbio_system.yaml `pwd`/config.yaml -t local -n \$BCBIO_CORES &&
|
|
0
|
14 zcat final/*/Batch1-ensemble-annotated.vcf.gz > $output_vcf &&
|
|
|
15 python $__tool_directory__/make_html.py --input-dir $everything_else.files_path --root-html $everything_else
|
|
|
16 ]]></command>
|
|
|
17 <configfiles>
|
|
|
18 <configfile name='main'>
|
|
|
19 details:
|
|
|
20 {% for sample in samples -%}
|
|
|
21 - algorithm:
|
|
|
22 aligner: $aligner
|
|
|
23 mark_duplicates: ${mark_duplicates}
|
|
|
24 remove_lcr: ${remove_lcr}
|
|
|
25 variantcaller: [$variantcallers]
|
|
|
26 {% if bed_file_path -%} variant_regions: {{ bed_file_path }}{% endif %}
|
|
|
27 ensemble:
|
|
|
28 numpass: 2
|
|
|
29 align_split_size: false
|
|
|
30 analysis: $analysis
|
|
|
31 lane: {{ loop.index }}
|
|
|
32 description: {{ sample['description'] }}
|
|
|
33 files: [{{ sample['forward'] }}, {{ sample['reverse'] }}]
|
|
|
34 genome_build: $build
|
|
|
35 metadata:
|
|
|
36 phenotype: {{ sample['phenotype'] }}
|
|
|
37 batch: Batch1
|
|
|
38 upload:
|
|
|
39 dir: ./final
|
|
|
40 {% endfor %}
|
|
|
41 </configfile>
|
|
|
42 </configfiles>
|
|
|
43 <inputs>
|
|
|
44 <param name="build" type="select" label="Genome Build">
|
|
|
45 <option value="GRCh37">GRCh37</option>
|
|
|
46 <option value="hg19">hg19</option>
|
|
|
47 <option value="GRCm38">GRCm38</option>
|
|
|
48 <option value="mm10">mm10</option>
|
|
|
49 </param>
|
|
|
50 <param name="analysis" type="select" label="Analysis">
|
|
|
51 <option value="variant2">variant2</option>
|
|
|
52 <option value="RNA-seq">RNA-seq</option>
|
|
|
53 <option value="smallRNA-seq">smallRNA-seq</option>
|
|
|
54 </param>
|
|
|
55 <param name="aligner" type="select" label="Aligner">
|
|
|
56 <option value="bwa">bwa</option>
|
|
|
57 <option value="bowtie">bowtie</option>
|
|
|
58 <option value="bowtie2">bowtie2</option>
|
|
|
59 <option value="hisat2">hisat2</option>
|
|
|
60 <option value="minimap2">minimap2</option>
|
|
|
61 <option value="novoalign">novoalign</option>
|
|
|
62 <option value="snap">snap</option>
|
|
|
63 <option value="star">star</option>
|
|
|
64 <option value="tophat2">tophat2</option>
|
|
|
65 <option value="false">false</option>
|
|
|
66 </param>
|
|
|
67 <param name="variantcallers" type="select" multiple="true" label="Variantcallers">
|
|
|
68 <option value="false">false</option>
|
|
|
69 <option value="freebayes">freebayes</option>
|
|
|
70 <option value="gatk-haplotype">gatk-haplotype</option>
|
|
|
71 <option value="haplotyper">haplotyper</option>
|
|
|
72 <option value="platypus">platypus</option>
|
|
|
73 <!--<option value="mutect">mutect</option>-->
|
|
|
74 <option value="mutect2">mutect2</option>
|
|
|
75 <option value="scalpel">scalpel</option>
|
|
|
76 <option value="tnhaplotyper">tnhaplotyper</option>
|
|
|
77 <option value="tnscope">tnscope</option>
|
|
|
78 <option value="vardict">vardict</option>
|
|
|
79 <option value="varscan">varscan</option>
|
|
|
80 <option value="samtools">samtools</option>
|
|
|
81 <option value="gatk">gatk</option>
|
|
|
82 </param>
|
|
|
83 <param name="mark_duplicates" type="select">
|
|
|
84 <option value="true" selected="true">true</option>
|
|
|
85 <option value="false">false</option>
|
|
|
86 </param>
|
|
|
87 <param name="remove_lcr" type="select">
|
|
|
88 <option value="true" selected="true">true</option>
|
|
|
89 <option value="false">false</option>
|
|
|
90 </param>
|
|
|
91 <repeat name="samples" title="Samples" min="1">
|
|
|
92 <param name="phenotype" type="select" label="Phenotype">
|
|
|
93 <option value="normal">Normal</option>
|
|
|
94 <option value="tumor">Tumor</option>
|
|
|
95 </param>
|
|
|
96 <param type="data" name="forward" format="fastq,fastq.gz,fastqsanger.gz" />
|
|
|
97 <param type="data" name="reverse" format="fastq,fastq.gz,fastqsanger.gz" />
|
|
|
98 </repeat>
|
|
|
99 <conditional name="region_file">
|
|
|
100 <param name="has_region_file_select" type="select" label="Do you have a region file">
|
|
|
101 <option value="yes">Yes</option>
|
|
|
102 <option value="no" selected="true">No</option>
|
|
|
103 </param>
|
|
|
104 <when value="yes">
|
|
|
105 <param type="data" name="bed" format="bed" />
|
|
|
106 </when>
|
|
|
107 </conditional>
|
|
|
108 </inputs>
|
|
|
109 <outputs>
|
|
|
110 <data name="output_vcf" format="vcf"/>
|
|
|
111 <data name="everything_else" format="html"/>
|
|
|
112 </outputs>
|
|
|
113 <help><![CDATA[
|
|
|
114 TODO: Fill in help.
|
|
|
115 ]]></help>
|
|
|
116 </tool> |
