Question: Fst calculation between my own population and 1000 genome populations
0
gravatar for seta
12 weeks ago by
seta1.1k
Sweden
seta1.1k wrote:

Hi all and Merry Christmas!

I would like to calculate Fst for some selected SNPs between my own population (the variants came from whole-genome sequencing of about 800 individuals) and other populations in 1000 genome. To this end, I extracted the variants of interest using VCFtools with --recode --keep-INFO-all in the command from each chromosome VCF file, then I concatenated all extracted vcf files using bcftools to have the single vcf file with all chromosomes within it. So, now, I have the single vcf file containing variants of interest in all chromosomes from my population and another similar file from 1000 genome. Now, for calculating Fst between my population and 1000 genome population, these two files should be merged, yes? or is there a possibility to calculate Fst between two separate vcf file? However, I am not sure how they should be appropriately merged for Fst calculation by the available tools, say vcftools. Please share me your suggestions and comments.

Thank you in advance

fst 1000 genome population • 264 views
ADD COMMENTlink modified 12 weeks ago by cetin.m20 • written 12 weeks ago by seta1.1k
0
gravatar for cetin.m
12 weeks ago by
cetin.m20
cetin.m20 wrote:

When I did something similar I used a python code which did not involve merging the files.

ADD COMMENTlink written 12 weeks ago by cetin.m20

I wanted use vcftools for fst calculation, I merged vcf files using bcftools, but couldn't successfully calculate fst on the merged vcf file. Could you please share your script?

ADD REPLYlink modified 12 weeks ago • written 12 weeks ago by seta1.1k
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: 1394 users visited in the last hour