Mercurial > repos > fubar > microsatbed

comparison microsatbed.xml @ 6:43d5f89e508c 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/microsatbed commit 7ceb6658309a7ababe622b5d92e729e5470e22f0-dirty
author fubar
date Sat, 13 Jul 2024 12:54:13 +0000
parents 1085e094cf5f
children
comparison
equal deleted inserted replaced
5:d54fa75d7d77 6:43d5f89e508c
1 <tool name="STR to bed" id="microsatbed" version="1.3.0" profile="22.05">
2 <description>Short Tandem Repeats to bed features from fasta</description>
3 <requirements>
4 <requirement version="3.12.3" type="package">python</requirement>
5 <requirement version="2.1.0" type="package">pyfastx</requirement>
6 <requirement version="1.3.0" type="package">pytrf</requirement>
7 </requirements>
8 <required_files>
9 <include path="find_str.py"/>
10 </required_files>
11 <version_command><![CDATA[python -c "import pytrf; from importlib.metadata import version; print(version('pytrf'))"]]></version_command>
12 <command><![CDATA[
13 #if $mode_cond.mode == "NATIVE":
14 pytrf findstr -f '$mode_cond.outformat' -o $bed -r $mode_cond.monomin $mode_cond.dimin $mode_cond.trimin $mode_cond.tetramin $mode_cond.pentamin $mode_cond.hexamin '${reference_genome.fasta}'
15 #else:
16 python '${__tool_directory__}/find_str.py'
17 --fasta '${reference_genome.fasta}'
18 --bed '$bed'
19 #if $mode_cond.mode == "SPECIFIC":
20 --specific '$mode_cond.specific'
21 --minreps '$mode_cond.minreps'
22 #else:
23 #if "MONO" in $mode_cond.subset:
24 --mono
25 #end if
26 #if "DI" in $mode_cond.subset:
27 --di
28 #end if
29 #if "TRI" in $mode_cond.subset:
30 --tri
31 #end if
32 #if "TETRA" in $mode_cond.subset:
33 --tetra
34 #end if
35 #if "PENTA" in $mode_cond.subset:
36 --penta
37 #end if
38 #if "HEXA" in $mode_cond.subset:
39 --hexa
40 #end if
41 --monomin '$mode_cond.monomin'
42 --dimin '$mode_cond.dimin'
43 --trimin '$mode_cond.trimin'
44 --tetramin '$mode_cond.tetramin'
45 --pentamin '$mode_cond.pentamin'
46 --hexamin '$mode_cond.hexamin'
47 #end if
48 #end if
49 ]]></command>
50 <inputs>
51 <conditional name="reference_genome">
52 <param name="genome_type_select" type="select" label="Select a source for fasta sequences to be searched for STRs" help="Options are to choose a built-in genome, or choose any history fasta file">
53 <option value="indexed">Use a Galaxy server built-in reference genome fasta</option>
54 <option selected="True" value="history">Use any fasta file from the current history</option>
55 </param>
56 <when value="indexed">
57 <param name="fasta" type="select" optional="false" multiple="false" label="Choose a built-in, or custom reference genome"
58 help="If the genome you need is not on the list, add a custom genome or choose a genome fasta from the current history">
59 <options from_data_table="all_fasta">
60 <filter column="2" type="sort_by"/>
61 <validator message="No genomes are available" type="no_options"/>
62 </options>
63 </param>
64 </when>
65 <when value="history">
66 <param name="fasta" type="data" format="fasta,fasta.gz" optional="false" multiple="false" label="Choose a fasta file from the current history"/>
67 </when>
68 </conditional>
69 <conditional name="mode_cond">
70 <param name="mode" type="select" label="Select patterns by motif length; or provide a specific motif pattern to report?" help="Choose *By length:* or *By pattern:* to configure STR selection mode">
71 <option selected="True" value="ALL">By length: Report all motifs of one or more specified lengths (1-6nt) as bed features</option>
72 <option value="SPECIFIC">By motif: Report one or more specific motifs (such as TCA,GC) as bed features</option>
73 <option value="NATIVE">All exact STR: use the pytrf findstr native command to a create csv, tsv or gtf output</option>
74 </param>
75 <when value="ALL">
76 <param name="subset" type="select" multiple="true" optional="false" label="Select at least 1 specific motif length to report"
77 help="Bed features will be output for every motif of the selected length(s) with the minimum required repeats or more">
78 <option value="DI" selected="true">All dimers (AC,AG,AT,...)</option>
79 <option value="TRI">All trimers (ACG,..)</option>
80 <option value="TETRA">All tetramers (ACGT,..)</option>
81 <option value="PENTA">All pentamers (ACGTC,..)</option>
82 <option value="HEXA">All hexamers (ACGTCG,..)</option>
83 <option value="MONO">All monomers (A,C...). Warning! Can produce overwhelming numbers of bed features</option>
84 </param>
85 <param name="dimin" type="integer" value="2" min="1" label="Minimum number of repeats for dimers"/>
86 <param name="trimin" type="integer" value="2" min="2" label="Minimum number of repeats for trimers"/>
87 <param name="tetramin" type="integer" value="2" min="2" label="Minimum number of repeats for tetramers"/>
88 <param name="pentamin" type="integer" value="2" min="2" label="Minimum number of repeats for pentamers"/>
89 <param name="hexamin" type="integer" value="2" min="2" label="Minimum number of repeats for hexamers"/>
90 <param name="monomin" type="integer" value="2" min="2" label="Minimum number of repeats for monomers"/>
91 </when>
92 <when value="SPECIFIC">
93 <param name="specific" type="text" label="Supply a specific motif pattern. Separate multiple patterns with commas such as GA,GC"
94 help="Make bed features only for the nominated specific motifs." optional="false"/>
95 <param name="minreps" type="integer" value="2" min="2" label="Minimum number of repeats of each of these motifs to report"
96 help="Short tandem repeats require 2 or more consecutive motifs by definition. WARNING: If monomers are included, every single matching base will be reported as a STR if minimum repeats = 1!"
97 optional="false"/>
98 </when>
99 <when value="NATIVE">
100 <param name="outformat" type="select" optional="false" label="Select the required output format"
101 help="Pytrf can create GFF, CSV or TSV output files">
102 <option value="gff" selected="true">GFF</option>
103 <option value="csv">Comma separated values</option>
104 <option value="tsv">Tab separated values</option>
105 </param>
106 <param name="dimin" type="integer" value="2" min="1" label="Minimum number of repeats for dimers"/>
107 <param name="trimin" type="integer" value="2" min="2" label="Minimum number of repeats for trimers"/>
108 <param name="tetramin" type="integer" value="2" min="2" label="Minimum number of repeats for tetramers"/>
109 <param name="pentamin" type="integer" value="2" min="2" label="Minimum number of repeats for pentamers"/>
110 <param name="hexamin" type="integer" value="2" min="2" label="Minimum number of repeats for hexamers"/>
111 <param name="monomin" type="integer" value="2" min="2" label="Minimum number of repeats for monomers"/>
112 </when>
113 </conditional>
114 </inputs>
115 <outputs>
116 <data name="bed" format="bed" label="STR from $fasta.element_identifier" hidden="false">
117 <change_format>
118 <when input="mode_cond.outformat" value="gff" format="gff" />
119 <when input="mode_cond.outformat" value="csv" format="csv" />
120 <when input="mode_cond.outformat" value="tsv" format="tabular" />
121 </change_format>
122 </data>
123 </outputs>
124 <tests>
125 <test expect_num_outputs="1">
126 <conditional name="reference_genome">
127 <param name="genome_type_select" value="history"/>
128 <param name="fasta" value="humsamp.fa"/>
129 </conditional>
130 <conditional name="mode_cond">
131 <param name="mode" value="ALL"/>
132 <param name="subset" value="DI,TRI,TETRA,PENTA,HEXA"/>
133 <param name="dimin" value="2"/>
134 <param name="trimin" value="2"/>
135 <param name="tetramin" value="2"/>
136 <param name="pentamin" value="2"/>
137 <param name="hexamin" value="2"/>
138 </conditional>
139 <output name="bed" value="bed_sample" compare="diff" lines_diff="0"/>
140 </test>
141 <test expect_num_outputs="1">
142 <conditional name="reference_genome">
143 <param name="genome_type_select" value="history"/>
144 <param name="fasta" value="humsamp.fa"/>
145 </conditional>
146 <conditional name="mode_cond">
147 <param name="mode" value="SPECIFIC"/>
148 <param name="specific" value="GC"/>
149 <param name="minreps" value="2"/>
150 </conditional>
151 <output name="bed" value="dibed_sample" compare="diff" lines_diff="0"/>
152 </test>
153 <test expect_num_outputs="1">
154 <conditional name="reference_genome">
155 <param name="genome_type_select" value="history"/>
156 <param name="fasta" value="mouse.fa"/>
157 </conditional>
158 <conditional name="mode_cond">
159 <param name="mode" value="NATIVE"/>
160 <param name="outformat" value="gff"/>
161 <param name="monomin" value="20"/>
162 <param name="dimin" value="10"/>
163 <param name="trimin" value="5"/>
164 <param name="tetramin" value="4"/>
165 <param name="pentamin" value="4"/>
166 <param name="hexamin" value="2"/>
167 </conditional>
168 <output name="bed" value="nativegff_sample" compare="diff" lines_diff="0"/>
169 </test>
170 </tests>
171 <help><![CDATA[
172
173 **Convert short repetitive sequences to bed features**
174
175 Microsatellites are usually defined as repeated short DNA patterns in an unbroken sequence.
176 A microsatellite pattern or *motif* can be any combination nucleotides, typically from 1 to 6nt in length.
177
178 This tool allows microsatellite and related features to be selected from a fasta sequence input file, and output into a single bed track, suitable for viewing in a genome browser such as JBrowse2.
179
180 All motifs of selected lengths can be reported as individual features in the output bed file, or specific motifs can be provided and all
181 others will be ignored. In all cases, a minimum required number of repeats can be specified. For example, requiring 2 or more repeats of the trimer *ACG* will report
182 every sequence of *ACGACG* or *ACGACGACG* or *ACGACGACGACG* and so on, as individual bed features. Similarly, requiring 3 repeats of any trimer will
183 report every distinct 3 nucleotide pattern, including *ACGACGACG* as well as every other unique 3 nucleotide pattern with 3 sequential repeats or more such, as "CTCCTCCTC*.
184
185 For other output formats, the pytrf native command line *findstr* can be used to produce a gff, csv or tsv output containing all exact short tandem repeats, as
186 described at the end of https://pytrf.readthedocs.io/en/latest
187
188 A fasta file must be supplied for processing. A built in genome can be selected, or a fasta file of any kind can be selected from the current history. Note that all
189 symbols are treated as valid nucleotides by pytrf, so extraneous characters such as *-* or *N* in the input fasta may appear as unexpected bed features. Lower case fasta symbols will be converted
190 to uppercase, to prevent them being reported as distinct motifs.
191
192
193 **Select motifs by length**
194
195 The default tool form setting is to select all dimer motif patterns.
196
197 Additional motif lengths from 1 to 6nt can be selected in the multiple-select drop-down list. All features will be returned in a single bed file. For each selected motif length,
198 the minimum number of repeats required for reporting can be adjusted. **Tandem repeats** are defined as at least 2 of any pattern. This tool allows singleton motifs to be reported,
199 so is not restricted to short tandem repeats (STR)
200
201 **Select motifs by pattern**
202
203 This option allows a motif pattern to be specified as a text string such as *CG* or *ATC*. Multiple motifs can be specified as a comma separated string such as *CG,ATC*.
204 All features will be returned as a single bed file.
205
206 The minimum number of repeats for all motifs can be set to match specific requirements.
207
208 For example, technical sequencing read bias may be influenced by the density of specific dimers, whether they are repeated or not
209 such as in https://github.com/arangrhie/T2T-Polish/tree/master/pattern
210
211 **Select all perfect STR using pytrf findstr in csv, tsv or gff output format**
212
213 This selection runs the pytrf *findstr* option to create gff/csv/tsv outputs as described at the end of https://pytrf.readthedocs.io/en/latest/.
214
215 Quoted here:
216
217 *A Tandem repeat (TR) in genomic sequence is a set of adjacent short DNA sequence repeated consecutively. The core sequence or repeat unit is generally called motif.
218 According to the motif length, tandem repeats can be classified as microsatellites and minisatellites. Microsatellites are also known as simple sequence repeats (SSRs)
219 or short tandem repeats (STRs) with motif length of 1-6 bp. Minisatellites are also sometimes referred to as variable number of tandem repeats (VNTRs) has longer motif length than microsatellites.
220 Pytrf is a lightweight Python C extension for identification of tandem repeats. The pytrf enables to fastly identify both exact or perfect SSRs.
221 It also can find generic tandem repeats with any size of motif, such as with maximum motif length of 100 bp. Additionally, it has capability of finding approximate or imperfect tandem repeats*
222
223 ]]></help>
224 <citations>
225 <citation type="bibtex">@misc{pytrf,
226 title = {{pytrf} Short tandem repeat finder, Accessed on July 10 2024},
227 howpublished = {\url{https://github.com/lmdu/pytrf}},
228 note = {Accessed on July 10 2024}
229 }</citation>
230 </citations>
231 </tool>
232