exclude variants based on position in vcf file
1
0
Entering edit mode
4.5 years ago

I have a list of positions (POS) in a text file I would like to remove from a .vcf file. What is the simplest way to do this? I have tried:

vcftools --vcf cluster_1_private.vcf --exclude remove_poly_1.txt --recode --recode-INFO-all --out cluster_1_fixed

but this doesnt seem to work, but note the text file only contains a list of numbers such as:

3416

14547

55489

56775

66521

91512

Maybe I need to specify the POS somehow?

SNP • 2.2k views
ADD COMMENT
1
Entering edit mode
4.5 years ago

Figure it out:

vcftools --vcf cluster_1_private.vcf --exclude-positions remove_poly_1.txt --recode --recode-INFO-all --out cluster_1_fixed

text file needs to be a tab-separated list of chromosome and positions

ADD COMMENT

Login before adding your answer.

Traffic: 1922 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