Question: vcf filter of lowQual
0
gravatar for genomics Newbie
20 months 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.2k views
ADD COMMENTlink modified 20 months ago by Kevin Blighe56k • written 20 months 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 20 months ago by RamRS26k
2
gravatar for Kevin Blighe
20 months ago by
Kevin Blighe56k
Kevin Blighe56k wrote:

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

bcftools view -f PASS

Kevin

ADD COMMENTlink modified 20 months ago • written 20 months ago by Kevin Blighe56k
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: 1073 users visited in the last hour