Mercurial > repos > xuebing > sharplabtool
diff tools/ilmn_pacbio/soap_denovo.xml @ 0:9071e359b9a3
Uploaded
author | xuebing |
---|---|
date | Fri, 09 Mar 2012 19:37:19 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tools/ilmn_pacbio/soap_denovo.xml Fri Mar 09 19:37:19 2012 -0500 @@ -0,0 +1,73 @@ +<tool id="soap_denovo" name="SOAPdenovo" version="1.0.0"> + <description>Short-read de novo assembly</description> + <!-- + # SOAPdenovo-127mer all -s ${soap_config} -o assembly -K ${k} -p 8 -d -D + # cat ${soap_config} > ${output1} + # cp ${soap_config} ${output1} && + --> + <command> + SOAPdenovo-127mer all -s ${soap_config} -o assembly -K ${k} -p 24 -d -D -R + </command> + <inputs> + <conditional name="inputs"> + <param name="read_type" type="select" label="Illumina read type"> + <option value="single">Single fragment</option> + <option value="paired">Paired-end</option> + </param> + <when value="single"> + <param name="input1" format="fastq" type="data" label="FASTQ file for reads"/> + </when> + <when value="paired"> + <param name="input1" format="fastq" type="data" label="FASTQ file for forward reads"/> + <param name="input2" format="fastq" type="data" label="FASTQ file for reverse reads"/> + <param name="d" type="integer" value="500" label="Estimated insert size for paired-end reads" /> + </when> + </conditional> + <param name="k" type="integer" value="23" label="Size of k for forming the de Bruijn overlap graph" /> + </inputs> + <configfiles> + <configfile name="soap_config">max_rd_len=105 +[LIB] +#if $inputs.read_type == "single" +q=${inputs.input1.file_name} +#else +avg_ins=${inputs.d} +asm_flags=3 +reverse_seq=0 +q1=${inputs.input1.file_name} +q2=${inputs.input2.file_name} +#end if + </configfile> + </configfiles> + <outputs> + <data name="assembled_contigs" format="fasta" from_work_dir="assembly.scafSeq" label="Assembled contigs from ${on_string}" /> + </outputs> + <help> + +**What it does** + +Runs SOAPdenovo_ to generate a genome assembly +using single-fragment or paired-end short reads. + +Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, Li S, Yang H, Wang J, Wang J. +"De novo assembly of human genomes with massively parallel short read sequencing." +*Genome Res.* 2010 Feb;20(2):265-72. + +.. _SOAPdenovo: http://soap.genomics.org.cn/soapdenovo.html + +**Parameter list** + +k + k-mer size for constructing the de Bruijn graph. The appropriate size of k is genome and data set dependent, but a good starting choice might be 75% of the read length. + +Insert size + For paired-end libraries, the expected insert size. + +**Output** + +assembly + + </help> +</tool> + +