Filter VCF file by list of sites.
1
0
Entering edit mode
8 months ago
selplat21 ▴ 20

I currently have a list of a few million sites I would like to extract from a larger VCF file. I have been using the following format:

module load bcftools tabix -R sites.list Sample.vcf.gz > Sample.subset.vcf.gz

This works, but unfortunately has been taking over a week to run. Does anyone have a faster solution?

vcf filtering genomics • 247 views
ADD COMMENT
2
Entering edit mode
8 months ago
bcftools view --threads 10 --regions-file sites.bed -O b -o samples.subset.bcf samples.vcf.gz
ADD COMMENT

Login before adding your answer.

Traffic: 1501 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6