This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filtering a vcf based on known variants

I have a variant table with hg19 genomic coordinates that I want to use to filter a vcf. What is a good way to go about filtering the vcf assuming I am only looking for lines where the genomic coordinates match between the vcf file and my known variants table?

vcf
• 2,687 views • link
written by idedios • idedios

1 answer

What are your file formats? You can use bedtools to intersect variants with coordinates, it can work with both -bed and -vcf files.

• 1 views • link
written by Alice Alice

As per Alice, you can use BEDTools; however, this will produce output that is not in VCF format.

If you still want to have the VCF format after you filter, then use bcftools filter --regions-file FILE (see Here for further information).

In both situations, your regions to filter should ideally be in BED format.

• 1 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this