SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Getting FQ from bcftools 1.2 (http://seqanswers.com/forums/showthread.php?t=60440)

Greg_owens 06-17-2015 09:31 AM

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

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!


All times are GMT -8. The time now is 07:04 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.