Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Sam Nooij
Jovian screener
Commits
84041347
Commit
84041347
authored
Jun 24, 2020
by
Sam Nooij
Browse files
Remove unnecessary whitespaces
parent
9773f70d
Changes
1
Hide whitespace changes
Inline
Side-by-side
Snakefile
View file @
84041347
"""
"""
Authors: Sam Nooij
Authors: Sam Nooij
Organisation: Leiden University Medical Centre (LUMC)
Organisation: Leiden University Medical Centre (LUMC)
Department: Medical Microbiology (MM) Research,
Department: Medical Microbiology (MM) Research,
Center for Microbiome Analyses and Therapeutics (CMAT)
Center for Microbiome Analyses and Therapeutics (CMAT)
Date: 2020-03-20
Date: 2020-03-20
...
@@ -10,7 +10,7 @@ Specifically for screening for genes of interest and their species
...
@@ -10,7 +10,7 @@ Specifically for screening for genes of interest and their species
of origin.
of origin.
Input: Jovian's output (classified scaffolds and trimmed reads)
Input: Jovian's output (classified scaffolds and trimmed reads)
Output: report on presence and abundance of species and
Output: report on presence and abundance of species and
genes of interest
genes of interest
Example use:
Example use:
...
@@ -18,7 +18,7 @@ Example use:
...
@@ -18,7 +18,7 @@ Example use:
Or on the Shark cluster:
Or on the Shark cluster:
$ snakemake --use-conda --latency-wait 60 --jobs 16 \
$ snakemake --use-conda --latency-wait 60 --jobs 16 \
--cluster "qsub -q all.q -pe BWA {threads} -l h_vmem={cluster.vmem} -cwd -j Y -o log/cluster/ -V" \
--cluster "qsub -q all.q -pe BWA {threads} -l h_vmem={cluster.vmem} -cwd -j Y -o log/cluster/ -V" \
--cluster-config config/cluster_config.yaml
--cluster-config config/cluster_config.yaml
Or with SLURM:
Or with SLURM:
$ snakemake -p --use-conda --latency-wait 60 --jobs 10 \
$ snakemake -p --use-conda --latency-wait 60 --jobs 10 \
--cluster "sbatch --parsable -J Snakejob-{name}.{jobid} -N 1 -n {threads} --mem={cluster.vmem} -t {cluster.time} -D . -e log/SLURM/{name}-{jobid}.err -o log/SLURM/{name}-{jobid}.out" \
--cluster "sbatch --parsable -J Snakejob-{name}.{jobid} -N 1 -n {threads} --mem={cluster.vmem} -t {cluster.time} -D . -e log/SLURM/{name}-{jobid}.err -o log/SLURM/{name}-{jobid}.out" \
...
@@ -74,7 +74,7 @@ else:
...
@@ -74,7 +74,7 @@ else:
rule all:
rule all:
input:
input:
expand("data/tmp/all_taxClassified-{species}.tsv",
expand("data/tmp/all_taxClassified-{species}.tsv",
species = SPECIES),
species = SPECIES),
#Filtered table for species of interest
#Filtered table for species of interest
...
@@ -92,7 +92,7 @@ rule all:
...
@@ -92,7 +92,7 @@ rule all:
species = SPECIES, ext = ["png", "pdf"]),
species = SPECIES, ext = ["png", "pdf"]),
#pyANI Average Nucleotide Identity plots for each species
#pyANI Average Nucleotide Identity plots for each species
expand(REFERENCE_DIR + "{ref}.blastdb.{ext}",
expand(REFERENCE_DIR + "{ref}.blastdb.{ext}",
ref = REFERENCE_NAMES, ext = ["nhr", "nin", "nsq"]),
ref = REFERENCE_NAMES, ext = ["nhr", "nin", "nsq"]),
#Blast database files for each reference
#Blast database files for each reference
...
@@ -138,7 +138,7 @@ rule all:
...
@@ -138,7 +138,7 @@ rule all:
"data/processed/Depth_of_coverage.tsv",
"data/processed/Depth_of_coverage.tsv",
#Concatenated tables with depth of coverage
#Concatenated tables with depth of coverage
### Step 3: include processing steps to generate output ###
### Step 3: include processing steps to generate output ###
## First, filter out species of interest ##
## First, filter out species of interest ##
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment