Skip to content

Commit 4d73e2d

Browse files
committed
Added assembly paramter in the config file
1 parent 6bc3e1b commit 4d73e2d

File tree

2 files changed

+3
-2
lines changed

2 files changed

+3
-2
lines changed

main.nf

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -266,7 +266,7 @@ process variant_effect_predictor {
266266
INFILE="$infile"
267267
OUTFILE="\${INFILE%.vcf}.vep.vcf"
268268
VEP_CACHE="/sw/data/uppnex/vep/84"
269-
ASSEMBLY="GRCh37"
269+
ASSEMBLY="$params.assembly"
270270
271271
case "\$INFILE" in
272272
*vcf) FORMAT="vcf" ;;
@@ -342,7 +342,7 @@ process snpEff {
342342
sed 's/ID=AD,Number=./ID=AD,Number=R/' "\$INFILE" \
343343
| vt decompose -s - \
344344
| vt normalize -r $params.ref_fasta - \
345-
| java -Xmx7G -jar "\$SNPEFFJAR" -formatEff -classic GRCh37.75 \
345+
| java -Xmx7G -jar "\$SNPEFFJAR" -formatEff -classic ${params.assembly}.75 \
346346
> "\$OUTFILE"
347347
"""
348348
}

nextflow.config

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -5,6 +5,7 @@ params {
55
prefix = ''
66

77
ref_fasta = "/sw/data/uppnex/ToolBox/ReferenceAssemblies/hg38make/bundle/2.8/b37/human_g1k_v37.fasta"
8+
assembly = 'GRCh37'
89

910
// Modules and their versions on the HPC-system
1011
modules {

0 commit comments

Comments
 (0)