bcftools filter
1
0
Entering edit mode
6 weeks ago
Gregor Rot ▴ 490

Dear all,

I am doing some work on filtering variants, using bcftools for now:

bcftools view -i '%QUAL>=20' calls.bcf

Why the % at %QUAL? Trying to figure out looking at manuals for bcftools expressions:

https://samtools.github.io/bcftools/howtos/variant-calling.html
https://samtools.github.io/bcftools/howtos/filtering.html

Thanks & Cheers,
Gregor

filter bcftools • 142 views
ADD COMMENT
1
Entering edit mode
6 weeks ago

I thinK '%ID' and '%QUAL' are used to make them distinct from INFO/ID and INFO/QUAL

and in fact, you can use both: https://github.com/samtools/bcftools/blob/682949a368ea0d4844b7b1791caf7f3df91a9a50/filter.c#L2641

        else if ( !strncasecmp(str,"QUAL",len) || !strncmp(str,"%QUAL",len) /* for backward compatibility */ )
ADD COMMENT
0
Entering edit mode

Thanks, this is great. Just to be sure: when you write INFO/ID and INFO/QUAL, you mean that in the 8th field (INFO field) where we can have several values separated by ; we could have a value with name INFO or ID? (distinct from the 6th field QUAL and 3rd ID). Looking at VCF ref here: https://samtools.github.io/hts-specs/VCFv4.2.pdf

ADD REPLY

Login before adding your answer.

Traffic: 1889 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