Question: How to replicate vcftool's sliding window SNP density?
0
gravatar for serpalma.v
12 months ago by
serpalma.v40
Germany
serpalma.v40 wrote:

Hello!

vcftools can calculate statistics such as pi and Fst using sliding windows.

I've noticed that windows in which the alternative allele is fixed are not taken into account. However, those appear in the VCF file.

For example, N_VARIANTS for pi at window_x is 20, but the region in the vcf file reports 25 SNPs, 5 of which are hom-alt in all samples.

I would like to find a way to calculate the SNP density as it comes in the vcf file using the same boundaries reported by vcftools.

Is there a way to make vcftools to do this?

Thanks!

sequencing snp vcftools • 559 views
ADD COMMENTlink modified 12 months ago by Vitis2.4k • written 12 months ago by serpalma.v40
1
gravatar for Vitis
12 months ago by
Vitis2.4k
New York
Vitis2.4k wrote:

Create a window file in BED format and use 'bedtools intersect' to intersect VCF and the window file. I believe there is an option to report number of overlaps, which in your case, would be number of SNPs in the windows. Please see:

https://bedtools.readthedocs.io/en/latest/

ADD COMMENTlink written 12 months ago by Vitis2.4k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1198 users visited in the last hour