Question: vcf SNP file minor allele frequency filter?
1
gravatar for sun.nation
3.7 years ago by
sun.nation120
United States
sun.nation120 wrote:

I want to filter vcf file with minor allele frequency of 30% = 0.3. I tried with: --maf 0.3

vcftools --vcf ps-pt.bcf.vcf --out ps-pt.bcf.vcf-5dp-miss0 --remove-filtered-all --maf .3 --max-missing-count 0 --minDP 5 --minQ 10 --minGQ 10 --recode --recode-INFO-all

But, when I checked the alignment and some het positions. Het was called such as G/T (1 read with G and 9 reads with T). I want like this (3 out of 10 reads should be G to be called GT in minimum)

I want to filter per individual not by combining all.Thanks for your suggestions.

 

 

vcftools minorallele snp allele vcf • 3.1k views
ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by sun.nation120

I think you are confusing minor allele frequency (the frequency of the rare allele observed among different samples) with allele balance (the ratio of reads showing the reference sequence versus the reads showing the alt sequence, at a given site, in one individual).

I am unsure if vcftools has a filter option for allele balance. I know GATK does:

https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_annotator_AlleleBalanceBySample.php

 

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by rbagnall1.4k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1867 users visited in the last hour