浏览代码

Update: add annotation

master
YaqingLiu 3 年前
父节点
当前提交
3acd950c74
共有 1 个文件被更改,包括 5 次插入2 次删除
  1. +5
    -2
      tasks/VEP.wdl

+ 5
- 2
tasks/VEP.wdl 查看文件



awk -F'\t' '{if(($1~"^#")||($1!~"^#" && $7=="PASS")){print $0}}' ${vcf} > ${sample_id}.vcf awk -F'\t' '{if(($1~"^#")||($1!~"^#" && $7=="PASS")){print $0}}' ${vcf} > ${sample_id}.vcf


# Judge the SAMPLE info of vcf file # Judge the SAMPLE info of vcf file
tumor_id=`awk -F'\t' '{if($1~"^#CHROM"){print $10}}' ${vcf}` tumor_id=`awk -F'\t' '{if($1~"^#CHROM"){print $10}}' ${vcf}`
normal_id=`awk -F'\t' '{if($1~"^#CHROM"){print $11}}' ${vcf}` normal_id=`awk -F'\t' '{if($1~"^#CHROM"){print $11}}' ${vcf}`
fi fi


# vcf2vcf: transfer into a standardized format # vcf2vcf: transfer into a standardized format
echo "Transfer the VCF file into a standardized format..."
perl ${vcf2maf_path}/vcf2vcf.pl \ perl ${vcf2maf_path}/vcf2vcf.pl \
--input-vcf ${sample_id}.vcf --output-vcf ${basename}.norm.vcf \ --input-vcf ${sample_id}.vcf --output-vcf ${basename}.norm.vcf \
$SAMPLE_vcf2vcf \ $SAMPLE_vcf2vcf \
--ref-fasta ${ref_dir}/${fasta} --ref-fasta ${ref_dir}/${fasta}


# VEP annotation # VEP annotation
echo "VEP annotation..."
perl ${vep_path}/vep --format vcf --vcf \ perl ${vep_path}/vep --format vcf --vcf \
--assembly ${ncbi_build} \ --assembly ${ncbi_build} \
--species ${species} \ --species ${species} \
$buffer_size \ $buffer_size \
--input_file ${basename}.norm.vcf --output_file ${basename}.vep.vcf --input_file ${basename}.norm.vcf --output_file ${basename}.vep.vcf


# vcf2maf # vcf2maf
echo "vcf2maf..."
perl ${vcf2maf_path}/vcf2maf.pl \ perl ${vcf2maf_path}/vcf2maf.pl \
--inhibit-vep \ --inhibit-vep \
--input-vcf ${basename}.vep.vcf --output-maf ${basename}.maf \ --input-vcf ${basename}.vep.vcf --output-maf ${basename}.maf \

正在加载...
取消
保存