|
1
|
1 <tool id="Codon_analysis" name="Codon_density">
|
|
13
|
2 <description>To compare Ribo-seq alignments between two sets of conditions, to determine codon occupancy</description>
|
|
1
|
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>
|
|
11
|
35 <param name="annot" type="data" label="Reference annotation file (GFF)" format="gff" />
|
|
1
|
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
|
|
11
|
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" />
|
|
1
|
47 </when>
|
|
|
48 <when value="no">
|
|
11
|
49 <param name="file1" type="data" label="First Bam file" format="bam" />
|
|
|
50 <param name="file2" type="data" label="Second Bam File" format="bam" />
|
|
1
|
51 </when>
|
|
|
52 </conditional>
|
|
11
|
53 <param name="site" type="select" label="Please choose a ribosomal site for codon analysis">
|
|
1
|
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" />
|
|
11
|
59 <param name="kmer" type="integer" label="Lenght of the best phasing footprints" value ="28" />
|
|
13
|
60 <param name="cond1" type="text" size="25" label="Condition one" help="Required even if no replicate" value="cond1"/>
|
|
|
61 <param name="cond2" type="text" size="25" label="Condition two" help="Required even if no replicate" value="cond2"/>
|
|
1
|
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 -------
|
|
11
|
77 This tool uses Ribo-seq (BAM file) to determine whether codon occupancy differs between two sets of conditions. For each footprint, the codons at the chosen site are recorded and a histogram displaying all the normalised codon numbers is plotted for both sets of conditions. A second histogram groups together all the codons corresponding to a given amino acid. A chi-squared test is then carried out to determine whether the distribution of the codons used is the same in both sets of conditions.
|
|
|
78
|
|
1
|
79
|
|
|
80 Output
|
|
|
81 -------
|
|
11
|
82 This tool provides an html file containing graphs and a statistical test result. An additional csv file with codon numbers is provided.
|
|
1
|
83
|
|
|
84
|
|
|
85 Dependances
|
|
|
86 ------------
|
|
|
87
|
|
|
88 .. class:: warningmark
|
|
|
89
|
|
|
90 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.
|
|
|
91
|
|
|
92
|
|
|
93
|
|
|
94
|
|
|
95
|
|
|
96
|
|
|
97
|
|
|
98
|
|
|
99
|
|
|
100 </help>
|
|
|
101 </tool>
|
