Question: vcf filter of lowQual
0
gravatar for genomics Newbie
2.0 years ago by
genomics Newbie20 wrote:

How do I filter my vcf so that all "lowQual" entries are excluded from the variant file.

When I ran the following command I included the -e flag which I thought would exclude all LowQual, but they appear in my vcf and I need to remove them so that I can use my vcf as further input into another pipeline package.

bcftools mpileup -Ou -f /path/to/ucsc.hg19.fasta example.exome.bam | \
    bcftools call -Ou -mv | \
    bcftools filter -s LowQual -e '%QUAL<20 || DP>100' > var.flt.vcf

It did not appear that I can use vcftools to only include "Passed" in the vcf file.

Thank you.

vcf filter • 2.8k views
ADD COMMENTlink modified 2.0 years ago by Kevin Blighe63k • written 2.0 years ago by genomics Newbie20

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

ADD REPLYlink written 2.0 years ago by RamRS28k
2
gravatar for Kevin Blighe
2.0 years ago by
Kevin Blighe63k
Kevin Blighe63k wrote:

To only include those with PASS set in the FILTER column, use:

bcftools view -f PASS

Kevin

ADD COMMENTlink modified 2.0 years ago • written 2.0 years ago by Kevin Blighe63k
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: 1356 users visited in the last hour