Skip to content
Snippets Groups Projects
symlinks_renamming.smk 2.48 KiB
Newer Older
nicolas.fernandez_ird.fr's avatar
nicolas.fernandez_ird.fr committed
###I###R###D######U###2###3###3#######T###R###A###N###S###V###I###H###M###I####
###                                                                         ###
###    /\  ______      ___ ____ _  _ __   ____ __   ____     ______  /\     ###
###    ||  \ \ \ \    / __( ___( \/ )__\ (  _ (  ) (_  _)   / / / /  ||     ###
###    ||   > > > >  ( (_-.)__) \  /(__)\ )   /)(__ _)(_   < < < <   ||     ###
###    ||  /_/_/_/    \___(____) \(__)(__(_)\_(____(____)   \_\_\_\  ||     ###
###    \/                                                            \/     ###
###                                                                         ###
###I###R###D######U###2###3###3#######T###R###A###N###S###V###I###H###M###I####
# Name ___________________ symlinks_renamming.smk
# Version ________________ v.2025.01
# Author _________________ Nicolas Fernandez
# Affiliation ____________ IRD_U233_TransVIHMI
# Aim ____________________ Create fastq symlinks
# Date ___________________ 2025.01.31
# Latest modifications ___ 2025.03.12
# Use ____________________ snakemake -s Snakefile --use-conda -j
###############################################################################

nicolas.fernandez_ird.fr's avatar
nicolas.fernandez_ird.fr committed
###############################################################################
rule config:
    # Aim: Load configuration file
    # Use: yq -r . <CONFIG_FILE>
    message:
        """
        ~ Configuration ∞ Show analyses settings ~
        """
    conda:
        WORKFLOW
    input:
        config_file = "config/config.yaml"
    output:
        done = temp("done.temp")
    shell:
        "bash workflow/scripts/settings.sh "
        "&& "
        "touch {output.done}"

nicolas.fernandez_ird.fr's avatar
nicolas.fernandez_ird.fr committed
###############################################################################
rule symlinks:
    # Aim: Create fastq files symlinks
    # Use: ln -s <SOURCE> <LINK>
    message:
        """
        ~ Symlinks ∞ Creating symbolic links ~
        Sample: __________ {wildcards.sample}
        Reads: ___________ R{wildcards.mate}
        """
    input:
nicolas.fernandez_ird.fr's avatar
nicolas.fernandez_ird.fr committed
        valid_fastq = lambda wildcards: os.path.abspath(VALID_FASTQ[wildcards.sample][wildcards.mate]),
        done = "done.temp"
nicolas.fernandez_ird.fr's avatar
nicolas.fernandez_ird.fr committed
    output:
        symlink = temp("results/symlinks/{sample}_R{mate}.fastq.gz")
    log:
        "results/10_Reports/tools-log/symlinks/{sample}_R{mate}.log"
    shell:
        "mkdir -p $(dirname {output.symlink}) "
        "&& "
        "ln -sf "
        "{input.valid_fastq} "
        "{output.symlink} "
        "&> {log}"

###############################################################################