Skip to content

Commit 52e95ea

Browse files
committed
Avoid hard-codig reference name
1 parent 4a892a2 commit 52e95ea

2 files changed

Lines changed: 3 additions & 1 deletion

File tree

workflow/rules/vaf.smk

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -127,6 +127,8 @@ rule filter_tsv:
127127
rule tsv_to_vcf:
128128
threads: 1
129129
conda: "../envs/biopython.yaml"
130+
params:
131+
ref_name = config["ALIGNMENT_REFERENCE"],
130132
input:
131133
tsv = OUTDIR/f"{OUTPUT_NAME}.masked.prefiltered.tsv",
132134
output:

workflow/scripts/tsv_to_vcf.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -17,7 +17,7 @@ def tsv_to_vcf(tsv_file, vcf_file):
1717
# Process each row in the TSV file
1818
for index, row in tsv_df.iterrows():
1919
# Extract fields from the TSV row
20-
chrom = 'NC_045512.2'
20+
chrom = snakemake.params.ref_name
2121
pos = row['POS']
2222
ref = row['REF']
2323
alt = row['ALT']

0 commit comments

Comments
 (0)