annotate RPKM_count.xml @ 29:907d4b021ff6

Uploaded
author nilesh
date Thu, 11 Jul 2013 12:31:33 -0400
parents
children cc5eaa9376d8
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
29
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
1 <tool id="RPKM_count" name="RPKM Count">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
2 <description>calculates raw count and RPKM values for transcript at exon, intron, and mRNA level</description>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
3 <requirements>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
4 <requirement type="package" version="0.1.18">samtools</requirement>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
5 <requirement type="package" version="2.3.7">rseqc</requirement>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
6 </requirements>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
7 <command interpreter="python"> samtoolshelper.py RPKM_count.py -i $input -o output -r $refgene
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
8
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
9 #if $nx
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
10 -x
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
11 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
12
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
13 #if str($strand_type.strand_specific) == "pair"
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
14 -d
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
15 #if str($strand_type.pair_type) == "sd"
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
16 '1++,1--,2+-,2-+'
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
17 #else
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
18 '1+-,1-+,2++,2--'
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
19 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
20 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
21
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
22 #if str($strand_type.strand_specific) == "single"
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
23 -d
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
24 #if str($strand_type.single_type) == "s"
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
25 '++,--'
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
26 #else
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
27 '+-,-+'
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
28 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
29 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
30
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
31 #if $skiphits
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
32 -u
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
33 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
34
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
35 #if $onlyexonic
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
36 -e
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
37 #end if
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
38
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
39 </command>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
40 <inputs>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
41 <param name="input" type="data" format="bam" label="input bam/sam file" />
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
42 <param name="refgene" type="data" format="bed" label="Reference gene model" />
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
43 <conditional name="strand_type">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
44 <param name="strand_specific" type="select" label="Strand-specific?" value="None">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
45 <option value="none">None</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
46 <option value="pair">Pair-End RNA-seq</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
47 <option value="single">Single-End RNA-seq</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
48 </param>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
49 <when value="pair">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
50 <param name="pair_type" type="select" display="radio" label="Pair-End Read Type (format: mapped --> parent)" value="sd">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
51 <option value="sd"> read1 (positive --> positive; negative --> negative), read2 (positive --> negative; negative --> positive)</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
52 <option value="ds">read1 (positive --> negative; negative --> positive), read2 (positive --> positive; negative --> negative)</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
53 </param>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
54 </when>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
55 <when value="single">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
56 <param name="single_type" type="select" display="radio" label="Single-End Read Type (format: mapped --> parent)" value="s">
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
57 <option value="s">positive --> positive; negative --> negative</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
58 <option value="d">positive --> negative; negative --> positive</option>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
59 </param>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
60 </when>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
61 <when value="none"></when>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
62 </conditional>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
63 <param name="skiphits" type="boolean" value="false" label="Skip Multiple Hit Reads" />
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
64 <param name="onlyexonic" type="boolean" value="false" label="Only use exonic (UTR exons and CDS exons) reads, otherwise use all reads" />
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
65 </inputs>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
66 <outputs>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
67 <data format="xls" name="outputxls" from_work_dir="output_read_count.xls"/>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
68 </outputs>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
69 <help>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
70 .. image:: https://code.google.com/p/rseqc/logo?cct=1336721062
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
71
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
72 -----
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
73
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
74 About RSeQC
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
75 +++++++++++
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
76
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
77 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.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
78
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
79 The RSeQC package is licensed under the GNU GPL v3 license.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
80
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
81 Inputs
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
82 ++++++++++++++
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
83
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
84 Input BAM/SAM file
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
85 Alignment file in BAM/SAM format.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
86
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
87 Reference gene model
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
88 Gene model in BED format.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
89
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
90 Strand sequencing type (default=none)
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
91 See Infer Experiment tool if uncertain.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
92
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
93 Options
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
94 ++++++++++++++
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
95
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
96 Skip Multiple Hit Reads
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
97 Use Multiple hit reads or use only uniquely mapped reads.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
98
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
99 Only use exonic reads
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
100 Renders program only used exonic (UTR exons and CDS exons) reads, otherwise use all reads.
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
101
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
102 Sample Output
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
103 ++++++++++++++
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
104
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
105 ===== ===== === ========= ===== =========== ============= ============= ======== ========
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
106 chrom start end accession score gene strand tag count (+) tag count (-) RPKM (+) RPKM (-)
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
107 ===== ===== === ========= ===== =========== ============= ============= ======== ========
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
108 chr1 29213722 29313959 NM_001166007_intron_1 0 + 431 4329 0.086 0.863
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
109 chr1 29314417 29319841 NM_001166007_intron_2 0 + 31 1 0.114 0.004
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
110 chr1 29320054 29323726 NM_001166007_intron_3 0 + 32 0 0.174 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
111 chr1 29323831 29338376 NM_001166007_intron_4 0 + 33 2 0.045 0.003
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
112 chr1 29338419 29342203 NM_001166007_intron_5 0 + 7 0 0.037 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
113 chr1 29342279 29344735 NM_001166007_intron_6 0 + 35 4 0.285 0.033
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
114 chr1 29344954 29356911 NM_001166007_intron_7 0 + 34 2 0.057 0.003
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
115 chr1 29356999 29359604 NM_001166007_intron_8 0 + 19 1 0.146 0.008
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
116 chr1 29359757 29362337 NM_001166007_intron_9 0 + 31 0 0.24 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
117 chr1 29362435 29365765 NM_001166007_intron_10 0 + 11 1 0.066 0.006
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
118 chr1 29365938 29379615 NM_001166007_intron_11 0 + 63 0 0.092 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
119 chr1 29379824 29391493 NM_001166007_intron_12 0 + 383 8 0.656 0.014
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
120 chr1 29391670 29424318 NM_001166007_intron_13 0 + 817 10 0.5 0.006
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
121 chr1 29424447 29435847 NM_001166007_intron_14 0 + 28 0 0.049 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
122 chr1 29435949 29438879 NM_001166007_intron_15 0 + 12 0 0.082 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
123 chr1 29438960 29442210 NM_001166007_intron_16 0 + 22 2 0.135 0.012
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
124 chr1 29442315 29443330 NM_001166007_intron_17 0 + 9 0 0.177 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
125 chr1 29213602 29213722 NM_001166007_exon_1 0 + 164 0 27.321 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
126 chr1 29313959 29314417 NM_001166007_exon_2 0 + 1699 4 74.158 0.175
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
127 chr1 29319841 29320054 NM_001166007_exon_3 0 + 528 1 49.554 0.094
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
128 chr1 29323726 29323831 NM_001166007_exon_4 0 + 168 0 31.985 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
129 chr1 29338376 29338419 NM_001166007_exon_5 0 + 88 0 40.911 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
130 chr1 29342203 29342279 NM_001166007_exon_6 0 + 114 3 29.986 0.789
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
131 chr1 29344735 29344954 NM_001166007_exon_7 0 + 290 10 26.472 0.913
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
132 chr1 29356911 29356999 NM_001166007_exon_8 0 + 146 1 33.166 0.227
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
133 chr1 29359604 29359757 NM_001166007_exon_9 0 + 404 11 52.786 1.437
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
134 chr1 29362337 29362435 NM_001166007_exon_10 0 + 85 7 17.339 1.428
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
135 chr1 29365765 29365938 NM_001166007_exon_11 0 + 198 2 22.88 0.231
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
136 chr1 29379615 29379824 NM_001166007_exon_12 0 + 306 5 29.269 0.478
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
137 chr1 29391493 29391670 NM_001166007_exon_13 0 + 243 7 27.445 0.791
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
138 chr1 29424318 29424447 NM_001166007_exon_14 0 + 298 7 46.18 1.085
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
139 chr1 29435847 29435949 NM_001166007_exon_15 0 + 396 8 77.611 1.568
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
140 chr1 29438879 29438960 NM_001166007_exon_16 0 + 307 0 75.767 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
141 chr1 29442210 29442315 NM_001166007_exon_17 0 + 138 0 26.273 0
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
142 chr1 29443330 29446558 NM_001166007_exon_18 0 + 2434 84 15.074 0.52
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
143 chr1 29213602 29446558 NM_001166007_mRNA 0 + 8006 150 27.704 0.519
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
144 ===== ===== === ========= ===== =========== ============= ============= ======== ========
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
145
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
146 </help>
907d4b021ff6 Uploaded
nilesh
parents:
diff changeset
147 </tool>