view reactivity_cal/reactivity_calculation.xml~ @ 33:1a92d934f8d1 draft

Uploaded

<tool id="react_cal_pipeline" name="Reactivity calculation" version="1.0">
	<description></description>
	<command interpreter="python">react_cal.py $dist_file1 $dist_file2 $seq_file $nt_spec $flag_in $threshold $output </command>
        <requirements>
                <requirement type="package" version="1.61">biopython</requirement>
        </requirements>
	<inputs>
                <param name="dist_file1" type="data" format="txt" label="RTSC file for (+) library"/>
		<param name="dist_file2" type="data" format="txt" label="RTSC file for (-) library"/>
                <param name="seq_file" type="data" format="fasta" label="Reference library"/>
                <param name="nt_spec" type="text" value = "AC" label="Nucleotide specificity"/>
                <param name="flag_in" type="boolean" checked="true" truevalue = "1" falsevalue = "0" label="Normalization is performed if checked"/>
                <param name="threshold" type="float" value = "7" label="Value to cap the reactivities"/>
	</inputs>
	<outputs>
		<data name="output" format="txt"/>
	</outputs>

	<help>


**TIPS**:

-----

**Input**:

* 1. RTSC files (Output of II) for (+) and (-) library
* 2. Reference file (fasta) used to map the reads
* 3. Nucleotide Specificity (Type of nucleotide to have reactivity, e.g. AC for DMS and ACTG for SHAPE)
* [Optional]:
* 1. A threshold to cap the structural reactivities. {Default: 7}
* 2. Flag that determines whether to perform 2%-8% normalization {Default: yes}

-----

**Output**:

A text file with structural reactivity for each nucleotide (Reactivity file)



	</help>
</tool>