annotate geneBody_coverage.xml @ 30:b5d2f575ccb6

Uploaded
author nilesh
date Thu, 11 Jul 2013 12:33:27 -0400
parents ce0995c33291
children cc5eaa9376d8
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
17
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
1 <tool id="geneBody_coverage" name="Gene Body Converage (BAM)">
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
2 <description>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
3 Read coverage over gene body.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
4 </description>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
5 <requirements>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
6 <requirement type="package" version="2.15.1">R</requirement>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
7 <requirement type="package" version="2.3.7">rseqc</requirement>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
8 </requirements>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
9 <command interpreter="python">
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
10 geneBody_coverage.py -i $input -r $refgene -o output
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
11 </command>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
12 <inputs>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
13 <param name="input" type="data" label="Input .bam file" format="bam" />
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
14 <param name="refgene" type="data" label="Reference Genome" format="bed" />
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
15 </inputs>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
16 <outputs>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
17 <data name="outputpdf" format="pdf" from_work_dir="output.geneBodyCoverage.pdf" />
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
18 <data name="outputr" format="r" from_work_dir="output.geneBodyCoverage_plot.r" />
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
19 <data name="outputtxt" format="txt" from_work_dir="output.geneBodyCoverage.txt" />
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
20 </outputs>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
21 <help>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
22 .. image:: https://code.google.com/p/rseqc/logo?cct=1336721062
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
23
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
24 -----
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
25
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
26 About RSeQC
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
27 +++++++++++
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
28
ce0995c33291 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.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
30
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
31 The RSeQC package is licensed under the GNU GPL v3 license.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
32
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
33 Inputs
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
34 ++++++++++++++
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
35
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
36 Input BAM/SAM file
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
37 Alignment file in BAM/SAM format.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
38
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
39 Reference gene model
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
40 Gene Model in BED format.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
41
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
42
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
43 Outputs
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
44 ++++++++++++++
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
45
ce0995c33291 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.
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
47
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
48 Example output:
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
49 .. image:: http://dldcc-web.brc.bcm.edu/lilab/liguow/RSeQC/figure/geneBody_coverage.png
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
50
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
51
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
52
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
53 </help>
ce0995c33291 Uploaded
nilesh
parents:
diff changeset
54 </tool>