SNPs that are at least "N" bp away from each other

[Rhewter]

Hello,

I got a set of SNPs in a vcf file and would like to discard SNPs many nearby. The idea is to take only SNPs that are at least 200 bp away from each other. How can I do this?

Thanks in advance.

[dpryan]

You could do this with awk, assuming the VCF file is sorted in a useful way. The general idea is to keep track of two variables: the chromosome of the last entry and the position of the last entry. If the chromosome of an entry and the previous are different or their distances are above 200bp then print $0 (the variables should be updated in any case). Have a go at writing such a script in awk (or python, or perl, or whatever else you know) and post a question when you run into problems.

[Rhewter]

You could do this with awk, assuming the VCF file is sorted in a useful way. The general idea is to keep track of two variables: the chromosome of the last entry and the position of the last entry. If the chromosome of an entry and the previous are different or their distances are above 200bp then print $0 (the variables should be updated in any case). Have a go at writing such a script in awk (or python, or perl, or whatever else you know) and post a question when you run into problems.

Hi...

Thank you for your answer. I thought that possibility before but would like to know whether there was a ready-made option. But I think your suggestion is the best option.

Thank you!

[SNPsaurus]

vcftools has a --thin option

--thin <int>

Thin sites so that no two sites are within the specified distance from one another.

[Rhewter]

vcftools has a --thin option

--thin <int>

Thin sites so that no two sites are within the specified distance from one another.

Hi SNPsaurus,

I tried using this function but it only generates a log file (with the number of unique variants that range) and not a new VCF. Is that so?


Thank you for your answer!

[SNPsaurus]

I think you need to add the --recode option, which tells vcftools to actually make a new vcf based on the filtering.

[Rhewter]

I think you need to add the --recode option, which tells vcftools to actually make a new vcf based on the filtering.

It really is that. Thank you very much once again.

[[email protected]]

Hi SNPsaurus,
I have applied various filters from VCFfilter tool i.e. DP>10, AF>0.1 etc.
This filter information is shown in top of the VCF file also. However, when I used to work with VCFtools for minimum SNPs distance (--Thin 10000) parametere, it does not show the information in vcf file.
Can you make me understand regarding this?

Regards
Prakash Thakor