annotate dante_gff_to_dna.xml @ 0:77d9f2ecb28a draft

Uploaded
author petr-novak
date Wed, 03 Jul 2019 02:45:00 -0400
parents
children ed4d9ede9cb4
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
1 <tool id="domains_extract" name="Extract Domains Nucleotide Sequences" version="1.0.0">
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
2 <description> Tool to extract nucleotide sequences of protein domains found by DANTE </description>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
3 <command>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
4 TEMP_DIR_LINEAGES=\$(mktemp -d) &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
5 python3 ${__tool_directory__}/dante_gff_to_dna.py --domains_gff ${domains_gff} --input_dna ${input_dna} --out_dir \$TEMP_DIR_LINEAGES
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
6
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
7 #if $extend_edges:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
8 --extended True
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
9 #else:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
10 --extended False
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
11 #end if
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
12 --classification ${__tool_data_path__ }/protein_domains/${db_type}_class
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
13 &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
14
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
15 cat \$TEMP_DIR_LINEAGES/domains_counts.txt \$TEMP_DIR_LINEAGES/*fasta > $out_fasta &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
16 rm -rf \$TEMP_DIR_LINEAGES
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
17 </command>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
18 <inputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
19 <param format="fasta" type="data" name="input_dna" label="Input DNA" help="Choose input DNA sequence(s) to extract the domains from" />
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
20 <param format="gff" type="data" name="domains_gff" label="Protein domains GFF" help="Choose filtered protein domains GFF3 (DANTE's output)" />
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
21 <param name="db_type" type="select" label="Select taxon and protein domain database version (REXdb)" help="">
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
22 <options from_file="rexdb_versions.txt">
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
23 <column name="name" index="0"/>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
24 <column name="value" index="1"/>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
25 </options>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
26 </param>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
27
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
28 <param name="extend_edges" type="boolean" truevalue="True" falsevalue="False" checked="True" label="Extend sequence edges" help="Extend extracted sequence edges to the full length of database domains sequences"/>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
29 </inputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
30 <outputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
31 <data format="fasta" name="out_fasta" label="Concatenated fasta domains NT sequences from ${input_dna.hid}" />
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
32 </outputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
33
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
34 <help>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
35
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
36 **WHAT IT DOES**
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
37
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
38 This tool extracts nucleotide sequences of protein domains from reference DNA based on DANTE's output. It can be used e.g. for deriving phylogenetic relations of individual mobile elements within a species. This can be done separately for individual protein domains types.
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
39 In this case, prior running this tool use DANTE on input DNA:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
40
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
41 1. Protein Domains Finder
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
42 2. Protein Domains Filter (quality filter + type of domain, e.g. RT)
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
43
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
44 INPUTS:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
45 * original DNA sequence in multifasta format to extract the domains from
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
46 * DANTE's output GFF3 file (preferably filtered for quality and specific domain type)
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
47
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
48 OUTPUT:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
49
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
50 * concatenated fasta file of nucleotide sequences for individual transposons lineages
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
51
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
52 By default sequences will be EXTENDED if the alignment reported by LASTAL does not cover the whole protein sequence from the database.
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
53 As the result, the corresponding nucleotide region of the WHOLE aligned database domain will be reported. For every record in the GFF3 file the sequence is reported for the BEST HIT within the domain region under following conditions:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
54
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
55 * The domain cannot be ambiguous, i.e. the FINAL CLASSIFICATION of the domains region corresponds to the last classification level
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
56 * The extracted sequences are not reported in the case they contain any Ns within the extracted region
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
57
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
58 </help>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
59 </tool>