Calculating Fst from low coverage data using genotype likelihoods
0
1
Entering edit mode
3.6 years ago
Jautis ▴ 580

Hello, I have low coverage (~1x) sequencing data from ~50 individuals from two populations and would like to calculate Fst for each biallelic variant. I can calculate allele frequencies based on the genotype likelihoods and estimate Fst following Wright 1951, but this assumes infinite populations which accurately sample allele frequencies. Alternately, I can use vcftools to estimate Fst following Weir & Cockerman 1984, but this implementation uses 0/1/2 genotypes while accounting for finite sample sizes (i.e. doesn't incorporate genotyping error). Is there a method that incorporates both genotype error and sample sizes? Alternately, does it seem kosher to use genotype likelihoods as the allele estimate per individual and the method of Weir & Cokerman?

Thank you in advance!

vcftools fst population_differentiation • 675 views
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