Mercurial > repos > iuc > fermikit_variants

comparison fermikit_run_calling.xml @ 0:96bf82df80c6 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/fermikit commit 16dcfc0fb84fad80fcf18417ae46c5499c96147a
author iuc
date Thu, 05 Jan 2017 08:35:17 -0500
parents
children 4c52b134a42f
comparison
equal deleted inserted replaced
-1:000000000000 0:96bf82df80c6
1 <tool id="fermikit_variants" name="fermikit-variants" version="0.14.dev1">
2 <description>call variants from genome-aligned contigs</description>
3 <requirements>
4 <requirement type="package" version="0.14.dev1">fermikit</requirement>
5 <requirement type="package" version="0.6.5">sambamba</requirement>
6 </requirements>
7 <command detect_errors="aggressive"><![CDATA[
8 #import re
9 #set escaped_element_identifier = re.sub('[^\w\-\s]', '_', str($bam.element_identifier))
10 #set ref = $reference_genome.fasta_item.fields.path if $reference_genome.reference_genome_source == "reference" else $reference_genome.history_item
11 export ROOT=\$(dirname \$(type -P k8)) &&
12 ln -f -s '$bam' '$escaped_element_identifier' &&
13 htsbox pileup -cuf '$ref' '$escaped_element_identifier' | gzip -1 > raw.vcf.gz &&
14 k8 "\$ROOT"/hapdip.js deovlp raw.vcf.gz | k8 "\$ROOT"/hapdip.js anno | gzip -1 > tmp.vcf.gz 2> flt.vcf.log &&
15 k8 "\$ROOT"/hapdip.js filter -q3 tmp.vcf.gz > flt.vcf 2>> flt.vcf.log &&
16 htsbox abreak -bcuf '$ref' <(sambamba sort -n -o /dev/stdout '$escaped_element_identifier') > sv.vcf
17 ]]></command>
18 <inputs>
19 <param name="bam" type="data" label="aligned contigs" help="To generate aligned contigs align fermi2 contigs with BWA mem options -B9 -O16 -L5" format="bam"/>
20 <conditional name="reference_genome" label="Reference genome to call variants against">
21 <param name="reference_genome_source" type="select">
22 <option value="reference">Use a built-in genome to call variants</option>
23 <option value="history">Use a genome from history to call variants</option>
24 </param>
25 <when value="history">
26 <param format="fasta" label="Select a reference genome" name="history_item" type="data" />
27 </when>
28 <when value="reference">
29 <param label="Select a reference genome" name="fasta_item" type="select">
30 <options from_data_table="all_fasta">
31 <filter column="2" type="sort_by"/>
32 </options>
33 </param>
34 </when>
35 </conditional>
36 </inputs>
37 <outputs>
38 <data name="structural_variants" format="vcf" label="fermikit SV on ${on_string}" from_work_dir="sv.vcf"/>
39 <data name="snps_indels" format="vcf" label="fermikit SNPs and short INDELs on ${on_string}" from_work_dir="flt.vcf"/>
40 </outputs>
41 <tests>
42 <test>
43 <param name="reference_genome_source" value="history"/>
44 <param name="history_item" value="small.fa.gz" ftype="fasta"/>
45 <param name="bam" value="aligned_contigs.bam"/>
46 <output name="structural_variants" file="sv.vcf" lines_diff="2"/>
47 <output name="snps_indels" file="flt.vcf" lines_diff="2"/>
48 </test>
49 </tests>
50 <help><![CDATA[
51
52 FermiKit is a de novo assembly based variant calling pipeline for deep Illumina
53 resequencing data. This galaxy wrapper can be used to call variants from contigs
54 generated by fermi2 that have subsequently been aligned to a reference genome
55 using bwa (options -B9 -O16 -L5 or -x intractg).
56
57 ]]></help> <citations> <citation
58 type="doi">10.1093/bioinformatics/btv440</citation> </citations> </tool>