annotate bam2wig.xml @ 1:d9ac0d9cb48e

Uploaded
author nilesh
date Mon, 08 Jul 2013 22:03:22 -0400
parents e73f30d518f9
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
1 <tool id="bam2wig" name="BAM to Wiggle">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
2 <description>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
3 converts all types of RNA-seq data from .bam to .wig
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
4 </description>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
5 <requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
6 <requirement type="package" version="2.15.1">R</requirement>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
7 </requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
8 <requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
9 <requirement type="package" version="0.1.18">samtools</requirement>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
10 </requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
11 <requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
12 <requirement type="package" version="2.3.7">rseqc</requirement>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
13 </requirements>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
14 <command interpreter="python">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
15 samtoolshelper.py /home/nilesh/RSeQC-2.3.3/scripts/bam2wig.py -i $input -s $chromsize -o outfile
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
16
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
17 #if str($strand_type.strand_specific) == "pair"
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
18 -d
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
19 #if str($strand_type.pair_type) == "sd"
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
20 '1++,1--,2+-,2-+'
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
21 #else
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
22 '1+-,1-+,2++,2--'
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
23 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
24 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
25
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
26 #if str($strand_type.strand_specific) == "single"
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
27 -d
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
28 #if str($strand_type.single_type) == "s"
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
29 '++,--'
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
30 #else
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
31 '+-,-+'
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
32 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
33 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
34
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
35 #if $wigsum.wigsum_type
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
36 -t $wigsum.totalwig
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
37 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
38
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
39 #if $skipmultihits
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
40 -u
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
41 #end if
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
42 </command>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
43 <inputs>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
44 <param name="input" type="data" label="Input .bam File" format="bam" />
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
45 <param name="chromsize" type="data" label="Chromosome size file (tab or space separated)" format="txt,tabular" />
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
46 <param name="skipmultihits" type="boolean" label="Skip Multiple Hit Reads/Only Use Uniquely Mapped Reads" value="false" />
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
47 <conditional name="wigsum">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
48 <param name="wigsum_type" type="boolean" label="Specify wigsum?" value="false">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
49 </param>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
50 <when value="true">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
51 <param name="totalwig" value="0" type="integer" label="specified wigsum" />
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
52 </when>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
53 <when value="false"></when>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
54 </conditional>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
55 <conditional name="strand_type">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
56 <param name="strand_specific" type="select" label="Strand-specific?" value="None">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
57 <option value="none">None</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
58 <option value="pair">Pair-End RNA-seq</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
59 <option value="single">Single-End RNA-seq</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
60 </param>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
61 <when value="pair">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
62 <param name="pair_type" type="select" display="radio" label="Pair-End Read Type (format: mapped --> parent)" value="sd">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
63 <option value="sd"> read1 (positive --> positive; negative --> negative), read2 (positive --> negative; negative --> positive)</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
64 <option value="ds">read1 (positive --> negative; negative --> positive), read2 (positive --> positive; negative --> negative)</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
65 </param>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
66 </when>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
67 <when value="single">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
68 <param name="single_type" type="select" display="radio" label="Single-End Read Type (format: mapped --> parent)" value="s">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
69 <option value="s">positive --> positive; negative --> negative</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
70 <option value="d">positive --> negative; negative --> positive</option>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
71 </param>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
72 </when>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
73 <when value="none"></when>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
74 </conditional>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
75 </inputs>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
76 <outputs>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
77 <data format="wig" name="output" from_work_dir="outfile.wig">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
78 <filter>strand_type['strand_specific'] == 'none'</filter>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
79 </data>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
80 <data format="wig" name="outputfwd" from_work_dir="outfile_Forward.wig">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
81 <filter>strand_type['strand_specific'] != 'none'</filter>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
82 </data>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
83 <data format="wig" name="outputrv" from_work_dir="outfile_Reverse.wig">
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
84 <filter>strand_type['strand_specific'] != 'none'</filter>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
85 </data>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
86 </outputs>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
87 <help>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
88 .. image:: https://code.google.com/p/rseqc/logo?cct=1336721062
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
89
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
90 -----
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
91
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
92 About RSeQC
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
93 +++++++++++
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
94
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
95 The RSeQC package provides a number of useful modules that can comprehensively evaluate high throughput sequence data especially RNA-seq data. “Basic modules” quickly inspect sequence quality, nucleotide composition bias, PCR bias and GC bias, while “RNA-seq specific modules” investigate sequencing saturation status of both splicing junction detection and expression estimation, mapped reads clipping profile, mapped reads distribution, coverage uniformity over gene body, reproducibility, strand specificity and splice junction annotation.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
96
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
97 The RSeQC package is licensed under the GNU GPL v3 license.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
98
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
99 Inputs
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
100 ++++++++++++++
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
101
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
102 Input BAM file
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
103 Alignment file in BAM format (SAM is not supported). BAM file will be sorted and indexed using samTools.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
104
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
105 Chromosome size file
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
106 Tab or space separated text file with 2 columns: first column is chromosome name, second column is size of the chromosome. Chromosome names (such as "chr1") should be consistent between this file and BAM file.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
107
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
108 Specified wigsum (default=none)
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
109 Specified wigsum. Wigsum of 100000000 equals to coverage achieved by 1 million 100nt reads. Ignore this option to disable normalization.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
110
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
111 Skip multiple Hit reads
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
112 skips multiple hit reads or only use uniquely mapped reads
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
113
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
114 Strand-specific (default=none)
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
115 How read(s) were stranded during sequencing. If you are not sure about the strand rule, run infer_experiment.py
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
116
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
117 Outputs
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
118 ++++++++++++++
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
119
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
120 If RNA-seq is not strand specific, one wig file will be generated, if RNA-seq
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
121 is strand specific, two wig files corresponding to Forward and Reverse will be generated.
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
122
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
123
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
124 </help>
e73f30d518f9 Uploaded
nilesh
parents:
diff changeset
125 </tool>