SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools "is recognized as '*'" "truncated file" error axiom7 Bioinformatics 3 11-26-2014 02:53 AM
strategy to convert low GQ calls to "no call" manducasexta Bioinformatics 0 08-27-2012 08:49 AM
Specifying the region with "samtools mpileup -r " changes the output calls vplagnol Bioinformatics 0 06-07-2012 04:26 AM
"Bcftools mpileup" removes marked duplicates? angerusso Bioinformatics 0 03-05-2012 06:43 AM
comparing two samples using "bcftools -T" option? NanYu Bioinformatics 0 10-25-2011 07:01 AM

Reply
 
Thread Tools
Old 01-10-2014, 10:07 AM   #1
molecules
MU Informatics Research Core F
 
Location: Columbia, MO

Join Date: Oct 2010
Posts: 6
Default Samtools and bcftools: Why aren't some of the calls "missing data".

I'm using samtools mpileup and bcftools view to find variants. I'm perplexed why none of the alleles for any of the samples are shown as "missing data" (e.g. "./."). The coverage shouldn't be deep enough nor the samples similar enough for there not to be at least some "missing data". Here's my sample code:
Code:
samtools mpileup -uf ref.nfa A.bam B.bam C.bam  | bcftools view -vcg - > ABC.vcf
Is there a flag I'm missing or misusing?

This is using samtools version 0.1.19-44428cd and bcftools version 0.1.17-dev (r973:277) on CentOS 6.5 (64-bit).

Last edited by molecules; 01-10-2014 at 01:46 PM. Reason: Trying to fit code onto one line
molecules is offline   Reply With Quote
Reply

Tags
bam file, samtools bcftools, vcf file

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:06 PM.


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