Przeglądaj źródła

Fixbug: vcf2vcf need --vcf-normal-id

master
YaqingLiu 3 lat temu
rodzic
commit
b04b9423ab
1 zmienionych plików z 14 dodań i 5 usunięć
  1. +14
    -5
      tasks/VEP.wdl

+ 14
- 5
tasks/VEP.wdl Wyświetl plik

if [ $normal_id ]; then if [ $normal_id ]; then
SAMPLE_vcf2maf="--tumor-id $tumor_id --normal-id $normal_id" SAMPLE_vcf2maf="--tumor-id $tumor_id --normal-id $normal_id"
SAMPLE_vcf2vcf="--vcf-tumor-id $tumor_id --vcf-normal-id $normal_id" SAMPLE_vcf2vcf="--vcf-tumor-id $tumor_id --vcf-normal-id $normal_id"
${vcf} > raw.vcf
else else
SAMPLE_vcf2maf="--tumor-id $tumor_id" SAMPLE_vcf2maf="--tumor-id $tumor_id"
SAMPLE_vcf2vcf="--vcf-tumor-id $tumor_id"
# Add a column and remove it after vcf2vcf
SAMPLE_vcf2vcf="--vcf-tumor-id $tumor_id --vcf-normal-id $tumor_id"
awk -F'\t' 'OFS="\t" {if($1!~"^##" && length($11)==0) $11=$10; print $0}' ${vcf} > raw.vcf
fi fi




# 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 ${vcf} --output-vcf ${basename}.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
${basename}.norm.vcf=norm.vcf
else
cut -f 1,2,3,4,5,6,7,8,9,10 norm.vcf > ${basename}.norm.vcf

# VEP annotation # VEP annotation
echo "VEP annotation..." echo "VEP annotation..."
perl ${vep_path}/vep --format vcf --vcf \ perl ${vep_path}/vep --format vcf --vcf \
--input_file ${basename}.norm.vcf --output_file ${basename}.vep.vcf \
--assembly ${ncbi_build} \ --assembly ${ncbi_build} \
--species ${species} \ --species ${species} \
--everything --af_exac \ --everything --af_exac \
--cache --dir_cache ${cache} \ --cache --dir_cache ${cache} \
--fasta ${ref_dir}/${fasta} \ --fasta ${ref_dir}/${fasta} \
$buffer_size \ $buffer_size \
--fork $nt \
--input_file ${basename}.norm.vcf --output_file ${basename}.vep.vcf
--fork $nt


# vcf2maf # vcf2maf

Ładowanie…
Anuluj
Zapisz