Skip to content

Commit

Permalink
[FIX] Have FASTQ directory of one sample be input for cellranger_count_8
Browse files Browse the repository at this point in the history
Not the base directory with all samples.
  • Loading branch information
Anne Bertolini committed May 27, 2024
1 parent d932e73 commit 8bb9c85
Showing 1 changed file with 11 additions and 1 deletion.
12 changes: 11 additions & 1 deletion workflow/rules/scAmpi_basic_rules.smk
Original file line number Diff line number Diff line change
Expand Up @@ -48,9 +48,19 @@ rule cellranger_count:
'ln -sr "{params.cr_out}{params.mySample}/outs/metrics_summary.csv" "{params.metrics_summary}" '


# Retrieve the fastqs directory name (ie. uses cellranger sample name) corresponding to a given sample
def get_fastq_dir(wildcards):
"return fastq directory of one sample"
sample = wildcards.sample
fastqs_dir = config["inputOutput"]["input_fastqs"]
sample_fastq_dir = fastqs_dir + sample
return sample_fastq_dir


# Run cellranger v8. Some new parameters are required (e.g. --create-bam)
rule cellranger_count_8:
input:
fastqs_dir=config["inputOutput"]["input_fastqs"],
fastqs_dir=get_fastq_dir,
reference=config["resources"]["reference_transcriptome"],
output:
# the cellranger output cannot be specified directly because this would trigger Snakemake
Expand Down

0 comments on commit 8bb9c85

Please sign in to comment.