Mercurial > repos > dereeper > pangenome_explorer
comparison PanExplorer_workflow/PanExplorer.xml @ 1:032f6b3806a3 draft
Uploaded
author | dereeper |
---|---|
date | Thu, 30 May 2024 11:16:08 +0000 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
0:3cbb01081cde | 1:032f6b3806a3 |
---|---|
1 <tool id="PanExplorer2" name="PanExplorer2" version="2.0"> | |
2 <description> Bacterial pan-genome analysis </description> | |
3 <requirements> | |
4 <!-- | |
5 <requirement type="package" version="2.2.26">blast-legacy</requirement> | |
6 <requirement type="package" version="7.480">mafft</requirement> | |
7 <requirement type="package" version="14.137">mcl</requirement> | |
8 <requirement type="package" version="3.697">phylip</requirement> | |
9 <requirement type="package" version="1.7.2">perl-bioperl</requirement> | |
10 --> | |
11 <requirement type="package" version="2.4.2">singularity</requirement> | |
12 | |
13 </requirements> | |
14 | |
15 <command><![CDATA[ | |
16 module load system/singularity/3.6.0; | |
17 export PANEX_PATH=${__tool_directory__}; | |
18 | |
19 perl ${__tool_directory__}/Perl/generateConfig.pl '$private_genomes' '$input' config.yaml '$private_genomes_fasta'; | |
20 | |
21 cat config.yaml >$logfile; | |
22 | |
23 | |
24 | |
25 #if $software=="pgap" | |
26 sed "s/identity\=80/identity\=$min_identity/g" \$PANEX_PATH/Snakemake_files/Snakefile_wget_PGAP_heatmap_upset_COG >snakefile; | |
27 #else if $software=="roary" | |
28 sed "s/identity\=80/identity\=$min_identity/g" \$PANEX_PATH/Snakemake_files/Snakefile_wget_roary_heatmap_upset_COG >snakefile; | |
29 #else if $software=="orthofinder" | |
30 sed "s/identity\=80/identity\=$min_identity/g" \$PANEX_PATH/Snakemake_files/Snakefile_orthofinder_heatmap_upset >snakefile; | |
31 #else if $software=="cactus" | |
32 cp -rf \$PANEX_PATH/Snakemake_files/Snakefile_wget_cactus_heatmap_upset_COG2 snakefile; | |
33 #else if $software=="pggb" | |
34 sed "s/identity\=30/identity\=$min_identity/g" \$PANEX_PATH/Snakemake_files/Snakefile_wget_pggb_heatmap_upset_COG >snakefile; | |
35 #else | |
36 sed "s/identity\=80/identity\=$min_identity/g" \$PANEX_PATH/Snakemake_files/Snakefile_wget_panacota_heatmap_upset_COG >snakefile; | |
37 #end if | |
38 | |
39 singularity exec \$PANEX_PATH/panexplorer.sif snakemake --cores 1 -s snakefile >>$logfile 2>&1; | |
40 | |
41 | |
42 cp -rf outputs/upsetr.svg $upset; | |
43 cp -rf outputs/pav_matrix.tsv $output; | |
44 cp -rf outputs/heatmap.svg.gz $heatmap; | |
45 cp -rf outputs/heatmap.svg.heatmap_plotly.html $heatmap_html; | |
46 cp -rf outputs/rarefaction_curves.txt $rarefaction_curves; | |
47 cp -rf outputs/rarefaction_curves.svg $rarefaction_curves_svg; | |
48 cp -rf outputs/heaps.tsv $heaps; | |
49 cp -rf outputs/heatmap.svg.complete.pdf.distance_matrix.txt $distance_matrix; | |
50 cp -rf outputs/heatmap.svg.complete.pdf.distance_matrix.hclust.newick $njtree; | |
51 | |
52 #if $software!="orthofinder" | |
53 cp -rf outputs/GCskew.txt $gcfile; | |
54 cp -rf outputs/cog_output.txt $cogfile; | |
55 cp -rf outputs/cog_stats.txt $outcog_stat; | |
56 cp -rf outputs/cog_stats2.txt $outcog_stat2; | |
57 cp -rf outputs/cog_of_clusters.txt $outcog_clusters; | |
58 cp -rf outputs/genomes/genes.txt $genes; | |
59 cp -rf outputs/fastani.out.matrix.complete $fastani; | |
60 cp -rf outputs/fastani.out.svg $ani_svg; | |
61 #end if | |
62 | |
63 #if $software=="pggb" | |
64 cp -rf outputs/all_genomes.vcf $vcf; | |
65 #else | |
66 touch $vcf; | |
67 #end if | |
68 ]]></command> | |
69 | |
70 | |
71 <inputs> | |
72 <param name="input" type="text" multiple="true" label="List of genbank identifiers" help="Coma separated list (ex: CP000235.1,CP001079.1,CP001759.1,CP015994.2)"/> | |
73 <param name="private_genomes" type="data" format="zip" label="Zip of genbank or GFF files" optional="true"/> | |
74 <param name="private_genomes_fasta" type="data" format="zip" label="Zip of Fasta files" optional="true"/> | |
75 <param name="min_identity" type="text" value="80" label="Minimum percentage identity for BlastP" /> | |
76 <param type="select" name="software" label="Choose the pan-genome software"> | |
77 <option value="pgap">PGAP</option> | |
78 <option value="roary">Roary</option> | |
79 <option value="panacota">PanACoTA</option> | |
80 <option value="orthofinder">OrthoFinder</option> | |
81 <option value="cactus">Minigraph-Cactus</option> | |
82 <option value="pggb">PanGenome Graph Builder (PGGB)</option> | |
83 </param> | |
84 </inputs> | |
85 | |
86 <outputs> | |
87 <data format="txt" name="output" label="Pangenome presence absence matrix"/> | |
88 <data format="newick" name="njtree" label="PanBased NJ tree"/> | |
89 <data format="txt" name="genes" label="Genes"/> | |
90 <data format="txt" name="cogfile" label="COG assignation"/> | |
91 <data format="txt" name="gcfile" label="GC_percent"/> | |
92 <data format="svg" name="upset" label="Upset Diagram"/> | |
93 <data format="svg" name="heatmap" label="Presence Absence Heatmap"/> | |
94 <data format="html" name="heatmap_html" label="Presence Absence Heatmap interactive"/> | |
95 <data format="tabular" name="outcog_stat" label="COG category counts"/> | |
96 <data format="tabular" name="outcog_stat2" label="COG category 2 counts"/> | |
97 <data format="tabular" name="outcog_clusters" label="COG of clusters"/> | |
98 <data format="tabular" name="fastani" label="ANI" /> | |
99 <data format="svg" name="ani_svg" label="ANI heatmap" /> | |
100 <data format="txt" name="rarefaction_curves" label="Rarefaction curves data"/> | |
101 <data format="svg" name="rarefaction_curves_svg" label="Rarefaction curves"/> | |
102 <data format="txt" name="heaps" label="Heaps law alpha"/> | |
103 <data format="txt" name="distance_matrix" label="Accessory based distance matrix"/> | |
104 <data format="vcf" name="vcf" label="VCF file"/> | |
105 <data format="txt" name="logfile" label="Logfile"/> | |
106 <data format="txt" name="roary_log" label="Roary Logfile"/> | |
107 </outputs> | |
108 | |
109 </tool> |