diff --git a/snakecnv/detection.snk b/snakecnv/detection.snk index 60a8f0015e2e7207d6febc50f101c2895bc26be0..785c50af9a5196a34dd7ef3ccd4d4778692691bb 100644 --- a/snakecnv/detection.snk +++ b/snakecnv/detection.snk @@ -91,8 +91,9 @@ tools = config['tools'] # list of chromosomes chromosomes = config["chromosomes"] chr_batches = {} +ref_chr = "reference_raw.fasta" if not os.path.exists(REFERENCE) and os.path.exists("reference_raw.fasta") else REFERENCE for chromosome in chromosomes: - chr_batches[chromosome] = get_chr_batches(REFERENCE, chromosome) + chr_batches[chromosome] = get_chr_batches(ref_chr, chromosome) # list of variants to be detected varianttypes = ['DEL', 'INV', 'DUP', 'mCNV']