melodyeeeee

ERROR: You're either using an outdated samtools, or --ref-fasta is not the same genome build as your --input-vcf.

2

Hi, When I running this command : `vcf2maf=/data/jiaf/wes_cancer/biosoft/vcf2maf/vcf2maf.pl vep_path=~/miniconda3/envs/vep/bin/ GENOME=/data/jiaf/wes_cancer/biosoft/gsutil/gatk/hg19/v0/Homo_sapiens_assembly19.fasta vep_data=/data/jiaf/wes_cancer/data/vep/homo_sapiens/100_GRCh37 cat config | while read id do arr=(${id}) tumor=${arr[0]} perl ${vcf2maf} \ --input-vcf filter/clean/${tumor}.clean.vcf \ --output-maf annotation/vep/${tumor}_vep.maf \ --ref-fasta...