How to filter out SNPs with Minor allele frequency less than 5% in VCF file

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)

See less

jdpr_100

Member

Join Date: Jun 2014

Posts: 14
- Share
- Tweet
#1

How to filter out SNPs with Minor allele frequency less than 5% in VCF file

09-30-2014, 06:40 AM

Hello,
I would like to filter out SNPs with allele frequencies in VCF file from multi-sample variant calling.
I tried VCFtools with the following commands, but it doesn't work.

vcftools --vcf input.vcf --recode --recode-INFO-all --stdout --maf 0.05 --max-maf 1 --out outpur.vcf

Here is the vcf file format generated from samtools mpileup.

#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 788 712 908 544 541 581 801 381 479 487 379 570
Chr01 11356 . C A 124 . DP=88;VDB=1.624097e-04;RPB=-2.346857e+00;AF1=0.1553;AC1=3;DP4=75,3,10,0;MQ=30;FQ=124;PV4=1,1,1,1
GT:PLP:GQ 0/0:0,0,0:0:5 0/0:0,6,59:2:11 0/0:0,0,0:0:5 0/1:28,0,123:8:24 1/1:157,24,0:8:14 0/0:0,39,221:13:43 0/0:0,12,85:4
:17 0/0:0,36,164:12:40 0/0:0,39,147:13:43 0/0:0,9,71:3:14 0/0:0,48,163:16:52 0/0:0,27,126:9:31
Chr01 11384 . T C 15.9 . DP=123;VDB=5.617680e-03;RPB=1.820871e+00;AF1=0.04721;AC1=1;DP4=114,3,5,0;MQ=31;FQ=16.1;PV4=1,1,0.05,1 GT:PLP:GQ 0/0:0,12,79:4:22 0/0:0,6,60:2:16 0/0:0,0,0:0:10 0/0:0,30,163:10:40 0/0:0,27,159:9:37 0/0:0,60,244:20:70 0/0:0,12,85:4:22 0/0:0,48,163:16:58 0/1:56,0,116:15:46 0/0:0,18,122:6:28 0/0:0,69,185:23:79 0/0:0,39,135:13:49

I wonder if anyone could help me
Thanks
Tags: None
swbarnes2

Senior Member

Join Date: May 2008

Posts: 910
- Share
- Tweet
#2

10-01-2014, 02:14 PM

The DP4 part spells out how many reads aligned to reference and alternate allele. So you can parse based on that. You could do that with some finagling in Excel, if you had to.
Comment

Previous template Next

A Review of Next-Generation Sequencing Methods for Studying Epigenetics—Part 1

by seqadmin

Developments in sequencing technologies and methodologies have transformed the field of epigenetics, giving researchers a better way to understand the complex world of gene regulation and heritable modifications. This article explores some of the diverse sequencing methods employed in the study of epigenetics, ranging from classic techniques to cutting-edge innovations while providing a brief overview of their processes, applications, and advances.

Methylation Detect...
- Channel: Articles
05-31-2023, 10:46 AM
Differential Expression and Data Visualization: Recommended Tools for Next-Level Sequencing Analysis

by seqadmin

After covering QC and alignment tools in the first segment and variant analysis and genome assembly in the second segment, we’re wrapping up with a discussion about tools for differential gene expression analysis and data visualization. In this article, we include recommendations from the following experts: Dr. Mark Ziemann, Senior Lecturer in Biotechnology and Bioinformatics, Deakin University; Dr. Medhat Mahmoud Postdoctoral Research Fellow at Baylor College of Medicine;...
- Channel: Articles
05-23-2023, 12:26 PM
Variant Analysis and Genome Assembly: Recommended Tools for Next-Level Sequencing Analysis

by seqadmin

Continuing from our previous article, we share variant analysis and genome assembly tools recommended by our experts Dr. Medhat Mahmoud, Postdoctoral Research Fellow at Baylor College of Medicine, and Dr. Ming "Tommy" Tang, Director of Computational Biology at Immunitas and author of From Cell Line to Command Line.

Variant detection and analysis tools
Mahmoud classifies variant detection work into two main groups: short variants (<50...
- Channel: Articles
05-19-2023, 10:03 AM

Topics	Statistics	Last Post
Genetic Discoveries Transform Our Understanding of Primate Diversity and Behavior by seqadmin Started by seqadmin, 06-01-2023, 08:56 PM	0 responses 11 views 0 likes	Last Post by seqadmin 06-01-2023, 08:56 PM
Deep Sequencing Unearths Novel Genetic Variants: Enhancing Precision Medicine for Vascular Anomalies by seqadmin Started by seqadmin, 06-01-2023, 07:33 AM	0 responses 14 views 0 likes	Last Post by seqadmin 06-01-2023, 07:33 AM
Unveiling Genetic Associations Through Transcription Factor Binding Quantitative Trait Loci by seqadmin Started by seqadmin, 05-31-2023, 07:50 AM	0 responses 52 views 0 likes	Last Post by seqadmin 05-31-2023, 07:50 AM
Exploring French-Canadian Ancestry: Insights into Migration, Settlement Patterns, and Genetic Structure by seqadmin Started by seqadmin, 05-26-2023, 09:22 AM	0 responses 56 views 0 likes	Last Post by seqadmin 05-26-2023, 09:22 AM

Header Leaderboard Ad

How to filter out SNPs with Minor allele frequency less than 5% in VCF file

Announcement

SEQanswers June Challenge Has Begun!

How to filter out SNPs with Minor allele frequency less than 5% in VCF file

Comment

Latest Articles

ad_right_rmr

News