Question: How to remove regions from a vcf file ?
2
d0ct0r • 20 wrote:
I need to remove specific regions (in a bed file) from a vcf file. I am aware of the bcftools view options to filter regions. But It only subsets the regions. I need to do the opposite. I need to drop those regions from the vcf file. My solution is by subtracting the bed regions from vcf file regions using bedtools and then using the resulting bed file to filter the vcf file. But I think there will be some straightforward way to do it.
ADD COMMENT
• link
•
modified 4.5 years ago
by
andrew.j.skelton73 ♦ 5.9k
•
written
4.5 years ago by
d0ct0r • 20
Are you looking for something like this? With this command, you'll get all the variants that do not overlap with the regions specified in the bed file.
What is
-voption? It says unrecognized parameterSorry. My mistake. Try now the command. -v is the argument for printing all the vcf variants that DO NOT overlap with the given region, without -v, you'd get the overlapping variants.
intersectBed