0
|
1 <tool id="identify_chimeric_seqs" name="identify_chimeric_seqs" version="2.0.0">
|
|
2 <description>Identify chimeric sequences in input FASTA file</description>
|
|
3 <requirements>
|
|
4 <requirement type="binary">identify_chimeric_seqs.py</requirement>
|
|
5 </requirements>
|
|
6 <command interpreter="python">
|
|
7 qiime_wrapper.py
|
|
8 identify_chimeric_seqs.py
|
|
9 --input_fasta_fp=$input_fasta_fp
|
|
10 #if $pick.chimera_detection_method == 'ChimeraSlayer':
|
|
11 --chimera_detection_method=$pick.chimera_detection_method
|
|
12 --aligned_reference_seqs_fp=$pick.aligned_reference_seqs_fp
|
|
13 #if $pick.min_div_ratio.__str__ != '0.0':
|
|
14 --min_div_ratio=$pick.min_div_ratio
|
|
15 #end if
|
|
16 #elif $pick.chimera_detection_method == 'blast_fragments':
|
|
17 --chimera_detection_method=$pick.chimera_detection_method
|
|
18 --id_to_taxonomy_fp=$pick.id_to_taxonomy_fp
|
|
19 #if $pick.blast_db != None and $pick.blast_db.__str__ != "":
|
|
20 --blast_db=$pick.blast_db
|
|
21 #else:
|
|
22 --reference_seqs_fp=$pick.reference_seqs_fp
|
|
23 #end if
|
|
24 --num_fragments=$pick.num_fragments
|
|
25 --taxonomy_depth=$pick.taxonomy_depth
|
|
26 --max_e_value=$pick.max_e_value
|
|
27 #end if
|
|
28 --output_fp=$output_fp
|
|
29 </command>
|
|
30 <inputs>
|
|
31 <param name="input_fasta_fp" type="data" format="fasta" label="input_fasta_fp"
|
|
32 help="path to the input fasta file (usually result from pick_rep_set) [REQUIRED]"/>
|
|
33 <conditional name="pick">
|
|
34 <param name="chimera_detection_method" type="select" label="chimera_detection_method"
|
|
35 help="Chimera detection method. Choices: blast_fragments or ChimeraSlayer. [default:ChimeraSlayer]">
|
|
36 <option value="blast_fragments">blast_fragments</option>
|
|
37 <option value="ChimeraSlayer" selected="true">ChimeraSlayer</option>
|
|
38 </param>
|
|
39 <when value="ChimeraSlayer">
|
|
40 <param name="aligned_reference_seqs_fp" type="data" format="txt" label="aligned_reference_seqs_fp"
|
|
41 help="Path to (Py)Nast aligned reference sequences. [REQUIRED]"/>
|
|
42 <param name="min_div_ratio" type="float" value="0.0" label="min_div_ratio"
|
|
43 help="min divergence ratio (passed to ChimeraSlayer). If set to None uses ChimeraSlayer default value. [default: None]"/>
|
|
44 </when>
|
|
45 <when value="blast_fragments">
|
|
46 <param name="blast_db" type="text" label="blast_db"
|
|
47 help="Database to blast against. Must provide either blast_db or reference_seqs_fp when method is blast_fragments [default: None]"/>
|
|
48 <param name="reference_seqs_fp" type="data" format="txt" label="reference_seqs_fp"
|
|
49 help="Path to reference sequences (used to build a blast db when method blast_fragments). [default: None; REQUIRED if no blast_db is provided;]"/>
|
|
50 <param name="id_to_taxonomy_fp" type="data" format="tabular" label="id_to_taxonomy_fp"
|
|
51 help="Path to tab-delimited file mapping sequences to assigned taxonomy. Each assigned taxonomy is provided as a comma-separated list. [REQUIRED]"/>
|
|
52 <param name="taxonomy_depth" type="integer" value="4" label="taxonomy_depth"
|
|
53 help="Number of taxonomic divisions to consider when comparing taxonomy assignments [default: 4]"/>
|
|
54 <param name="num_fragments" type="integer" value="3" label="num_fragments"
|
|
55 help="Number of fragments to split sequences into (i.e., number of expected breakpoints + 1) [default: 3]"/>
|
|
56 <param name="max_e_value" type="float" value="1e-30" label="max_e_value"
|
|
57 help="Max e-value to assign taxonomy [default: 1e-30]"/>
|
|
58 </when>
|
|
59 </conditional> <!-- pick -->
|
|
60 </inputs>
|
|
61 <outputs>
|
|
62 <data format="txt" name="output_fp" label="${tool.name} on ${on_string}: chimeras"/>
|
|
63 </outputs>
|
|
64 <tests>
|
|
65 </tests>
|
|
66 <help>For more information, see identify_chimeric_seqs_ in the Qiime documentation.
|
|
67
|
|
68 Updated and validated 01/19/12 by Amanda Zuzolo, Microbiome Analysis Center, George Mason University, Fairfax, VA
|
|
69 Qiime integration courtesy Jim Johnson, Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN
|
|
70
|
|
71 .. _identify_chimeric_seqs: http://qiime.org/scripts/identify_chimeric_seqs.html</help>
|
|
72 </tool>
|
|
73
|