SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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
Impact on quality of SNP calls using samtools mpileup myi Bioinformatics 3 03-03-2014 08:18 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

Reply
 
Thread Tools
Old 01-20-2012, 02:20 PM   #1
agel
Junior Member
 
Location: USA

Join Date: Dec 2011
Posts: 5
Default Samtools mpileup calls drastically more SNPs with -I

The command:
Code:
samtools mpileup -b bams.txt -uDSf human_b36_male.fa | bcftools view -vcg - > var.raw.vcf
where bam.txt is a list of BAM files, outputs a var.raw.vcf file with about 200 lines, and the command finishes almost immediately. The same command with the -I flag has so far outputted 3000 lines and is still going. -I is supposed to suppress indel calling, why is it affecting the number of SNPs called so drastically?
agel is offline   Reply With Quote
Reply

Tags
indel, mpileup, samtools

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


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