Skip to content
Snippets Groups Projects
Commit acb110be authored by mariabernard's avatar mariabernard
Browse files

split_vcf need that input vcf is indexed with tabix

parent 445992e0
No related branches found
No related tags found
No related merge requests found
......@@ -18,10 +18,21 @@ def get_input_vcf(wildcards) :
for Set in config['vqsr_resources']:
if os.path.basename(Set['file']).startswith(wildcards.prefix):
return Set['file']
def get_input_index(wildcards) :
if wildcards.prefix == "common" :
return "results/intersect/common_variants.vcf.gz.tbi"
if wildcards.prefix == "VQSR_trainingSet_untrusted_uniquely_called_variants":
return "results/intersect/VQSR_trainingSet_untrusted_uniquely_called_variants.vcf.gz.tbi"
else:
for Set in config['vqsr_resources']:
if os.path.basename(Set['file']).startswith(wildcards.prefix):
return Set['file'] + ".tbi"
rule split_vcf:
input :
vcf = get_input_vcf ,
tbi = get_input_index ,
ref = config["reference_genome"]
output :
vcf = temp("results/split_SNP_INDEL/{prefix}_{var}.vcf.gz")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment