view RPKM_count.xml @ 7:bc65a465caed default tip

Merge heads at 5:09a576e512cc and 6:e994caf2503d which were created as a result of a recently fixed bug.

<tool id="RPKM_count" name="RPKM Count">
	<description>calculates raw count and RPKM values for transcript at exon, intron, and mRNA level</description>
	<requirements>
		<requirement type="package" version="0.1.18">samtools</requirement>
		<requirement type="package" version="2.3.7">rseqc</requirement>
	</requirements>
	<command interpreter="python"> samtoolshelper.py RPKM_count.py -i $input -o output -r $refgene

		#if $nx
			-x
		#end if
		
		#if str($strand_type.strand_specific) == "pair"
			-d
			#if str($strand_type.pair_type) == "sd"
				'1++,1--,2+-,2-+'
			#else
				'1+-,1-+,2++,2--'
			#end if
		#end if

		#if str($strand_type.strand_specific) == "single"
			-d
			#if str($strand_type.single_type) == "s"
				'++,--'
			#else
				'+-,-+'
			#end if
		#end if

		#if $skiphits
			-u
		#end if

		#if $onlyexonic
			-e
		#end if

	</command>
	<inputs>
		<param name="input" type="data" format="bam" label="input bam/sam file" />
		<param name="refgene" type="data" format="bed" label="Reference gene model" />
		<conditional name="strand_type">
			<param name="strand_specific" type="select" label="Strand-specific?" value="None">
				<option value="none">None</option>
				<option value="pair">Pair-End RNA-seq</option>
				<option value="single">Single-End RNA-seq</option>
			</param>
			<when value="pair">
				<param name="pair_type" type="select" display="radio" label="Pair-End Read Type (format: mapped --> parent)" value="sd">
					<option value="sd"> read1 (positive --> positive; negative --> negative), read2 (positive --> negative; negative --> positive)</option>
					<option value="ds">read1 (positive --> negative; negative --> positive), read2 (positive --> positive; negative --> negative)</option>
				</param>
			</when>
			<when value="single">
				<param name="single_type" type="select" display="radio" label="Single-End Read Type (format: mapped --> parent)" value="s">
					<option value="s">positive --> positive; negative --> negative</option>
					<option value="d">positive --> negative; negative --> positive</option>
				</param>
			</when>
			<when value="none"></when>
		</conditional>
		<param name="skiphits" type="boolean" value="false" label="Skip Multiple Hit Reads" />
		<param name="onlyexonic" type="boolean" value="false" label="Only use exonic (UTR exons and CDS exons) reads, otherwise use all reads" />
	</inputs>
	<outputs>
		<data format="xls" name="outputxls" from_work_dir="output_read_count.xls"/>
	</outputs>
	<tests>
		<test>
			<param name="input" value="Pairend_nonStrandSpecific_36mer_Human_hg19.bam" />
			<param name="refgene" value="hg19_RefSeq.bed" />
			<output name="outputxls" file="rpkmctout_read_count.xls" />
		</test>
	</tests>
	<help>
.. image:: https://code.google.com/p/rseqc/logo?cct=1336721062

-----

About RSeQC
+++++++++++

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.

The RSeQC package is licensed under the GNU GPL v3 license.

Inputs
++++++++++++++

Input BAM/SAM file
	Alignment file in BAM/SAM format.

Reference gene model
	Gene model in BED format.

Strand sequencing type (default=none)
	See Infer Experiment tool if uncertain.

Options
++++++++++++++

Skip Multiple Hit Reads
	Use Multiple hit reads or use only uniquely mapped reads.

Only use exonic reads 
	Renders program only used exonic (UTR exons and CDS exons) reads, otherwise use all reads.
	
	</help>
</tool>