Skip to content

Commit a769ca8

Browse files
authored
Merge pull request #72 from microbiomedata/tahoma
Update Binning Workflows to run on Tahoma
2 parents 48cf6ce + 91c992e commit a769ca8

3 files changed

Lines changed: 68 additions & 66 deletions

File tree

input.json

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
{
2-
"nmdc_mags.proj_name": "nmdc_wfmgan-xx-xxxxxxxx",
3-
"nmdc_mags.contig_file": "/path/to/Assembly/nmdc_wfmgan-xx-xxxxxxx_contigs.fna",
2+
"nmdc_mags.proj": "nmdc_wfmgan-xx-xxxxxxxx",
3+
"nmdc_mags.contig_file": "/path/to/Assembly/nmdc_wfmgas-xx-xxxxxxx_contigs.fna",
44
"nmdc_mags.sam_file": "/path/to/Assembly/nmdc_wfmgas-xx-xxxxxxx_pairedMapped_sorted.bam",
55
"nmdc_mags.gff_file": "/path/to/Annotation/nmdc_wfmgan-xx-xxxxxxx_functional_annotation.gff",
66
"nmdc_mags.proteins_file": "/path/to/Annotation/nmdc_wfmgan-xx-xxxxxxx_proteins.faa",

mbin_nmdc.wdl

Lines changed: 65 additions & 63 deletions
Original file line numberDiff line numberDiff line change
@@ -2,20 +2,20 @@ version 1.0
22
workflow nmdc_mags {
33
input {
44
String proj
5-
String contig_file
6-
String sam_file
7-
String gff_file
8-
String proteins_file
9-
String cog_file
10-
String ec_file
11-
String ko_file
12-
String pfam_file
13-
String tigrfam_file
14-
String crispr_file
15-
String product_names_file
16-
String gene_phylogeny_file
17-
String lineage_file
18-
String? map_file
5+
File contig_file
6+
File sam_file
7+
File gff_file
8+
File proteins_file
9+
File cog_file
10+
File ec_file
11+
File ko_file
12+
File pfam_file
13+
File tigrfam_file
14+
File crispr_file
15+
File product_names_file
16+
File gene_phylogeny_file
17+
File lineage_file
18+
File? map_file
1919
String? scratch_dir
2020
Int cpu=32
2121
Int threads=64
@@ -68,24 +68,26 @@ workflow nmdc_mags {
6868
mbin_container = container
6969
}
7070
call package {
71-
input: proj = proj,
72-
bins=flatten([mbin_nmdc.hqmq_bin_fasta_files,mbin_nmdc.lq_bin_fasta_files]),
73-
json_stats=mbin_nmdc.stats_json,
74-
gff_file=stage.gff,
75-
proteins_file=stage.proteins,
76-
cog_file=stage.cog,
77-
ec_file=stage.ec,
78-
ko_file=stage.ko,
79-
pfam_file=stage.pfam,
80-
tigrfam_file=stage.tigrfam,
81-
crispr_file=stage.crispr,
82-
gene_phylogeny_file=stage.gene_phylogeny,
83-
product_names_file=stage.product_names,
84-
container=package_container
71+
input:
72+
proj = proj,
73+
bins=flatten([mbin_nmdc.hqmq_bin_fasta_files,mbin_nmdc.lq_bin_fasta_files]),
74+
json_stats=mbin_nmdc.stats_json,
75+
gff_file=stage.gff,
76+
proteins_file=stage.proteins,
77+
cog_file=stage.cog,
78+
ec_file=stage.ec,
79+
ko_file=stage.ko,
80+
pfam_file=stage.pfam,
81+
tigrfam_file=stage.tigrfam,
82+
crispr_file=stage.crispr,
83+
gene_phylogeny_file=stage.gene_phylogeny,
84+
product_names_file=stage.product_names,
85+
container=package_container
8586
}
8687
8788
call finish_mags {
88-
input: container="microbiomedata/workflowmeta:1.1.1",
89+
input:
90+
container="microbiomedata/workflowmeta:1.1.1",
8991
proj=proj,
9092
bacsum= mbin_nmdc.bacsum,
9193
arcsum = mbin_nmdc.arcsum,
@@ -127,18 +129,18 @@ workflow nmdc_mags {
127129

128130
task mbin_nmdc {
129131
input{
130-
File fna
131-
File aln
132-
File gff
133-
File lineage
134132
String name
135-
File? map_file
136-
Int? threads
137-
Int? pthreads
133+
File fna
134+
File aln
135+
File gff
136+
File lineage
137+
File? map_file
138+
Int? threads
139+
Int? pthreads
138140
String gtdbtk_env
139141
String checkm_env
140142
String? eukcc2_env
141-
String mbin_container
143+
String mbin_container
142144
}
143145

144146
command<<<
@@ -214,20 +216,20 @@ task mbin_nmdc {
214216
task stage {
215217
input{
216218
String container
217-
String contig_file
218-
String sam_file
219-
String gff_file
220-
String proteins_file
221-
String cog_file
222-
String ec_file
223-
String ko_file
224-
String pfam_file
225-
String tigrfam_file
226-
String crispr_file
227-
String product_names_file
228-
String gene_phylogeny_file
229-
String lineage_file
230-
String? map_file
219+
File contig_file
220+
File sam_file
221+
File gff_file
222+
File proteins_file
223+
File cog_file
224+
File ec_file
225+
File ko_file
226+
File pfam_file
227+
File tigrfam_file
228+
File crispr_file
229+
File product_names_file
230+
File gene_phylogeny_file
231+
File lineage_file
232+
File? map_file
231233
String contigs_out="contigs.fasta"
232234
String bam_out="pairedMapped_sorted.bam"
233235
String gff_out="functional_annotation.gff"
@@ -355,8 +357,8 @@ task stage {
355357
task check_id_map{
356358
input{
357359
String container
358-
File contig_file
359-
File proteins_file
360+
File contig_file
361+
File proteins_file
360362
String contig_file_name=basename(contig_file)
361363
}
362364
command<<<
@@ -450,22 +452,22 @@ task package{
450452
task finish_mags {
451453
input{
452454
String container
453-
File mbin_sdb
454-
File mbin_version
455+
File mbin_sdb
456+
File mbin_version
455457
String proj
456458
String prefix=sub(proj, ":", "_")
457-
File bacsum
458-
File arcsum
459-
File? short
460-
File? low
461-
File? unbinned
462-
File? checkm
459+
File bacsum
460+
File arcsum
461+
File? short
462+
File? low
463+
File? unbinned
464+
File? checkm
463465
Array[File] hqmq_bin_tarfiles
464466
Array[File] lq_bin_tarfiles
465467
File stats_json
466468
File stats_tsv
467-
Int n_hqmq=length(hqmq_bin_tarfiles)
468-
Int n_lq=length(lq_bin_tarfiles)
469+
Int n_hqmq=length(hqmq_bin_tarfiles)
470+
Int n_lq=length(lq_bin_tarfiles)
469471
File barplot
470472
File heatmap
471473
File kronaplot

version.txt

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1 +1 @@
1-
v1.3.14
1+
v1.3.15

0 commit comments

Comments
 (0)