This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bcftools / samtools, which database does the MAF function use?

Hi

I am working with bcftools view to filter variants based on a specific MAF (minor allele frequency) cutoff.

It works well but there's an information I am missing: which database is bcftools using to apply the MAF filter, that is to claim if a variant is rare or not?

Thanks!

-- Davide

bcftools samtools maf
• 2,950 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Davide Chicco Davide Chicco

1 answer

there is no database associated de facto with bcftools. It uses the AC/AN/AF from the genotypes in the VCF.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Hi Pierre

Thanks for the reply. Are you saying the information about if a variant is rare (or not) is present in the header of the VCF file we're using?

If yes, how does bcftools view retrieve it?

And, more important, how can we know which database did the people who generated the file used to state which variants are are and which variants are not?

Merci!

• 0 views • link
written by Davide Chicco Davide Chicco

Log in to answer this question.

More posts like this