Skip to content
Snippets Groups Projects
Commit d75cd8be authored by nicolas.fernandez_ird.fr's avatar nicolas.fernandez_ird.fr :shinto_shrine:
Browse files

Fix: l.544 pass 'reference' as input (not params)

parent 131ad259
No related branches found
No related tags found
No related merge requests found
Pipeline #75931 passed
......@@ -540,20 +540,19 @@ rule bedtools_masking:
"""
conda:
BEDTOOLS
params:
path = REF_PATH
input:
reference = "resources/genomes/{reference}.fasta",
low_cov_mask = "results/03_Coverage/{reference}/bed/{sample}_{aligner}_{min_cov}X_low-cov-mask.bed"
output:
masked_ref = "results/04_Variants/{reference}/{sample}_{aligner}_{min_cov}X_masked-ref.fasta"
log:
"results/10_Reports/tools-log/bedtools/{reference}/{sample}_{aligner}_{min_cov}X_masking.log"
shell:
"bedtools maskfasta " # Bedtools maskfasta, mask a fasta file based on feature coordinates
"-fi {params.path}{wildcards.reference}.fasta " # Input FASTA file
"-bed {input.low_cov_mask} " # BED/GFF/VCF file of ranges to mask in -fi
"-fo {output.masked_ref} " # Output masked FASTA file
"&> {log}" # Log redirection
"bedtools maskfasta " # Bedtools maskfasta, mask a fasta file based on feature coordinates
"-fi {input.reference} " # Input FASTA file
"-bed {input.low_cov_mask} " # BED/GFF/VCF file of ranges to mask in -fi
"-fo {output.masked_ref} " # Output masked FASTA file
"&> {log}" # Log redirection
###############################################################################
rule bedtools_merged_mask:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment