This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VCFtools filtration strategy --hwe

Hi everyone,

I'm currently confused by the option of VCFtools: --hwe. From my understanding, this option measure whether each site from every individual in the VCF file significantly departures from the HW balance. But what if individuals in the file are not from the same population or even species? They definitely cannot reach HW balance for they breed separately. In this case what's the meaning of it?

Thanks!

snp genome next-gen sequencing
• 2,876 views • link
updated by pufferfish pufferfish • written by chenyangkang24 • chenyangkang24

1 answer

You are correct. This option calculates the HWE score and then filters sites according to the threshold you provide.

It would not be a valid filter if any of the following:

  • different species (??)
  • different populations
  • recently admixed populations
  • a variant under strong selection

This is a common challenge in bioinformatics: just because an option exists doesn't mean it applies to your particular dataset.

• 0 views • link
written by pufferfish pufferfish

Thank you,

That's exactly what I thought. I asked this "stupid" question because some tutorial add this option without explanation, even if the samples comes from different species. In this case many SNPs were gone.

Yangkang

• 0 views • link
written by chenyangkang24 • chenyangkang24

Log in to answer this question.

More posts like this