Thanks for the answer.
If I do this, I would need a dataset with genes name and their coordinates wouldn't I?
I have a VCF files with big deletions and I want to know if these deletion are in (completly or partially) of a particular gene. I want to basically know the gene names affected by these deletions
Simply annotate your VCF file using any of the program; ChIPseeker, MEME, bcftools ... the list is long. I am demonstrating bcftools here:
bcftools annotate -a genes.bed.gz -c CHROM,FROM,TO,GENE -h <(echo '##INFO=ID=GENE,Number=1,Type=String,Description="Gene name">') variants.vcf.gz
A few things to note:
what have you tried so far?
quick search should give you half the solution already (bedtools intersect)
Log in to answer this question.