Question: vcf filter of lowQual
0
gravatar for genomics Newbie
8 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 • 844 views
ADD COMMENTlink modified 8 months ago by Kevin Blighe41k • written 8 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 8 months ago by RamRS21k
2
gravatar for Kevin Blighe
8 months ago by
Kevin Blighe41k
London, England
Kevin Blighe41k wrote:

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

bcftools view -f PASS

Kevin

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