This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to determine the number of SNPs captured by a list of genes in VCF file

I have a bed file called mybedfile.bed which has a defined region for a list of genes. I have VCF files split into 24 chromosomes. I want to determine the number of SNPs captured by the regions defined in mybedfile.bed. What would be the easiest way to get this done (preferably in R)?

snps

Why not bedtools intersect?

0 answers

No answers yet.

Log in to answer this question.