Mercurial > repos > petr-novak > dante

comparison dante_gff_to_dna.xml @ 0:77d9f2ecb28a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author petr-novak
date Wed, 03 Jul 2019 02:45:00 -0400
parents
children ed4d9ede9cb4
comparison
equal deleted inserted replaced
-1:000000000000 0:77d9f2ecb28a
1 <tool id="domains_extract" name="Extract Domains Nucleotide Sequences" version="1.0.0">
2 <description> Tool to extract nucleotide sequences of protein domains found by DANTE </description>
3 <command>
4 TEMP_DIR_LINEAGES=\$(mktemp -d) &amp;&amp;
5 python3 ${__tool_directory__}/dante_gff_to_dna.py --domains_gff ${domains_gff} --input_dna ${input_dna} --out_dir \$TEMP_DIR_LINEAGES
6
7 #if $extend_edges:
8 --extended True
9 #else:
10 --extended False
11 #end if
12 --classification ${__tool_data_path__ }/protein_domains/${db_type}_class
13 &amp;&amp;
14
15 cat \$TEMP_DIR_LINEAGES/domains_counts.txt \$TEMP_DIR_LINEAGES/*fasta > $out_fasta &amp;&amp;
16 rm -rf \$TEMP_DIR_LINEAGES
17 </command>
18 <inputs>
19 <param format="fasta" type="data" name="input_dna" label="Input DNA" help="Choose input DNA sequence(s) to extract the domains from" />
20 <param format="gff" type="data" name="domains_gff" label="Protein domains GFF" help="Choose filtered protein domains GFF3 (DANTE's output)" />
21 <param name="db_type" type="select" label="Select taxon and protein domain database version (REXdb)" help="">
22 <options from_file="rexdb_versions.txt">
23 <column name="name" index="0"/>
24 <column name="value" index="1"/>
25 </options>
26 </param>
27
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"/>
29 </inputs>
30 <outputs>
31 <data format="fasta" name="out_fasta" label="Concatenated fasta domains NT sequences from ${input_dna.hid}" />
32 </outputs>
33
34 <help>
35
36 **WHAT IT DOES**
37
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.
39 In this case, prior running this tool use DANTE on input DNA:
40
41 1. Protein Domains Finder
42 2. Protein Domains Filter (quality filter + type of domain, e.g. RT)
43
44 INPUTS:
45 * original DNA sequence in multifasta format to extract the domains from
46 * DANTE's output GFF3 file (preferably filtered for quality and specific domain type)
47
48 OUTPUT:
49
50 * concatenated fasta file of nucleotide sequences for individual transposons lineages
51
52 By default sequences will be EXTENDED if the alignment reported by LASTAL does not cover the whole protein sequence from the database.
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:
54
55 * The domain cannot be ambiguous, i.e. the FINAL CLASSIFICATION of the domains region corresponds to the last classification level
56 * The extracted sequences are not reported in the case they contain any Ns within the extracted region
57
58 </help>
59 </tool>