SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Maximum read Depth in Samtools Anjali Bioinformatics 2 01-16-2012 02:15 AM
How to determine maximum read depth in Samtools VarFilter Muhua wang Bioinformatics 0 10-06-2011 03:09 PM
How to decide maximum read depth in Samtools VarFilter Muhua wang Illumina/Solexa 0 10-06-2011 03:02 PM
first read in a pair means F3 or R3? yhong General 1 08-24-2011 09:10 AM
About the read depth of coverage El Mariachi Illumina/Solexa 2 12-30-2010 12:22 AM

Reply
 
Thread Tools
Old 01-25-2011, 01:12 AM   #1
ketan_bnf
Member
 
Location: India

Join Date: Oct 2010
Posts: 59
Default What is read depth -D100 means in samtools?

Hi all,

I am working with exome reads, i extracted SNP and indels using,

# samtools mpileup -ugf ref.fa aln1.bam aln2.bam | bcftools view -bvcg - > var.raw.bcf
# bcftools view var.raw.bcf | vcfutils.pl varFilter -D100 > var.flt.vcf

i don't know what -D read depth means in vcfutils.pl?

i am having 101245 reads FYI.
ketan_bnf is offline   Reply With Quote
Old 07-28-2011, 06:54 PM   #2
Shambala
Junior Member
 
Location: hongkong

Join Date: Jul 2011
Posts: 1
Default

The -D option of varFilter controls the maximum read depth, which should be adjusted to about twice the average read depth.
Shambala 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 02:30 PM.


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