filter 1000 genome vcf file to keep common snps · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

filter 1000 genome vcf file to keep common snps

I want to filter "1000 genome vcf file" to keep only common SNPs. Can someone explain the difference between

1) Filtering by allele frequency AF greater than 0.05 (which I found yields more variants)

bcftools view -i 'AF>0.05'  variants.vcf -o common_variants.vcf

2) Filtering by MAF greater than 0.05

bcftools view -i 'MAF>0.05' variants.vcf -o common_variants.vcf

Which one is correct?

bcftools snps

• 1,007 views • link

updated 2.6 years ago by Ram

• written 2.6 years ago by shersky •

1 answer

AF is the frequency for an alternate allele. If AF < 0.5, then AF is equal to MAF. Rare variants generally has AF or MAF < 0.05. MAF (Minor allele frequency) refers to the minor allele (least frequent) frequency. Common SNPs have MAF > 0.05. Your filtering method 2 is the correct one to filter common SNPs.

bcftools view -i 'MAF>0.05' variants.vcf -o common_variants.vcf

• 0 views • link

written 2.6 years ago by bk11 •

Log in to answer this question.

More posts like this

The allele with index 2 is not defined in the REF/ALT columns" error after AF filtering with gnomAD…
written 1.2 years ago by erikatatianacs •

I'm running a germline variant calling pipeline with bcftools, and I annotate/filter variants based on allele frequency from gnomAD exomes (gnomad.exomes.r2.1.1.sites.vcf.bgz). After filtering, I get …
How to filter vcf file by MAF using bcftools?
written 2.9 years ago by _quantum_girl_ •

I saw this thread (https://github.com/samtools/bcftools/issues/357) and it seems sometimes bcftools recomputes MAF and sometimes it takes it from INFO. Which one is recommended to use …
Filtering and sorting of multiallelic sites
written 3.9 years ago by gernophil •

Hey everyone, since splitting of multiallelic sites the way bcftools does it, [doesn't really work for me](https://www.biostars.org/p/9537387/), I am looking for a good way to …
BCFtools allele frequency for specific population
written 4.4 years ago by Dean •

Hi, Apologies if this is a silly question or has an obvious answer. I am trying to output five VCF files using bcftools filtered by …
Minor Allele Frequency (MAF) Filtering differences: Gappet vs Tassel
written 5.1 years ago by waqaskhokhar999

Hi, My input genotypic file (Hapmap format) contains 48425 SNPs, I have tried to filter in Tassel V.5 based on Minor Allele frequency filtering at …
rare SNP or SNV
written 6.0 years ago by prasundutta87

Hi, Can a rare variant (point mutation) be called a rare SNP (single nucleotide polymorphism) or instead should be called a rare SNV (single nucleotide …
SNP subsetting using bcftools
written 6.3 years ago by alexcullsci •

I am trying to subset SNPs from 32 cultivars. I have the following exclusion criteria: 1. >20% missing data 2. SNPs present in less than …
bcftools / samtools, which database does the MAF function use?
written 7.3 years ago by Davide Chicco

Hi I am working with `bcftools view` to [filter variants][1] based on a specific MAF (minor allele frequency) cutoff. It works well but there's an …
G5/G5A flag, COMMON=Integer, and GMAF relationships from dbSNP
written 11.3 years ago by tlee23389 •

Off of dbSNP, it annotates with a couple of flags: ``` ##INFO=<ID=G5A,Number=0,Type=Flag,Description=">5% minor allele frequency in each and all populations"> ##INFO=<ID=G5,Number=0,Type=Flag,Description=">5% minor allele frequency in …
What is SIFT-score?
written 11.7 years ago by mangfu100

Hi. I am wondering a term with SIFT score. I think that SIFT refers to some measurement of SNPs, and while reading Annovar paper, I …