본문 바로가기

Bioinformatics(생정보학)

vep 특정 genome 및 gtf파일 사용하기

728x90
반응형

특정 gtf와 reference genome을 다운로드 받음


0. HTS lib다운로드

http://www.htslib.org/download/

tar -xvf htslib-1.5.tar.bz2

cd htslib-1.5

make

sudo make install # 또는 ./bgzip 을 쓰면 사용가능함.


GTF사용 준비

/home/sjoh/hg38/htslib/htslib-1.5


grep -v "#" data.gff | sort -k1,1 -k4,4n -k5,5n -t$'\t' | bgzip -c > data.gff.gz
tabix -p gff data.gff.gz
./vep -i input.vcf -gff data.gff.gz -fasta genome.fa.gz

728x90
반응형