|
0
|
1 <tool id="DimontDataExtractor" name="Dimont Data Extractor" version="0.1">
|
|
|
2 <description>prepares data in the format required by Dimont</description>
|
|
|
3 <command interpreter="perl">extract_data_single_galaxy.pl $genomefa $regions $chromcol $startcol $seccol $seccoord $width $statcol $extracted</command>
|
|
|
4 <inputs>
|
|
|
5 <param type="data" format="bed,gtf,txt,tabular" name="regions" label="Genomic regions" value="" optional="false" help="The genomic regions to be extracted in a BED-like file format, e.g., BED, GTF, narrowPeak." />
|
|
|
6 <param type="data" format="fasta" name="genomefa" label="Genome sequence" value="" optional="false" help="The input genome to which the genomic regions refer." />
|
|
|
7 <param type="data_column" data_ref="regions" force_select="true" name="chromcol" label="Chromosome column" value="1" optional="false" help="The column of the regions file, which contains the chromosome information." />
|
|
|
8 <param type="data_column" data_ref="regions" force_select="true" name="startcol" label="Start column" value="2" optional="false" help="The column of the Regions file containing the start position of the genomic region." />
|
|
|
9 <param type="select" name="seccol" label="Meaning of second coordinate" optional="false" help="The meaning of the second genomic coordinate. This may either be the position of the peak summit relative to the position in Start, or the end position of the peak.">
|
|
|
10 <option value="end">End of peak (in global coordinates)</option>
|
|
|
11 <option value="center">Center of peak (relative to start)</option>
|
|
|
12 </param>
|
|
|
13 <param type="data_column" data_ref="regions" force_select="true" name="seccoord" label="Second coordinate" value="3" optional="false" help="The second genomic coordinate with meaning specified by parameter "Meaning of second coordinate"" />
|
|
|
14 <param type="data_column" data_ref="regions" force_select="true" name="statcol" label="Statistics column" value="7" optional="false" help="The column containing the peak statistics information (or another measure of peak confidence)." />
|
|
|
15 <param type="integer" name="width" label="Width" value="1000" optional="false" help="The width of the genomic region to be extracted. Recommended values: 1000 for ChIP-seq and 100 for ChIP-exo." />
|
|
|
16 </inputs>
|
|
|
17 <outputs>
|
|
|
18 <data format="fasta" name="extracted">
|
|
|
19 </data>
|
|
|
20 </outputs>
|
|
|
21 <tests>
|
|
|
22 <test>
|
|
|
23 <param name="regions" value="mini.bed" />
|
|
|
24 <param name="genomefa" value="minigenome.fa" />
|
|
|
25 <param name="chromcol" value="1" />
|
|
|
26 <param name="startcol" value="2" />
|
|
|
27 <param name="seccol" value="end" />
|
|
|
28 <param name="seccoord" value="3" />
|
|
|
29 <param name="statcol" value="7" />
|
|
|
30 <param name="width" value="200" />
|
|
|
31 <output name="extracted" file="mini_extracted.fa" />
|
|
|
32 </test>
|
|
|
33 <test>
|
|
|
34 <param name="regions" value="mini2.bed" />
|
|
|
35 <param name="genomefa" value="minigenome.fa" />
|
|
|
36 <param name="chromcol" value="1" />
|
|
|
37 <param name="startcol" value="2" />
|
|
|
38 <param name="seccol" value="center" />
|
|
|
39 <param name="seccoord" value="3" />
|
|
|
40 <param name="statcol" value="7" />
|
|
|
41 <param name="width" value="200" />
|
|
|
42 <output name="extracted" file="mini2_extracted.fa" />
|
|
|
43 </test>
|
|
|
44 </tests>
|
|
|
45 <help>
|
|
|
46 **Dimont Data Extractor** prepares an annotated FastA file as required by Dimont from a genome (in FastA format) and a tabular file (e.g., BED, GTF, narrowPeak,...). The regions specified in the tabular file are used to determine the center of the extracted sequences. All extracted sequences have the same length as specified by parameter "Width".
|
|
|
47
|
|
|
48 In case of ChIP data, the center position could for instance be the peak summit.
|
|
|
49 An annotated FastA file for ChIP-exo data comprising sequences of length 100 centered around the peak summit might look like::
|
|
|
50
|
|
|
51 > peak: 50; signal: 515
|
|
|
52 ggccatgtgtatttttttaaatttccac...
|
|
|
53 > peak: 50; signal: 199
|
|
|
54 GGTCCCCTGGGAGGATGGGGACGTGCTG...
|
|
|
55 ...
|
|
|
56
|
|
|
57 where the center is given as 50 for the first two sequences, and the confidence amounts to 515 and 199, respectively.
|
|
|
58
|
|
|
59 We also provide an example_ input file and a stand alone Perl script_ for preparing data in the format required by Dimont_.
|
|
|
60
|
|
|
61
|
|
|
62 If you experience problems using Dimont Data Extractor, please contact_ us.
|
|
|
63
|
|
|
64 .. _example: http://www.jstacs.de/downloads/dimont-example.fa
|
|
|
65 .. _script: http://www.jstacs.de/index.php/Dimont#Data_preparation
|
|
|
66 .. _Dimont: http://jstacs.de/index.php/Dimont
|
|
|
67 .. _contact: mailto:grau@informatik.uni-halle.de
|
|
|
68 </help>
|
|
|
69 </tool>
|
|
|
70
|
