Skip to content

Commit

Permalink
Switch to version 1.20 from seqera
Browse files Browse the repository at this point in the history
  • Loading branch information
verku committed Jun 14, 2024
1 parent 801fd7d commit 1e9a163
Show file tree
Hide file tree
Showing 7 changed files with 32 additions and 32 deletions.
2 changes: 1 addition & 1 deletion utilities/mutational_load_snpeff/Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -60,7 +60,7 @@ rule extract_number_of_samples:
log:
os.path.abspath("logs/rule_logs/extract_number_of_samples.log"),
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
threads: 2
shell:
"""
Expand Down
2 changes: 1 addition & 1 deletion workflow/rules/11_ROH.smk
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,7 @@ rule compress_roh_vcf:
log:
"results/logs/11_ROH/{dataset}/" + REF_NAME + ".{dataset}_fmissing{fmiss}.{chr}_compress_roh_vcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools view -Oz -o {output.compressed} {input.vcf} 2> {log} &&
Expand Down
8 changes: 4 additions & 4 deletions workflow/rules/4_genotyping.smk
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ rule variant_calling:
log:
"results/logs/4_genotyping/{dataset}/" + REF_NAME + "/{sample}.{processed}_variant_calling.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools mpileup -Ou -Q 30 -q 30 -B -f {input.ref} {input.bam} | bcftools call -c -M -O b --threads {threads} -o {output.bcf} 2> {log}
Expand All @@ -43,7 +43,7 @@ rule sort_vcfs:
log:
"results/logs/4_genotyping/{dataset}/" + REF_NAME + "/{sample}.{processed}_sort_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools sort -O b -T {params.tmpdir} -o {output.sort} {input.bcf} 2> {log}
Expand All @@ -61,7 +61,7 @@ rule index_sorted_vcfs:
log:
"results/logs/4_genotyping/{dataset}/" + REF_NAME + "/{sample}.{processed}_index_sorted_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools index -o {output.index} {input.sort} 2> {log}
Expand All @@ -80,7 +80,7 @@ rule sorted_vcf_stats:
log:
"results/logs/4_genotyping/{dataset}/" + REF_NAME + "/{sample}.{processed}_sorted_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.sort} > {output.stats} 2> {log}
Expand Down
2 changes: 1 addition & 1 deletion workflow/rules/5_CpG_identification.smk
Original file line number Diff line number Diff line change
Expand Up @@ -76,7 +76,7 @@ rule sorted_bcf2vcf_CpG_id:
log:
"results/logs/5_CpG_identification/{dataset}/" + REF_NAME + "/{sample}.{processed}_sorted_bcf2vcf_CpG_id.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O z -o {output.vcf} {input.bcf} 2> {log}
Expand Down
8 changes: 4 additions & 4 deletions workflow/rules/8.1_vcf_CpG_filtering.smk
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,7 @@ rule sorted_bcf2vcf_CpG_removal:
log:
"results/logs/8.1_vcf_CpG_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_sorted_bcf2vcf_CpG_removal.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O z -o {output.vcf} {input.bcf} 2> {log}
Expand Down Expand Up @@ -124,7 +124,7 @@ rule CpG_vcf2bcf:
log:
"results/logs/8.1_vcf_CpG_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}.no{CpG_method}_CpG_vcf2bcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O b -o {output.bcf} {input.filtered} 2> {log}
Expand All @@ -142,7 +142,7 @@ rule index_CpG_bcf:
log:
"results/logs/8.1_vcf_CpG_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}.no{CpG_method}_index_CpG_bcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools index -o {output.index} {input.bcf} 2> {log}
Expand All @@ -161,7 +161,7 @@ rule CpG_filtered_vcf_stats:
log:
"results/logs/8.1_vcf_CpG_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}.no{CpG_method}_CpG_filtered_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.bcf} > {output.stats} 2> {log}
Expand Down
18 changes: 9 additions & 9 deletions workflow/rules/8.2_vcf_qual_repeat_filtering.smk
Original file line number Diff line number Diff line change
Expand Up @@ -208,7 +208,7 @@ rule remove_snps_near_indels:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_remove_snps_near_indels.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools filter -g 5 -O b --threads {threads} -o {output.snps} {input.bcf} 2> {log}
Expand All @@ -228,7 +228,7 @@ rule filter_vcfs_qual_dp:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_filter_vcfs_qual_dp.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
minDP=`head -n 1 {input.dp} | cut -d' ' -f 2`
Expand Down Expand Up @@ -256,7 +256,7 @@ rule filter_vcfs_allelic_balance:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_filter_vcfs_allelic_balance.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools view -e 'GT="0/1" & (DP4[2]+DP4[3])/(DP4[0]+DP4[1]+DP4[2]+DP4[3]) < 0.2' {input.bcf} | \
Expand All @@ -275,7 +275,7 @@ rule index_filtered_vcfs:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_index_filtered_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools index -o {output.index} {input.bcf} 2> {log}
Expand All @@ -294,7 +294,7 @@ rule filtered_vcf_stats:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_filtered_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.bcf} > {output.stats} 2> {log}
Expand Down Expand Up @@ -348,7 +348,7 @@ rule filtered_bcf2vcf:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_filtered_bcf2vcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O z -o {output.vcf} {input.bcf} 2> {log}
Expand Down Expand Up @@ -384,7 +384,7 @@ rule filtered_vcf2bcf:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_filtered_vcf2bcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O b -o {output.bcf} {input.filtered} 2> {log}
Expand All @@ -402,7 +402,7 @@ rule index_repmasked_vcfs:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_index_repmasked_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools index -o {output.index} {input.bcf} 2> {log}
Expand All @@ -421,7 +421,7 @@ rule repmasked_vcf_stats:
log:
"results/logs/8.2_vcf_qual_repeat_filtering/{dataset}/" + REF_NAME + "/{sample}.{processed}_repmasked_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.bcf} > {output.stats} 2> {log}
Expand Down
24 changes: 12 additions & 12 deletions workflow/rules/9_merge_vcfs.smk
Original file line number Diff line number Diff line change
Expand Up @@ -315,7 +315,7 @@ rule merge_all_vcfs:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + "_merge_all_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
files=`echo {input.bcf} | awk '{{print NF}}'`
Expand All @@ -340,7 +340,7 @@ rule index_merged_vcf:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_index_merged_vcfs.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools index -o {output.index} {input.bcf} 2> {log}
Expand All @@ -359,7 +359,7 @@ rule merged_vcf_stats:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_merged_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.merged} > {output.stats} 2> {log}
Expand Down Expand Up @@ -400,7 +400,7 @@ rule filter_vcf_biallelic:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_filter_vcf_biallelic.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools view -m2 -M2 -v snps -Ob -o {output.bcf} {input.bcf} 2> {log} &&
Expand All @@ -420,7 +420,7 @@ rule biallelic_filtered_vcf_stats:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_biallelic_filtered_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.bcf} > {output.stats} 2> {log}
Expand Down Expand Up @@ -462,7 +462,7 @@ rule filter_vcf_missing:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_fmissing{fmiss}_filter_vcf_missing.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
# only include sites with zero missing data
Expand Down Expand Up @@ -496,7 +496,7 @@ rule remove_chromosomes:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".all_fmissing{fmiss}.autos_remove_chromosomes.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools view {input.bcf} \
Expand Down Expand Up @@ -535,7 +535,7 @@ rule extract_historical_samples:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".historical_fmissing{fmiss}.{chr}_extract_historical_samples.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
params:
samples=hist_sm,
all_samples=ALL_SAMPLES,
Expand Down Expand Up @@ -568,7 +568,7 @@ rule extract_modern_samples:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".modern_fmissing{fmiss}.{chr}_extract_modern_samples.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
params:
samples=mod_sm,
all_samples=ALL_SAMPLES,
Expand Down Expand Up @@ -600,7 +600,7 @@ rule missingness_filtered_vcf_stats:
log:
"results/logs/9_merge_vcfs/" + REF_NAME + ".{dataset}_fmissing{fmiss}.{chr}_missingness_filtered_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.merged} > {output.stats} 2> {log}
Expand Down Expand Up @@ -642,7 +642,7 @@ rule repmasked_bcf2vcf:
log:
"results/logs/9_merge_vcfs/{dataset}/" + REF_NAME + "/{sample}.{processed}_repmasked_bcf2vcf.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools convert -O z -o {output.vcf} {input.bcf} 2> {log}
Expand Down Expand Up @@ -677,7 +677,7 @@ rule biallelic_missing_filtered_vcf_stats:
log:
"results/logs/9_merge_vcfs/{dataset}/" + REF_NAME + "/{sample}.{processed}_fmissing{fmiss}.{chr}_biallelic_missing_filtered_vcf_stats.log",
singularity:
"docker://quay.io/biocontainers/bcftools:1.9--h68d8f2e_9"
"oras://community.wave.seqera.io/library/bcftools:1.20--f18ddc693c3747b4"
shell:
"""
bcftools stats {input.filtered} > {output.stats} 2> {log}
Expand Down

0 comments on commit 1e9a163

Please sign in to comment.