SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
vcfutils.pl -specifying minimum read depth? Lspoor Bioinformatics 3 05-27-2013 01:20 AM
GATK excludes some samples for cohort variant calling liu_xt005 Bioinformatics 2 02-01-2012 11:58 AM
Samtools variant calling questions Chiel Bioinformatics 2 06-07-2011 09:10 AM
variant calling using samtools -v- bcftools ksc Bioinformatics 2 04-13-2011 06:44 AM
What's maximum depth and minimum depth? slny Bioinformatics 0 03-17-2011 10:13 AM

Reply
 
Thread Tools
Old 12-06-2011, 08:01 AM   #1
m_elena_bioinfo
Member
 
Location: Ospedali Riuniti di Bergamo, ITALY

Join Date: Oct 2009
Posts: 99
Default minimum depth variant calling samtools/gatk

Hi NGS user,
does anyone know how I can set a minimum coverage threshold in variant calling in samtools mpileup and in GATK Unified Genotyper. I only find parameter to set max per-BAM depth to avoid excessive memory usage.

Thanks a lot,
ME
m_elena_bioinfo is offline   Reply With Quote
Old 12-06-2011, 08:31 AM   #2
Heisman
Senior Member
 
Location: St. Louis

Join Date: Dec 2010
Posts: 535
Default

With SAMtools you can filter out calls under X coverage using:

/samtools-0.1.18/bcftools/vcfutils.pl varFilter -d X <input.vcf> > <output.vcf>
Heisman 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:07 AM.


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