Mercurial > repos > spanish_national_institue_of_bioinformatics > nucleosome_dynamics

comparison nucleR.xml @ 2:e88c806ddf3e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author spanish_national_institue_of_bioinformatics
date Fri, 12 Apr 2019 05:28:43 -0400
parents
children
comparison
equal deleted inserted replaced
1:707216744193 2:e88c806ddf3e
1 <tool id="nucleR" name="NucleR" version="0.1">
2 <description>: analyze aligned BAM files from MNAse (in RData format)</description>
3 <requirements>
4 <requirement type="binary">docker</requirement>
5 </requirements>
6 <command>
7 <![CDATA[
8 docker run -v $__root_dir__/database/files:$__root_dir__/database/files -v /tmp:/tmp -u `id -u`:`id -g` mmbirb/nucleosome-dynamics nucleR --input $rdata_file --output $output_gff_file --type $seq_type --width $width --minoverlap $minoverlap --dyad_length $dyad_length --hthresh $hthres --wthresh $wthres --pcKeepComp $pcKeepComp --fdrOverAmp $fdrOverAmp --trim $trim --fragmentLen $fragmentLen
9 #if $threshold.is == "absolute":
10 --thresholdValue ${threshold.thresholdValue}
11 #else if $threshold.is == "percentage":
12 --thresholdPercentage ${threshold.thresholdPercentage}
13 #end if
14 #if $chr
15 --chr $chr
16 #end if
17 #if $start
18 --start $start
19 #end if
20 #if $end
21 --end $end
22 #end if
23 #if $comp
24 --components $comp
25 #end if
26 ]]>
27 </command>
28 <inputs>
29 <param name="rdata_file" type="data" format="rdata" label="Input MNase-seq/ATAC-seq reads (RData)" help="Input BAM file in RData format as generated by readBAM."/>
30 <param name="seq_type" type="select" label="Type of sequence reads">
31 <option value="paired">Paired</option>
32 <option value="single">Single</option>
33 </param>
34 <param name="width" size="4" type="integer" value="147" label="Width (bp)" help="Size of each nucleosome."/>
35 <param name="minoverlap" type="integer" value="80" label="Minimum Overlap (bp)" help="Minimum overlap between two nucleosomes for merging them."/>
36 <param name="dyad_length" type="integer" value="50" label="Dyad length (bp)" help="Lenght of the reads that should be used for nucleosome calling to define the dyad of the nucleosomes keeping the number of bases around the center of the read."/>
37 <param name="fragmentLen" type="integer" value="170" label="Fragment Length (bp)" help="Maximum fragment length allowed (bp)" optional="True" />
38 <param name="chr" type="text" value="" label="Chromosome" optional="True" help="Chromosome to consider for the analysis in the given input file."/>
39 <param name="start" type="integer" value="" label="Start" optional="True" help="Start genomic position to consider for the analysis in the given input file."/>
40 <param name="end" type="integer" value="" label="End" optional="True" help="End genomic position to consider for the analysis in the given input file."/>
41 <param name="comp" type="integer" value="1" label="Components" optional="True" help="Number of negative binomials that will be used to filter duplicated reads."/>
42 <conditional name="threshold">
43 <param name="is" type="select" label="Background threshold" help="Minimum number of reads (Coverage) to call a nucleosome. Can be given as a percentage, or an absolute number of reads.">
44 <option value="percentage" selected="True">use a percentage</option>
45 <option value="absolute">use a number of reads</option>
46 </param>
47 <when value="absolute">
48 <param name="thresholdValue" type="integer" value="10" label="Bakground level (reads)" help="Absolute value to filter out nucleosome calls. It is the minimum number of reads (coverage) in a nucleosome call expressed as reads per million of mapped reads."/>
49 </when>
50 <when value="percentage">
51 <param name="thresholdPercentage" type="float" value="35" label="Background level (%)" help="Percentile of coverage in the experiment used as threshold to filter out nucleosome calls (i.e., 25% would mean that only peaks with coverage in the 1st quantile would be considered)"/>
52 </when>
53 </conditional>
54 <param name="hthres" type="float" value="0.4" label="Height Threshold" help="Height threshold (between 0 and 1) to classify a nucleosome as fuzzy (class=F) or well-positioned ( class=W) according to the number of reads at the dyad. Nucleosomes below this value will be defined as fuzzy."/>
55 <param name="wthres" type="float" value="0.6" label="Width Threshold" help="Width threshold (between 0 and 1) to classify a nucleosome as fuzzy (class=F) or well-positioned (class=W) according to the dispersion of the reads around the dyad."/>
56 <param name="fdrOverAmp" type="float" value="0.05" label="fdrOverAmp" help="Threshold to filter over-amplified reads , as defined in filterDuplReads function of htSetqTools R package."/>
57 <param name="pcKeepComp" type="float" value="0.02" label="Coverage Smoothing" help="Parameter used in the smoothing when Fourier transformation is applied. Number of components to select with respect to the total size of the sample. Allowed values are numeric (in range 0:1) for manual setting, or auto for automatic detection."/>
58 <param name="trim" type="integer" value="50" label="Trim" help="Number of basepairs to keep from each read (or to extend in case it is larger than the read width)."/>
59 </inputs>
60 <outputs>
61 <data format="gff" name="output_gff_file" label="NR__${os.path.splitext(($rdata_file.name.split('__'))[1])[0]}.gff" />
62 </outputs>
63 <tests>
64 <test>
65 <param name="rdata_file" value="readBAM__cellcycleM_chrII.RData" />
66 <param name="seq_type" value="paired" />
67 <param name="width" value="147" />
68 <param name="minoverlap" value="80" />
69 <param name="dyad_length" value="50" />
70 <param name="thresholdPercentage" value="35" />
71 <param name="hthres" value="0.4" />
72 <param name="wthres" value="0.6" />
73 <param name="pcKeepComp" value="0.02" />
74 <output name="output_gff_file" file="NR__cellcycleM_chrII.gff" />
75 </test>
76 </tests>
77 <help>
78 .. image:: ${static_path}/images/NucleosomeDynamicsLogo.png
79 :height: 80
80 :width: 200
81
82 -----
83
84 Nucleosome Dynamics is a set of tools that take MNase-seq and ATAC-seq aligned reads and performs a serie of nucleosome-related analyses on them.
85
86 .. class:: infomark
87
88 Visit the documentation of the original application for learning more about the accepted values and formats. http://mmb.irbbarcelona.org/NucleosomeDynamics/help/usage/nucleosome-dynamics
89 </help>
90 <citations>
91 <citation type="bibtex">
92 @misc{github,
93 author = {Buitrago D},
94 year = {2019},
95 title = {Nucleosome Dynamics suite: containerized installation},
96 publisher = {Github},
97 journal = {GitHub repository},
98 url = {https://github.com/nucleosome-dynamics/nucleosome_dynamics},
99 }</citation>
100 </citations>
101 </tool>