Mercurial > repos > iuc > colibread_takeabreak

comparison takeabreak.xml @ 0:e1fdaacdabc3 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/colibread commit a213833526146246d277ec7851165971449b501e
author iuc
date Fri, 20 Oct 2017 03:16:32 -0400
parents
children a018b8f663ea
comparison
equal deleted inserted replaced
-1:000000000000 0:e1fdaacdabc3
1 <?xml version='1.0' encoding='utf-8'?>
2 <tool profile="16.0" id="takeabreak" name="TakeABreak" version="1.1.2">
3 <description>is a tool that can detect inversion breakpoints directly from raw NGS reads</description>
4 <requirements>
5 <requirement type="package" version="1.1.2">takeabreak</requirement>
6 </requirements>
7 <macros>
8 <import>macros.xml</import>
9 </macros>
10 <command><![CDATA[
11
12 ## if single reads
13 #if str( $input_type_option.input_type ) == "simple"
14 #for $input in $input_type_option.reads
15 #if $input
16 #set $filename = str($input)+ "." + $input.ext
17 ln -sf '${input}' '${filename}' &&
18 echo "${filename}" >> input.fof &&
19 #end if
20 #end for
21 #else
22 ## if paired reads in a list
23 #for $i, $list in enumerate( $input_type_option.reads_lists )
24 #for $read in $list.list_reads
25 #set $filename = str($read) + "." + $read.ext
26 ln -sf '${read}' '${filename}' &&
27 echo "${filename}" >> "indiv_${i}.fof" &&
28 #end for
29 echo "indiv_${i}.fof" >> input.fof &&
30 #end for
31 #end if
32
33 TakeABreak
34
35 -in input.fof
36 -kmer-size ${kmer_size}
37 -abundance-min '${abundance_min}'
38 -abundance-max ${abundance_max}
39 -solidity-kind ${solidity_kind}
40 -lct ${lct}
41 -max-sim ${max_sim}
42 -repeat ${repeat}
43
44 ]]></command>
45
46 <inputs>
47 <!-- Input data files -->
48 <conditional name="input_type_option">
49 <param name="input_type" type="select" label="Input option">
50 <option value="simple">one individual = one read </option>
51 <option value="list">one individual = several reads</option>
52 </param>
53 <when value="simple">
54 <param name='reads' argument="-in" format="fasta,fastq" type="data" multiple="true" label="Single read files" optional="True" />
55 </when>
56 <when value="list">
57 <repeat name="reads_lists" title="Individual" min="1">
58 <param name='list_reads' argument="-in" format="fasta,fastq" type="data_collection" collection_type="list" multiple="true" label="Individual read files list"/>
59 </repeat>
60 </when>
61 </conditional>
62 <param name="kmer_size" argument="-kmer-size" type="integer" label="Length of the k-mers" value="31" help="Set the length of used kmers" />
63 <param name="abundance_min" argument="-abundance-min" type="text" label="Minimal abundance threshold" value="auto" />
64 <param name="abundance_max" argument="-abundance-max" type="integer" label="Maximal abundance threshold" value="2147483647" />
65
66 <param name="solidity_kind" argument="-solidity-kind" type="select" label="Way to consider a solid kmer with several input datasets" >
67 <option value="one">One</option>
68 <option value="sum">Sum</option>
69 <option value="all">All</option>
70 </param>
71
72 <param argument="-lct" type="integer" label="local complexity threshold" value="100" help="Local complexity threshold" />
73 <param name="max_sim" argument="-max-sim" type="integer" label="max similarity percentage" value="80" help="Inversions with a and b' (or u and v') whose longest common subsequence size is bigger than k*(this value)/100 are discarded" />
74 <param argument="-repeat" type="integer" label="maximal repeat size" value="8" help="Maximal repeat size at the breakpoint. Longest common suffix between a and b" />
75 </inputs>
76
77 <outputs>
78 <data format="fasta" name="inversion" from_work_dir="TakeABreak_*.fasta" label="Inversions with ${tool.name} on ${on_string}"/>
79 </outputs>
80 <tests>
81 <test>
82 <conditional name="input_type_option">
83 <param name="input_type" value="list" />
84 <repeat name="reads_lists">
85 <param name="list_reads">
86 <collection type="list">
87 <element name="example_reads" value="takeabreak/toy_example_reads.fasta" ftype="fasta"/>
88 <element name="example_reads_with_inv" value="takeabreak/toy_example_with_inv_reads.fasta" ftype="fasta"/>
89 </collection>
90 </param>
91 </repeat>
92 </conditional>
93 <output name="inversion" file="takeabreak/inversion_list.fasta"/>
94 </test>
95 <test>
96 <conditional name="input_type_option">
97 <param name="input_type" value="simple" />
98 <param name="reads" value="takeabreak/toy_example_reads.fasta,takeabreak/toy_example_with_inv_reads.fasta" />
99 </conditional>
100 <output name="inversion" file="takeabreak/inversion.fasta"/>
101 </test>
102 </tests>
103 <help><![CDATA[
104
105 **Description**
106
107 TakeABreak is a tool that can detect inversion breakpoints directly from raw NGS reads, without the need of any reference genome and without de novo assembling the genomes. Its implementation has a very limited memory impact allowing its usage on common desktop computers and acceptable runtime (Illumina reads simulated at 2x40x coverage from human chromosome 22 can be treated in less than two hours, with less than 1GB of memory).
108
109 -------
110
111 **Web site**
112
113 http://colibread.inria.fr/takeabreak/
114
115
116 ]]></help>
117 <expand macro="citations">
118 <citation type="doi">10.1007/978-3-319-07953-0_10</citation>
119 </expand>
120
121 </tool>