Add "ALL" uniprot bed files to vcf files
0
0
Entering edit mode
7.2 years ago

Thanks to the help of all from this forum I know how to add the bed files from uniprot to my vcf file

ftp://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/genome_annotation_tracks/UP000005640_9606_beds/

I use this code

       bcftools annotate \
         -a ACTIV.bed.gz \
         -c CHROM,FROM,TO,UNIPROT  \
         -h <(echo '##INFO=<ID=UNIPROT,Number=1,Type=String,Description="name">') \
          sample.vcf.gz > sample_unipPROT.vcf

Easy to do with one file, but I would like to do with all the bed files from uniprot site. So at this point I have few questions...

  1. The best is to Merge all bed files toguether and annotate afther?
  2. Can I use bcftools to anotate all in one secuence? 3,Isnt there a bed file to download with all the info together?

Thanks

vcf uniprot bed • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 2076 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6