Difference between bcftools view and bcftools call?
Entering edit mode
6.8 years ago
AP ▴ 100

Hi all,

This post is related to a previous one I wrote yesterday. Would anyone be able to tell me what is the difference between bcftools view and bcftools call?

My understanding is that bcftools call is the new view for bcftools versions 1.2 and above (I am using 1.3). However, it seems that filter options from view are not compatible with call (while I thought some would).

For instance, I am trying to use the -g filter options from bcftools view with call by typing the following but without success:

bcftools call -v -c -g hom|het|miss output.bcf > output.vcf

Alternatively, doing the following does not work either:

bcftools view -v -c -g hom|het|miss output.bcf > output.vcf

Is it possible to use the -g filter option from view with call? Any thoughts would be much appreciated!

bcftools view call • 3.8k views
Entering edit mode
5.7 years ago
AP ▴ 100

13 months later:

Using bcftools view should work although view was replaced by call in recent versions of bcftools (1.2 and 1.3). The -g filter option triggers an error probably due to a conflict between versions. Using bcftools call just does the job. Otherwise you need to go back to older versions.


Login before adding your answer.

Traffic: 2624 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6