diff --git a/pipelines/annotations.snakefile b/pipelines/annotations.snakefile index 54b1e0ca..07d28df8 100644 --- a/pipelines/annotations.snakefile +++ b/pipelines/annotations.snakefile @@ -215,8 +215,8 @@ rule absplice: vcf=expand( [ absplice_repo_dir - / "example/data/resources/analysis_files/vcf_files" - / (vcf_pattern), + / "example/data/resources/analysis_files/input_files" + / (vcf_pattern + "_variants_header.vcf"), ], zip, chr=chromosomes, @@ -253,14 +253,13 @@ rule mod_config_absplice: rule link_files_absplice: input: - anno_tmp_dir / (vcf_pattern + "_variants_header.vcf.gz"), + anno_tmp_dir / (vcf_pattern + "_variants_header.vcf"), output: absplice_repo_dir - / "example/data/resources/analysis_files/vcf_files" - / (vcf_pattern), + / "example/data/resources/analysis_files/input_files" + / (vcf_pattern + "_variants_header.vcf"), shell: - " ".join(["ln", "-s", "-r", "{input}", "{output}"]) - + f"mkdir -p {absplice_repo_dir/'example/data/resources/analysis_files/input_files'} && ln -s -r {{input}} {{output}}" @@ -345,7 +344,7 @@ rule concat_deepSea: rule deepSea: input: - variants=anno_tmp_dir / (vcf_pattern + "_variants_header.vcf.gz"), + variants=anno_tmp_dir / (vcf_pattern + "_variants_header.vcf"), fasta=fasta_dir / fasta_file_name, output: anno_dir / (vcf_pattern + ".CLI.deepseapredict.diff.tsv"), @@ -451,7 +450,7 @@ rule extract_with_header: input: bcf_dir / (vcf_pattern + ".bcf"), output: - anno_tmp_dir / (vcf_pattern + "_variants_header.vcf.gz"), + anno_tmp_dir / (vcf_pattern + "_variants_header.vcf"), shell: ( load_bfc diff --git a/pipelines/resources/absplice_config.yaml b/pipelines/resources/absplice_config.yaml index 9bf28a39..ea35c1ad 100644 --- a/pipelines/resources/absplice_config.yaml +++ b/pipelines/resources/absplice_config.yaml @@ -3,8 +3,8 @@ # specify genome version hg19 or hg38 genome: hg38 -# store all the vcf file that you want to analyze to directory ../data/resources/vcf_files/ -vcf: ../data/resources/analysis_files/vcf_files/{vcf_id} +# store all the vcf file that you want to analyze to directory ../data/resources/input_files/ +vcf: ../data/resources/analysis_files/input_files/{vcf_id} # parameters for variant filtering variant_filtering: