SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SAMtools missing bcftools khadeeja Bioinformatics 8 06-19-2013 02:06 AM
Missing SNPs in samtools output eventhough it display VDB values k2bhide Bioinformatics 1 03-04-2013 12:54 AM
samtools variant calling for separate chromosomes tahamasoodi Bioinformatics 9 12-04-2012 10:45 AM
Samtools on Chromosomes Versus Genome ramirob Bioinformatics 9 05-22-2012 08:19 AM
samtools/bcftools missing obvious SNPs? SamH Bioinformatics 3 01-17-2012 02:16 PM

Reply
 
Thread Tools
Old 03-12-2013, 09:50 AM   #1
k2bhide
Junior Member
 
Location: West Lafayette

Join Date: Nov 2010
Posts: 5
Default Samtools - Missing SNPs in some chromosomes

Hello all,
I called SNPs using samtools(Version 0.1.18) on haploid genome and I found some interesting things. No SNPs were called for some chromosomes when SNP calling is carried out for entire genome however when the same parameters are used to call SNPs for those chromosomes individually, many SNPs are found.
Could you please let me know why would there be difference in number of SNPs when called for entire genome and individual chromosomes?

It would be really helpful if anyone of you could clarify how are SNPs called by samtools and explain why would there be difference (for SNPs called on entire genome and SNPs called on individual chromosomes). I have pasted samtools commands that I used:

1) Samtools command used to call SNPs from entire haploid genome
samtools mpileup -ABuf ref.fasta sample.bam | bcftools view -cvg -> sample.var.raw.vcf
No SNPs were called for chromosome1 using above command line.


2) Samtools command used to call SNPs for only chromosome1
samtools mpileup r chr1 -ABuf ref.fasta sample.bam | bcftools view -cvg -> sample_chr1.var.raw.vcf
Many SNPs were called for chromosome 1 using above command line.

Let me know your comments and suggestions.
k2bhide 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 09:02 PM.


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