Question: Fst calculation between my own population and 1000 genome populations
gravatar for seta
13 months ago by
seta1.2k 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 • 646 views
ADD COMMENTlink modified 12 months ago by cetin.m30 • written 13 months ago by seta1.2k
gravatar for cetin.m
12 months ago by
cetin.m30 wrote:

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

ADD COMMENTlink written 12 months ago by cetin.m30

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 months ago • written 12 months ago by seta1.2k

Hi.. Even I would like to have your script. It would be really great if you could share'

Thanks in Advance.

ADD REPLYlink written 8 months ago by shreyajha30
Please log in to add an answer.


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