浏览代码

test

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

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

# vcf2vcf: transfer into a standardized format # vcf2vcf: transfer into a standardized format
echo "Transfer the VCF file 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 raw.vcf --output-vcf norm.vcf \
--input-vcf $raw.vcf --output-vcf norm.vcf \
$SAMPLE_vcf2vcf \ $SAMPLE_vcf2vcf \
--ref-fasta ${ref_dir}/${fasta} --ref-fasta ${ref_dir}/${fasta}
if [ $normal_id ]; then if [ $normal_id ]; then
cp norm.vcf ${basename}.norm.vcf
cp $norm.vcf ${basename}.norm.vcf
else else
cut -f 1,2,3,4,5,6,7,8,9,10 norm.vcf > ${basename}.norm.vcf
cut -f 1,2,3,4,5,6,7,8,9,10 $norm.vcf > ${basename}.norm.vcf


# VEP annotation # VEP annotation
--fork $nt --fork $nt


# vcf2maf # vcf2maf
echo "vcf2maf..." echo "vcf2maf..."
perl ${vcf2maf_path}/vcf2maf.pl \ perl ${vcf2maf_path}/vcf2maf.pl \
$SAMPLE_vcf2maf \ $SAMPLE_vcf2maf \
--ref-fasta ${ref_dir}/${fasta} \ --ref-fasta ${ref_dir}/${fasta} \
--ncbi-build ${ncbi_build} \ --ncbi-build ${ncbi_build} \
--species ${species} \
--vep-fork $nt
--species ${species}
>>> >>>
runtime { runtime {

正在加载...
取消
保存