view cd_hit_est.xml @ 1:93f25d52cfa5 draft default tip

planemo upload for repository https://github.com/ASaiM/galaxytools/tree/master/tools/cdhit/ commit 4a97d14fadd86dda8abd088023c7f9f2dcab2a88

<tool id="cd_hit_est" name="CD-HIT-EST" version="1.3">
  <description>Cluster a nucleotide dataset into representative sequences</description>

  <macros>
    <import>cdhit_macros.xml</import>
  </macros>

  <expand macro="requirements"/>

  <command><![CDATA[
    cd-hit-est
      -i "$fasta_in" 
      -o rep_seq 
      -c $similarity 
      -n $wordsize $strand
   
      #include source=$common_cdhit_options#
  ]]></command>
  
  <inputs>
    <param name="fasta_in" type="data" format="fasta" label="EST Sequences to cluster" help="(-i)"/>

    <param name="similarity" type="float" value="0.9"  label="Similarity threshold" min=".75" max="1.0" help="(-c)"/>

    <param name="wordsize" type="integer" value="8" label="Word size" min="4" max="10" help="It is suggested to adjust word size in function of similarity threshold. 8,9 or 10  for threshold in [0.9;1.0] interval, 7 for [0.88;0.9], 6 for [0.85;0.88], 5 for [0.80;0.85], 4 for [0.75;0.8] (-n)"/>

    <param name="strand" type="boolean" truevalue="-r 1" falsevalue="" checked="false" label="Compare both strands?"/>

    <expand macro="common_cdhit_options" />
  </inputs>

  <outputs>
    <data format="txt" name="clusters_out" label="${tool.name} on ${on_string}: Clusters" from_work_dir="rep_seq.clstr"/>

    <data format="fasta" name="fasta_out" label="${tool.name} on ${on_string}: Representative sequences" from_work_dir="rep_seq"/>
  </outputs>
  
  <tests>
      <test>
        <param name="fasta_in" value="cd_hit_est_in.fa" />
        <param name="similarity" value="0.9"/>
        <param name="wordsize" value="8"/>
        <param name="strand" value="false"/>
        <param name="settings" value="no"/>
        <param name="tuning" value="default"/>
        <output name="clusters_out" file="est_clusters_output.txt"/>
        <output name="fasta_out" file="est_fasta_output.fasta"/>
      </test>
  </tests>

  <help><![CDATA[

**What it does**

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.

.. _CD-HIT: http://www.bioinformatics.org/cd-hit/

------

**Inputs**

cd-hit-est requires a fasta file as input. 

------

**Outputs**

The first output is a fasta file containing representative sequences.

The second output is a text file listing the mapping of sequences to the representative sequences::

	>Cluster 0
	0 2799aa, >PF04998.6|RPOC2_CHLRE/275-3073... *
	>Cluster 1
	0 2214aa, >PF06317.1|Q6Y625_9VIRU/1-2214... at 80%
	1 2215aa, >PF06317.1|O09705_9VIRU/1-2215... at 84%
	2 2217aa, >PF06317.1|Q6Y630_9VIRU/1-2217... *
	3 2216aa, >PF06317.1|Q6GWS6_9VIRU/1-2216... at 84%
	4 527aa, >PF06317.1|Q67E14_9VIRU/6-532... at 63%
	>Cluster 2
	0 2202aa, >PF06317.1|Q6UY61_9VIRU/8-2209... at 60%
	1 2208aa, >PF06317.1|Q6IVU4_JUNIN/1-2208... *
	2 2207aa, >PF06317.1|Q6IVU0_MACHU/1-2207... at 73%
	3 2208aa, >PF06317.1|RRPO_TACV/1-2208... at 69%
  ]]></help>

  <citations>
    <citation type="doi">10.1093/bioinformatics/btl158</citation>
    <citation type="doi">10.1093/bioinformatics/bts565</citation>
  </citations>
</tool>