Entering edit mode
9.4 years ago
MAPK
★
2.1k
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)?
Why not
bedtools intersect?