Entering edit mode
8.5 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
?