Entering edit mode
2.6 years ago
Filago
▴
100
Hello, I have a multi-sample VCF-file and want to filter out variants in single patients based on the ratio of certain quality values (which are sample specific and annotated in the FORMAT-field)
So I filtered my data like this:
bcftools view file.vcf.gz -i '(FORMAT/quality_value1)/(FORMAT/quality_value2) > 10' -o file_flt.vcf.gz -Oz
However, this step filtered whole variants instead of single sample variant calls. Any ideas about this issue?
Best, Andreas