SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Calculating allele frequency from VCF file bio_informatics Bioinformatics 3 10-25-2016 08:50 AM
Filter VCF using external Allele Frequency rdoan Bioinformatics 1 10-01-2014 02:08 AM
how to compute Minor Allele Frequency for a variant call? BhariD Bioinformatics 0 05-07-2013 11:50 AM
Question on Minor allele Frequency in dbSNP135 ngarg General 0 02-06-2012 10:55 AM
Minor Allele Frequency Cutoff/Threshold BertieWooster Bioinformatics 1 09-28-2009 02:50 PM

Reply
 
Thread Tools
Old 09-30-2014, 07:40 AM   #1
jdpr_100
Member
 
Location: Gaineville FL

Join Date: Jun 2014
Posts: 14
Default How to filter out SNPs with Minor allele frequency less than 5% in VCF file

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.
Quote:
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.
Quote:
#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
jdpr_100 is offline   Reply With Quote
Old 10-01-2014, 03:14 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

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.
swbarnes2 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:22 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO