Getting FQ from bcftools 1.2

Greg_owens · 06-17-2015, 09:31 AM

I'm using samtools and bcftools to call hapmap snps from WGS data. With the older version of samtools the vcf file included FQ per site, but when I use the latest version (samtools/bcftools 1.2) it doesn't. I need FQ for a downstream program, so is there a way to force bcftools to output it?

Here's what my code looks like now:

Code:

samtools mpileup -ugf hg19.fasta -l annotHapMap2L.txt $bam/$name.bam | bcftools call -mO v -o $vcf/$name.HM.vcf -V indels -P 0.5

jmarshall · 06-18-2015, 01:04 AM

Current bcftools 1.x contains two calling methods: call -c is the original bcftools consensus caller, while call -m is a more recent multi-allelic caller. Only the consensus caller outputs FQ.

If your downstream program is expecting to get old-bcftools-style input, you probably want to be using call -c.

Greg_owens · 06-18-2015, 08:18 AM

That works. Thanks!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
bcftools issue	willgrand	Bioinformatics	11	12-03-2015 01:20 PM
Bcftools stats	Nnogueira	General	0	04-17-2015 04:08 AM
SNV ignored by bcftools	behrj	Genomic Resequencing	0	11-14-2014 05:44 AM
samtools bcftools no FQ value	StefanP	Bioinformatics	2	06-20-2014 02:51 PM
bcftools error	Morpse	Bioinformatics	8	10-05-2012 11:24 AM

06-17-2015, 09:31 AM	#1
Greg_owens Junior Member Location: Vancouver, BC Join Date: Aug 2014 Posts: 3	Getting FQ from bcftools 1.2 I'm using samtools and bcftools to call hapmap snps from WGS data. With the older version of samtools the vcf file included FQ per site, but when I use the latest version (samtools/bcftools 1.2) it doesn't. I need FQ for a downstream program, so is there a way to force bcftools to output it? Here's what my code looks like now: Code: samtools mpileup -ugf hg19.fasta -l annotHapMap2L.txt $bam/$name.bam \| bcftools call -mO v -o $vcf/$name.HM.vcf -V indels -P 0.5

06-18-2015, 01:04 AM	#2
jmarshall Samtools maintainer Location: Cambridge, UK Join Date: Jul 2009 Posts: 39	Current bcftools 1.x contains two calling methods: call -c is the original bcftools consensus caller, while call -m is a more recent multi-allelic caller. Only the consensus caller outputs FQ. If your downstream program is expecting to get old-bcftools-style input, you probably want to be using call -c.

06-18-2015, 08:18 AM	#3
Greg_owens Junior Member Location: Vancouver, BC Join Date: Aug 2014 Posts: 3	That works. Thanks!