Mercurial > repos > cpt > cpt_prophage_relatedness

comparison prophage_relatedness.xml @ 0:7a23dda2e932 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 7ebbd0df0aea9e58c4df58b61d6da385ee0ebb49
author cpt
date Thu, 08 Aug 2024 03:09:32 +0000
parents
children 7ba8b1f0fdf0
comparison
equal deleted inserted replaced
-1:000000000000 0:7a23dda2e932
1 <tool id="edu.tamu.cpt2.phage.relatedProphages" name="Related Prophages tool" version="21.1.0.0">
2 <description/>
3 <macros>
4 <import>macros.xml</import>
5 </macros>
6 <requirements>
7 <requirement type="package" version="3.7">python</requirement>
8 <requirement type="package" version="1.77">biopython</requirement>
9 </requirements>
10 <command detect_errors="aggressive"><![CDATA[
11 python $__tool_directory__/prophageRelatedness.py
12 '${blastIn.blast}'
13 --outFile '$output'
14 --padding '$padding'
15 --cutoff '$cutoff'
16 --numReturn '$returns'
17 #if '$blastIn.blastType' == "TSV":
18 --isTSV
19 #end if
20 ]]></command>
21 <inputs>
22 <conditional name="blastIn">
23 <param name="blastType" type="select" label="Blastn Input Type">
24 <option value="XML" selected="true">Blast XML</option>
25 <option value="TSV">Blast 25-Column Tabular</option>
26 </param>
27 <when value="XML">
28 <param label="Blastn Results (Blast XML)" name="blast" type="data" format="blastxml"/>
29 </when>
30 <when value="TSV">
31 <param label="Blastn Results (Blast 25-Column Tabular)" name="blast" type="data" format="tsv, tabular"/>
32 </when>
33 </conditional>
34 <param label="Cluster Window" name="padding" type="text" value="2.0" help="Nucleotide window for HSPs to form a cluster (Multiplicative of length of query sequence)"/>
35 <param label="Score cutoff" name="cutoff" type="text" value=".3"/>
36 <param label="Number of Results to return" name="returns" type="integer" value="20"/>
37 </inputs>
38 <outputs>
39 <data format="tabular" name="output" label="Top related prophages"/>
40 </outputs>
41 <tests/>
42 <help><![CDATA[
43 **What it does**
44
45 Filters BLAST results for high-scoring clusters of HSPs.
46
47 The script first determines a window of X nucleotide bases, where X is the length of the query multiplied by the value supplied in the cluster window input. Then for each hit in the blast record, the high scoring pairs are sorted into sets based on the number of HSPs which would fall within that window.&ast; Finally, the total number of identities in the set is divided by the query length, and this score is used to return the top results. If a set would make it into the final results, any other HSP sets from that same hit will also be returned.
48
49 &ast; Sets are formed with a greedy method, where creating a set with the largest possible number of HSPs is the first action, then the next set is the largest possible number of HSPs that is also disjoint from the first, and so on until all HSPs are in a set.
50
51 ]]></help>
52 <expand macro="citations-2020-AJC-solo"/>
53 </tool>