annotate geneBody_coverage2.xml @ 30:b5d2f575ccb6

Uploaded
author nilesh
date Thu, 11 Jul 2013 12:33:27 -0400
parents 64dc986da596
children cc5eaa9376d8
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
64dc986da596 Uploaded
nilesh
parents:
diff changeset
1 <tool id="geneBody_coverage" name="Gene Body Converage (Bigwig)">
64dc986da596 Uploaded
nilesh
parents:
diff changeset
2 <description>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
3 Read coverage over gene body.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
4 </description>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
5 <requirements>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
6 <requirement type="package" version="2.15.1">R</requirement>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
7 <requirement type="package" version="2.3.7">rseqc</requirement>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
8 </requirements>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
9 <command interpreter="python">
64dc986da596 Uploaded
nilesh
parents:
diff changeset
10 geneBody_coverage2.py -i $input -r $refgene -o output
64dc986da596 Uploaded
nilesh
parents:
diff changeset
11 </command>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
12 <inputs>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
13 <param name="input" type="data" label="Input bigwig file" format="bigwig" />
64dc986da596 Uploaded
nilesh
parents:
diff changeset
14 <param name="refgene" type="data" label="Reference Genome" format="bed" />
64dc986da596 Uploaded
nilesh
parents:
diff changeset
15 </inputs>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
16 <outputs>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
17 <data name="outputpdf" format="pdf" from_work_dir="output.geneBodyCoverage.pdf" />
64dc986da596 Uploaded
nilesh
parents:
diff changeset
18 <data name="outputr" format="r" from_work_dir="output.geneBodyCoverage_plot.r" />
64dc986da596 Uploaded
nilesh
parents:
diff changeset
19 <data name="outputtxt" format="txt" from_work_dir="output.geneBodyCoverage.txt" />
64dc986da596 Uploaded
nilesh
parents:
diff changeset
20 </outputs>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
21 <help>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
22 .. image:: https://code.google.com/p/rseqc/logo?cct=1336721062
64dc986da596 Uploaded
nilesh
parents:
diff changeset
23
64dc986da596 Uploaded
nilesh
parents:
diff changeset
24 -----
64dc986da596 Uploaded
nilesh
parents:
diff changeset
25
64dc986da596 Uploaded
nilesh
parents:
diff changeset
26 About RSeQC
64dc986da596 Uploaded
nilesh
parents:
diff changeset
27 +++++++++++
64dc986da596 Uploaded
nilesh
parents:
diff changeset
28
64dc986da596 Uploaded
nilesh
parents:
diff changeset
29 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.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
30
64dc986da596 Uploaded
nilesh
parents:
diff changeset
31 The RSeQC package is licensed under the GNU GPL v3 license.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
32
64dc986da596 Uploaded
nilesh
parents:
diff changeset
33 Inputs
64dc986da596 Uploaded
nilesh
parents:
diff changeset
34 ++++++++++++++
64dc986da596 Uploaded
nilesh
parents:
diff changeset
35
64dc986da596 Uploaded
nilesh
parents:
diff changeset
36 Input BAM/SAM file
64dc986da596 Uploaded
nilesh
parents:
diff changeset
37 Alignment file in BAM/SAM format.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
38
64dc986da596 Uploaded
nilesh
parents:
diff changeset
39 Reference gene model
64dc986da596 Uploaded
nilesh
parents:
diff changeset
40 Gene Model in BED format.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
41
64dc986da596 Uploaded
nilesh
parents:
diff changeset
42
64dc986da596 Uploaded
nilesh
parents:
diff changeset
43 Outputs
64dc986da596 Uploaded
nilesh
parents:
diff changeset
44 ++++++++++++++
64dc986da596 Uploaded
nilesh
parents:
diff changeset
45
64dc986da596 Uploaded
nilesh
parents:
diff changeset
46 Read coverage over gene body. This module is used to check if reads coverage is uniform and if there is any 5’/3’ bias. This module scales all transcripts to 100 nt and calculates the number of reads covering each nucleotide position. Finally, it generates a plot illustrating the coverage profile along the gene body. NOTE: this module requires lots of memory for large BAM files, because it load the entire BAM file into memory. We add another script "geneBody_coverage2.py" into v2.3.1 which takes bigwig (instead of BAM) as input. It only use 200M RAM, but users need to convert BAM into WIG, and then WIG into BigWig.
64dc986da596 Uploaded
nilesh
parents:
diff changeset
47
64dc986da596 Uploaded
nilesh
parents:
diff changeset
48 Example output:
64dc986da596 Uploaded
nilesh
parents:
diff changeset
49 .. image:: http://dldcc-web.brc.bcm.edu/lilab/liguow/RSeQC/figure/geneBody_coverage.png
64dc986da596 Uploaded
nilesh
parents:
diff changeset
50
64dc986da596 Uploaded
nilesh
parents:
diff changeset
51
64dc986da596 Uploaded
nilesh
parents:
diff changeset
52
64dc986da596 Uploaded
nilesh
parents:
diff changeset
53 </help>
64dc986da596 Uploaded
nilesh
parents:
diff changeset
54 </tool>