This is a test version of Biostars. For the public version, visit https://www.biostars.org.
quick way to get samples that are NOT homozygous reference for a given record using bcftools?

I can do this in python, but just wondered if there is any way in bcftools to return a vcf with samples that are NOT homozygous reference for a given snp. I am guessing no, but thought I would ask anyways

bcftools
• 1,930 views • link
written by curious curious

1 answer 

 bcftools view -i 'count(GT=="RR")==0'
• 0 views • link
updated by GenoMax GenoMax • written by Pierre Lindenbaum Pierre Lindenbaum

[filter.c:1379 filters_init1] Error: the tag "INFO/count" is not defined in the VCF header ?

• 0 views • link
written by curious curious

what is your version of bcftools ? check you're using a recent version.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

that was it thanks

• 0 views • link
written by curious curious

Log in to answer this question.

More posts like this