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!
1 answer
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.
Log in to answer this question.
-
Only output SNPs in bcftools
written by NewtoBioinfo •
-
filtering on F_PASS(GT="het") in bcftools gives contradictory results
written by scottbrainard •
-
direct call of indels/snps with bcftools?
written by cristina_sabiers
-
Empty VCF file with bcftools call
written by AP
-
genotype calling using samtools
written by Karim •
-
bcftools "The argument to -g not recognised"
written by AP
-
How to Convert SNP's that produced from Varscan into 0,1,2 ?
written by bingnas •
-
How to Convert SNP's that produced from Varscan into 0,1,2 ?
written by bingnas •
-
Bcftools view options for SNPs calling.
written by Eddie_I •
-
Missing V600E Mutation In Wholeexome Analysis Of A375 Cell Line? Samtools Problem?
written by Angel
