SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics
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

Reply
 
Thread Tools
Old 06-17-2015, 09:31 AM   #1
Greg_owens
Junior Member
 
Location: Vancouver, BC

Join Date: Aug 2014
Posts: 3
Default 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
Greg_owens is offline   Reply With Quote
Old 06-18-2015, 01:04 AM   #2
jmarshall
Samtools maintainer
 
Location: Cambridge, UK

Join Date: Jul 2009
Posts: 39
Default

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.
jmarshall is offline   Reply With Quote
Old 06-18-2015, 08:18 AM   #3
Greg_owens
Junior Member
 
Location: Vancouver, BC

Join Date: Aug 2014
Posts: 3
Default

That works. Thanks!
Greg_owens is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off



All times are GMT -8. The time now is 06:21 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO