SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Question on calling SNPs using samtools/bcftools nkwuji Bioinformatics 6 02-19-2013 08:52 AM
Read Depth in vcf (samtools / bcftools) Marie_Noir Bioinformatics 1 04-17-2012 06:48 AM
samtools/bcftools missing obvious SNPs? SamH Bioinformatics 3 01-17-2012 02:16 PM
Conflict between mpileup/bcftools and GATK in VCF file ericarcher Bioinformatics 0 09-25-2011 04:33 PM
Predicting true SNPs from .vcf file swbarnes2 Bioinformatics 1 04-06-2011 03:29 PM

Reply
 
Thread Tools
Old 09-05-2012, 06:53 AM   #1
muenalan
Junior Member
 
Location: Hannover

Join Date: Sep 2012
Posts: 3
Default samtools/bcftools: How to inspect SNPs not made it into vcf file

Hi,

Is there any parameter to "log" the snp-calling for individual SNP's which did not make it into the .vcf ? Or any other way to look how are the parameters for that SNP, if called ?

Info: We generate our vcf's like this:

samtools mpileup -ugf 'human_build37_ucsc.fa' 'test.bam' | bcftools view -bvcg - > 'test.bam.bcf'
bcftools view 'test.bam.bcf' > 'test.bam.bcf.vcf'


Thx,
muenalan
muenalan is offline   Reply With Quote
Old 09-06-2012, 01:46 AM   #2
AWeller
Junior Member
 
Location: Germany

Join Date: Feb 2011
Posts: 4
Default

I'm not sure if I understand your question correctly.

In general, omitting the "-v" in bcftools will show you all positions, not just the variants.

A convenient way to manually inspect a single position that didn't get called is the "samtools mpileup -r" option, where eg "samtools mpileup -r CHROMOSOMEII:10-10 ..." will show you only this one location.

If you are talking about doing this in bulk for the whole genome, a possibility is to do a second pileup of all positions and then find the differences via "diff" or a custom script.

Hope to help,

best
Andreas
AWeller is offline   Reply With Quote
Old 09-06-2012, 05:13 AM   #3
muenalan
Junior Member
 
Location: Hannover

Join Date: Sep 2012
Posts: 3
Default [solved]

Yes,

Exactly the information I was looking for.

Thx,
M!
muenalan 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 10:23 PM.


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