|
11
|
1 <tool id="cd_hit_est" name="CD-HIT-EST" version="1.2">
|
|
0
|
2 <description>Cluster a nucleotide dataset into representative sequences</description>
|
|
|
3 <requirements>
|
|
2
|
4 <requirement type="package" version="4.6.1">cd-hit</requirement>
|
|
0
|
5 </requirements>
|
|
11
|
6 <macros>
|
|
|
7 <import>cdhit_macros.xml</import>
|
|
|
8 </macros>
|
|
0
|
9 <command>
|
|
11
|
10 cd-hit-est -i "$fasta_in" -o rep_seq -c $similarity -n $wordsize $strand
|
|
|
11 #include source=$common_cdhit_options#
|
|
|
12 #include source=$runtime_tuning#
|
|
0
|
13 </command>
|
|
|
14 <inputs>
|
|
|
15 <param name="fasta_in" type="data" format="fasta" label="EST Sequences to cluster"/>
|
|
|
16 <param name="similarity" type="float" value="0.9" label="similarity threshold: .75 - 1.0, default is .9">
|
|
|
17 <validator type="in_range" message="sequence similarity threshold should be .75 - 1.0" min=".75" max="1.0"/>
|
|
|
18 </param>
|
|
|
19 <param name="wordsize" type="integer" value="8" label="word size">
|
|
|
20 <help> Suggested word size:
|
|
|
21 8,9,10 for thresholds 0.90 ~ 1.0
|
|
|
22 7 for thresholds 0.88 ~ 0.9
|
|
|
23 6 for thresholds 0.85 ~ 0.88
|
|
|
24 5 for thresholds 0.80 ~ 0.85
|
|
|
25 4 for thresholds 0.75 ~ 0.8
|
|
|
26 </help>
|
|
|
27 <validator type="in_range" message="word size should be between 4 and 10" min="4" max="10"/>
|
|
|
28 </param>
|
|
|
29 <param name="strand" type="boolean" truevalue="-r 1" falsevalue="" checked="false" label="Compare both strands"/>
|
|
11
|
30 <expand macro="common_cdhit_options" />
|
|
|
31 <expand macro="runtime_tuning" />
|
|
0
|
32 </inputs>
|
|
|
33 <outputs>
|
|
|
34 <data format="txt" name="clusters_out" label="${tool.name} on ${on_string}: clusters" from_work_dir="rep_seq.clstr"/>
|
|
|
35 <data format="fasta" name="fasta_out" label="${tool.name} on ${on_string}: representatives.fasta" from_work_dir="rep_seq"/>
|
|
|
36 </outputs>
|
|
|
37 <tests>
|
|
3
|
38 <test>
|
|
11
|
39 <!-- Expect 3 clusters: 0,1,2 -->
|
|
3
|
40 <param name="fasta_in" value="cd_hit_est_in.fa" />
|
|
|
41 <param name="similarity" value="0.9"/>
|
|
|
42 <param name="wordsize" value="8"/>
|
|
11
|
43 <param name="strand" value="true"/>
|
|
|
44 <!-- conditionals in macros -->
|
|
|
45 <param name="settings" value="no"/>
|
|
|
46 <param name="tuning" value="default"/>
|
|
3
|
47 <output name="clusters_out">
|
|
|
48 <assert_contents>
|
|
11
|
49 <has_text text=">Cluster 0" />
|
|
|
50 <!-- There should not be a Cluster 3 -->
|
|
|
51 <not_has_text text="Cluster 3" />
|
|
|
52 <has_text_matching expression="F12Fcsw_481739" />
|
|
|
53 </assert_contents>
|
|
|
54 </output>
|
|
|
55 <output name="fasta_out">
|
|
|
56 <assert_contents>
|
|
|
57 <has_text_matching expression="^>[MF]\d\dFcsw_\d*" />
|
|
|
58 </assert_contents>
|
|
|
59 </output>
|
|
|
60 </test>
|
|
|
61 <test>
|
|
|
62 <!-- tighter constraints should yield more clusters -->
|
|
|
63 <param name="fasta_in" value="cd_hit_est_in.fa" />
|
|
|
64 <param name="similarity" value="0.95"/>
|
|
|
65 <param name="wordsize" value="9"/>
|
|
|
66 <param name="strand" value="true"/>
|
|
|
67 <!-- conditionals in macros -->
|
|
|
68 <param name="settings" value="no"/>
|
|
|
69 <param name="tuning" value="default"/>
|
|
|
70 <output name="clusters_out">
|
|
|
71 <assert_contents>
|
|
|
72 <has_text text=">Cluster 4" />
|
|
3
|
73 <has_text_matching expression=">F12Fcsw_481739" />
|
|
|
74 </assert_contents>
|
|
|
75 </output>
|
|
|
76 <output name="fasta_out">
|
|
|
77 <assert_contents>
|
|
|
78 <has_text_matching expression="^>[MF]\d\dFcsw_\d*" />
|
|
|
79 </assert_contents>
|
|
|
80 </output>
|
|
|
81 </test>
|
|
0
|
82 </tests>
|
|
|
83 <help>
|
|
|
84 **CD-HIT-EST**
|
|
|
85
|
|
|
86 CD-HIT_ stands for Cluster Database at High Identity with Tolerance. The program (cd-hit) takes a fasta format sequence database as input and produces a set of 'non-redundant' (nr) representative sequences as output. In addition cd-hit outputs a cluster file, documenting the sequence 'groupies' for each nr sequence representative. The idea is to reduce the overall size of the database without removing any sequence information by only removing 'redundant' (or highly similar) sequences. This is why the resulting database is called non-redundant (nr). Essentially, cd-hit produces a set of closely related protein families from a given fasta sequence database.
|
|
|
87
|
|
|
88 .. _CD-HIT: http://www.bioinformatics.org/cd-hit/
|
|
|
89
|
|
|
90 ------
|
|
|
91
|
|
|
92 **Inputs**
|
|
|
93
|
|
|
94 cd-hit-est requires a fasta dataset as input.
|
|
|
95
|
|
|
96 ------
|
|
|
97
|
|
|
98 **Outputs**
|
|
|
99
|
|
|
100 A fasta datasets containing representative sequences.
|
|
|
101
|
|
|
102 A text file listing the mapping of sequences to the representative sequences::
|
|
|
103
|
|
|
104 >Cluster 0
|
|
|
105 0 2799aa, >PF04998.6|RPOC2_CHLRE/275-3073... *
|
|
|
106 >Cluster 1
|
|
|
107 0 2214aa, >PF06317.1|Q6Y625_9VIRU/1-2214... at 80%
|
|
|
108 1 2215aa, >PF06317.1|O09705_9VIRU/1-2215... at 84%
|
|
|
109 2 2217aa, >PF06317.1|Q6Y630_9VIRU/1-2217... *
|
|
|
110 3 2216aa, >PF06317.1|Q6GWS6_9VIRU/1-2216... at 84%
|
|
|
111 4 527aa, >PF06317.1|Q67E14_9VIRU/6-532... at 63%
|
|
|
112 >Cluster 2
|
|
|
113 0 2202aa, >PF06317.1|Q6UY61_9VIRU/8-2209... at 60%
|
|
|
114 1 2208aa, >PF06317.1|Q6IVU4_JUNIN/1-2208... *
|
|
|
115 2 2207aa, >PF06317.1|Q6IVU0_MACHU/1-2207... at 73%
|
|
|
116 3 2208aa, >PF06317.1|RRPO_TACV/1-2208... at 69%
|
|
|
117
|
|
|
118
|
|
|
119 </help>
|
|
|
120 </tool>
|
