Question: vcftools: SNP filtering
1
gravatar for Parimala Devi
5.2 years ago by
India/Bangalore/Jubilant Biosys
Parimala Devi70 wrote:

I'm trying to filter my SNP vcf file. I'm looking for various filters that can be applied.

I've tried this, to remove the SNPs that have Genome quality < 50, but the number of files remain the same.

 vcftools --vcf input.vcf --recode --recode-INFO-all --out output.vcf --minGQ 50.00

How can I fix this? Also, I would like to know what other filters that I can use.

Thank you

tool snp vcftools filter • 6.3k views
ADD COMMENTlink modified 5.2 years ago by Dan D6.8k • written 5.2 years ago by Parimala Devi70
0
gravatar for Dan D
5.2 years ago by
Dan D6.8k
Tennessee
Dan D6.8k wrote:

To list here the possible filters you can pass into VCFTools would be an exercise in redundancy. The documentation does a good job of laying them all out.

You're specifying the genotype quality filter correctly--are you saying that you're getting the same number of lines in your output file? If you're not sure, here's a quick way to check the total number of lines in the input and output files:

wc -l input.vcf

wc -l output.vcf

You should also check to see if the "GQ" tag is actually in your VCF file. Depending on what software was used to generated the VCF, there might not be Genotype Quality data for all sites, which is what VCFTools requires..

EDIT: try changing --minGQ to --minQ

ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by Dan D6.8k

I did a word count before and after, it gives me the same number of SNPs.
I tried the same with the --minDP (depth), the number of SNPs remain a constant here as well.
I used bcftools to generate the VCF. Also, I checked for the GQ tag, it exists.

ADD REPLYlink written 5.2 years ago by Parimala Devi70

Would you be willing to post a dropbox (or google drive, or something similar) link to your VCF file?

ADD REPLYlink written 5.2 years ago by Dan D6.8k

Sure.

Here's a link to it. |

https://drive.google.com/file/d/0B5OyhrxgtH2yU0lFWUVTTE1kbk0/edit?usp=sharing

ADD REPLYlink written 5.2 years ago by Parimala Devi70

Thanks! Answer updated.

ADD REPLYlink written 5.2 years ago by Dan D6.8k
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: 1378 users visited in the last hour