Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
MiSeq gDNA reads still fail "Kmer content" and "per base seq content" after trimming" ysnapus Illumina/Solexa 4 11-12-2014 07:25 AM
Samtools and bcftools: Why aren't some of the calls "missing data". molecules Bioinformatics 0 01-10-2014 10:07 AM
"Bcftools mpileup" removes marked duplicates? angerusso Bioinformatics 0 03-05-2012 06:43 AM
BAM flags using "OR" RockChalkJayhawk Bioinformatics 1 02-06-2012 01:49 PM
comparing two samples using "bcftools -T" option? NanYu Bioinformatics 0 10-25-2011 07:01 AM

Thread Tools
Old 06-27-2016, 07:37 AM   #1
Junior Member
Location: USA

Join Date: Aug 2015
Posts: 7
Default confusion about BCFtools "EXPRESSIONS" flags

I'm hoping someone can help me understand the "EXPRESSIONS" flags used in bfctools.

I'm trying to filter variants from a vcf file produced using the 'mpileup' program from samtools. Essentially, I want to perform a simple filtering step that will exclude all reads that don't pass a threshold (say, > 30% of the reads sequenced at that position)

This is what I I'm using now:
bcftools filter -i'%MIN(DV/DP)>.3' --output-type v -o mpileup.vcf

However, I'm a bit confused about the proper usage of "MIN".

According the bfctools documentation page, "MIN" should be used with "on FORMAT tags (over samples) and INFO tags (over vector fields)" see

However, I have no clue what this means. I'm not sure if I should drop the MIN and just filter using -i'%(DV/DP)>0.3'.
lstbl is offline   Reply With Quote

bfctools, mpileup, samtools

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 06:57 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO