annotate dante_gff_to_dna.xml @ 10:d0431a839606 draft

Uploaded
author petr-novak
date Wed, 14 Aug 2019 11:24:15 -0400
parents ed4d9ede9cb4
children 039c45c01b47
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>
9
ed4d9ede9cb4 Uploaded
petr-novak
parents: 0
diff changeset
3 <requirements>
ed4d9ede9cb4 Uploaded
petr-novak
parents: 0
diff changeset
4 <requirement type="package">biopython</requirement>
ed4d9ede9cb4 Uploaded
petr-novak
parents: 0
diff changeset
5 </requirements>
0
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
6 <command>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
7 TEMP_DIR_LINEAGES=\$(mktemp -d) &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
8 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
9
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
10 #if $extend_edges:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
11 --extended True
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
12 #else:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
13 --extended False
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
14 #end if
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
15 --classification ${__tool_data_path__ }/protein_domains/${db_type}_class
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
16 &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
17
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
18 cat \$TEMP_DIR_LINEAGES/domains_counts.txt \$TEMP_DIR_LINEAGES/*fasta > $out_fasta &amp;&amp;
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
19 rm -rf \$TEMP_DIR_LINEAGES
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
20 </command>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
21 <inputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
22 <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
23 <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
24 <param name="db_type" type="select" label="Select taxon and protein domain database version (REXdb)" help="">
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
25 <options from_file="rexdb_versions.txt">
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
26 <column name="name" index="0"/>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
27 <column name="value" index="1"/>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
28 </options>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
29 </param>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
30
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
31 <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
32 </inputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
33 <outputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
34 <data format="fasta" name="out_fasta" label="Concatenated fasta domains NT sequences from ${input_dna.hid}" />
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
35 </outputs>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
36
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
37 <help>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
38
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
39 **WHAT IT DOES**
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
40
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
41 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
42 In this case, prior running this tool use DANTE on input DNA:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
43
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
44 1. Protein Domains Finder
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
45 2. Protein Domains Filter (quality filter + type of domain, e.g. RT)
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
46
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
47 INPUTS:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
48 * original DNA sequence in multifasta format to extract the domains from
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
49 * DANTE's output GFF3 file (preferably filtered for quality and specific domain type)
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
50
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
51 OUTPUT:
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
52
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
53 * concatenated fasta file of nucleotide sequences for individual transposons lineages
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
54
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
55 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
56 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
57
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
58 * 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
59 * The extracted sequences are not reported in the case they contain any Ns within the extracted region
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
60
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
61 </help>
77d9f2ecb28a Uploaded
petr-novak
parents:
diff changeset
62 </tool>