confusion about BCFtools "EXPRESSIONS" flags
0
0
Entering edit mode
7.8 years ago
lstbl ▴ 40

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 https://samtools.github.io/bcftools/...ml#expressions

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'.

samtools bfctools mpileup • 1.9k views
ADD COMMENT

Login before adding your answer.

Traffic: 2580 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6