|
1
|
1 <tool id="Codon_analysis" name="Codon_density">
|
|
|
2 <description> Analyse Ribo-seq alignments between two conditions to extract codon occupancy </description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package">samtools</requirement>
|
|
|
5 <requirement type="python-module">matplotlib</requirement>
|
|
|
6 <requirement type="python-module">numpy</requirement>
|
|
|
7 <requirement type="python-module">scipy</requirement>
|
|
|
8 <requirement type="python-module">Bio</requirement>
|
|
|
9 </requirements>
|
|
|
10 <command interpreter="python">
|
|
|
11 get_codon_frequency.py
|
|
|
12 --gff=$annot
|
|
|
13 --rep=$replicat_opt['rep']
|
|
|
14 --cond1 $cond1
|
|
|
15 --cond2 $cond2
|
|
|
16 --kmer=$kmer
|
|
|
17 --asite=$asite
|
|
|
18 --out $out
|
|
|
19 --site $site
|
|
|
20 --hex_col1 $color1
|
|
|
21 --hex_col2 $color2
|
|
|
22 --dirout $html_file,$html_file.files_path
|
|
|
23 #if str( $replicat_opt['rep'] ) == 'yes':
|
|
|
24 --file1=$file1,$file11
|
|
|
25 --file2=$file2,$file22
|
|
|
26 #else
|
|
|
27 --file1=$file1
|
|
|
28 --file2=$file2
|
|
|
29 #end if
|
|
|
30
|
|
|
31
|
|
|
32 </command>
|
|
|
33
|
|
|
34 <inputs>
|
|
|
35 <param name="annot" type="data" label="References Input Annotation File (gff)" format="gff" />
|
|
|
36 <conditional name="replicat_opt">
|
|
|
37 <param name="rep" type="select" label="Replicate option">
|
|
|
38 <option value="yes">Yes (only two replicates by condition)</option>
|
|
|
39 <option value="no">No</option>
|
|
|
40 </param>
|
|
|
41 <when value="yes">
|
|
|
42 ## Use conditional balise : if rep =yes : 4 files, else 2 files
|
|
|
43 <param name="file1" type="data" label="First replicate of first condition (bam file)" format="bam" />
|
|
|
44 <param name="file11" type="data" label="Second replicate of first condition (bam file)" format="bam" />
|
|
|
45 <param name="file2" type="data" label="First replicate of second condition (bam file)" format="bam" />
|
|
|
46 <param name="file22" type="data" label="First replicate of second condition (bam file)" format="bam" />
|
|
|
47 </when>
|
|
|
48 <when value="no">
|
|
|
49 <param name="file1" type="data" label="First bam file" format="bam" />
|
|
|
50 <param name="file2" type="data" label="Second bam File" format="bam" />
|
|
|
51 </when>
|
|
|
52 </conditional>
|
|
|
53 <param name="site" type="select" label="Please choose a ribosome site for codon analysis">
|
|
|
54 <option value="A">A</option>
|
|
|
55 <option value="P">P</option>
|
|
|
56 <option value="E">E</option>
|
|
|
57 </param>
|
|
|
58 <param name="asite" type="integer" label="Off-set from the 5'end of the footprint to the A-site" value ="15" />
|
|
|
59 <param name="kmer" type="integer" label="Size of the best phasing reads" value ="28" />
|
|
|
60 <param name="cond1" type="text" size="25" label="Condition one" help="Required even if no replicate" />
|
|
|
61 <param name="cond2" type="text" size="25" label="Condition two" help="Required even if no replicate" />
|
|
|
62 <param name="color1" type="text" size="50" label="Color for first condition" value ="SkyBlue" help="Enter standard name, hex color string, or rgb code. You cand find all colors here : http://pythonhosted.org/ete2/reference/reference_svgcolors.html" />
|
|
|
63 <param name="color2" type="text" size="50" label="Color for second condition" value ="Plum" help="Enter standard name, hex color string, or rgb code. You cand find all colors here : http://pythonhosted.org/ete2/reference/reference_svgcolors.html" />
|
|
|
64
|
|
|
65 </inputs>
|
|
|
66
|
|
|
67 <outputs>
|
|
|
68 <data format="csv" name="out" label="Codon analysis output file on ${on_string}"/>
|
|
|
69 <data format="html" name="html_file" label="Codon analysis results on ${on_string}"/>
|
|
|
70
|
|
|
71 </outputs>
|
|
|
72
|
|
|
73 <help>
|
|
|
74
|
|
|
75 Summary
|
|
|
76 -------
|
|
|
77 This tool uses Ribo-seq (bam file) to compare codon translation between two conditions.
|
|
|
78 For each footprint, codons at choosen site are saved and an histogram with all normalized codon numbers is plotted in both conditions.
|
|
|
79 A second histogram groups all codons corresponding to an amino acid.
|
|
|
80 A chisquare test is used for testing if distribution of used codons is the same in both conditions.
|
|
|
81
|
|
|
82 Output
|
|
|
83 -------
|
|
|
84 This tool provides an html output with graphical outputs and a statistical test result. An additionnal csv file with codon numbers is provided.
|
|
|
85
|
|
|
86
|
|
|
87 Dependances
|
|
|
88 ------------
|
|
|
89
|
|
|
90 .. class:: warningmark
|
|
|
91
|
|
|
92 This tool depends on Python (>=2.7) and following packages : numpy 1.8.0, Biopython 1.58, scipy 0.12.0 and matplotlib 1.3.1. Samtools is used for bam manipulation.
|
|
|
93
|
|
|
94
|
|
|
95
|
|
|
96
|
|
|
97
|
|
|
98
|
|
|
99
|
|
|
100
|
|
|
101
|
|
|
102 </help>
|
|
|
103 </tool>
|
