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

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

bcftools view -f PASS

Kevin

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