Go Back   SEQanswers > Applications Forums > Genomic Resequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools mpileup calls way too less SNPs TuA Bioinformatics 17 03-01-2018 05:17 PM
SAMTOOLS mpileup only INDELs, no SNPs? jkozubek Bioinformatics 3 09-10-2012 10:16 AM
samtools/mpileup heterozygous SNPs calling combiochem Bioinformatics 4 08-02-2011 08:05 AM
calling Heterozygous SNPs with samtools mpileup egatti Bioinformatics 1 07-21-2011 09:16 AM
samtools mpileup filter SNPs Hit Bioinformatics 3 05-25-2011 05:55 PM

Thread Tools
Old 11-29-2012, 07:37 PM   #1
Location: USA

Join Date: Oct 2010
Posts: 38
Default Problems when calling SNPs using samtools mpileup

Hi there,
Recently, I'm using samtools mpileup to call SNPs using low coverage (average 1x) resequencing data of 600+ populations. But I encountered a few problems.

My command line is as follows:
samtools mpileup -E -C50 -Q20 -q20 -DSuf ref.fasta -r chr4:21000001-24000000 -b bam_list | bcftools view -Nbcvg - > partE-chr4:
The individuals sequenced are of inbred lines. So we want to get homoSNPs. HeteroSNPs shold be wrong calling.

1) The SNPs printed by samtools mpileup seemed weird:
SNP1: 0/1:26,3,0:1:0:10
SNP2: 0/1:26,3,0:1:0:9
Both SNP1 and SNP2 have PLs with (26,3,0), but samtools- bcftools printed them as hetero. lh3 have give a answer at the post:
but finally which genotype should I give? First, I think I should compare the three values of PLs, and choose the maximum value as
the final genotype. If two of the values were very close, how to give the genotype? I'm not sure of this. Can someone give me some
suggestions and some explanations?

2) For large population, how to deal with QUAL columns (the sixth column)? How to choose the cutoff?


Last edited by zeam; 11-29-2012 at 07:44 PM. Reason: Add some information
zeam is offline   Reply With Quote

heterosnp, quality control, samtools mpileup

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 11:58 PM.

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