I am trying to calculate Pooled Heterozygosity (hp) by identifying nMAJ and nMIN from vcf file with a sliding window 150kb. I am confused after reading papers where they calculated using formula but no particular method to calculate nMAJ and nMIN. Can you please share the way to calculate the same. I will be grateful to you all.
Hp = 2ΣnMAJΣnMIN/(ΣnMAJ + ΣnMIN)2
Thankyou
•
written
by
tothepoint
1 answer
Hi Devarora,
As far as I understand, you have to find, for every SNP the most common allele (nMAJ). You sum the count of all these major alleles in your 150 kb window. Same goes for the least common allele at every SNP (nMIN). again, you sum them in the window.
You can extract the count of alleles from the vcf (https://www.internationalgenome.org/wiki/Analysis/vcf4.0/). For this, you will probably need to program something in bash, awk, perl, python, to extract the right column and retrieve the numbers that you need to sum.
Best,
Guenole
Log in to answer this question.
-
Exact steps to calculate the SNP index to identify locality of EMS induced mutation
written by n.anuragsharma •
-
Different insertion times from two methods
written by hashim.rana11 •
-
pooled-heterozygosity calculation
written by reza
-
The pooled heterozygosity (Hp) score
written by Siavash Salek Ardestani •
-
Question regarding deFuse and fusion calling
written by CY
-
What is signal to interference (S2I) and fraction quantile normalization in LC-MS/MS
written by agicict
-
Qst Values calculation
written by biobudhan •
-
Illumina Exome Chip And Pathway Analysis
written by alizaid12 •
-
Phenotypic Variation Explained By Snps In Gwas
written by kumar.vinod81
-
Overlapping Reads In Window Approach
written by Vikas Bansal
Hi to thepoint
I have same question. It's been about three months since you asked the question. If you have received an answer during this time, please share it with us so that we can use it as well.
Thankful