Question: Add "ALL" uniprot bed files to vcf files
0
gravatar for cristina_sabiers
2.8 years ago by
Spain
cristina_sabiers70 wrote:

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

uniprot bed vcf • 851 views
ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by cristina_sabiers70
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2096 users visited in the last hour