SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BEDTools coverageBed error JMorlan Bioinformatics 1 04-22-2018 11:17 PM
bedTools 2.20.1 coverageBed cmccabe Bioinformatics 1 06-18-2015 07:10 AM
Using coverageBed in bedtools DrCheymus Bioinformatics 11 11-01-2013 04:09 AM
differences in coverage BEDtools - samtools mpileup KNS Bioinformatics 3 12-01-2011 09:55 PM
BEDTools..coverageBed query Siva Bioinformatics 2 08-18-2010 11:10 AM

Reply
 
Thread Tools
Old 09-16-2015, 03:09 AM   #1
akvarel
Member
 
Location: Berlin

Join Date: Nov 2011
Posts: 11
Default bedtools coverageBed vs samtools mpileup, speed

I would like to use bedtools as it seems to have much more options for the input/output file. However, it seems that bedtools is way slower than the samtools.


I have: multiple bad files (sorted and indexed) and multiple regions and I want a coverage per base.
gene_coord_red.bed has only one region.
I ran bedtools coverageBed(v2.25.0) following:


Code:
time coverageBed -a gene_coord_red.bed -b reads.sort.bam -d > bedtools.txt
real    2m17.861s
user    1m57.572s
sys    0m19.926s
I ran samtools mpileup:

Code:
time samtools mpileup -Q 0 -l gene_coord_red.bed reads.sort.bam > mpileup.txt

real    0m23.595s
user    0m22.808s
sys    0m0.301s
​Why is the bedtools running so slow? Can I accelerate it somehow?
akvarel is offline   Reply With Quote
Old 09-16-2015, 03:19 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

See also on biostars.
dpryan 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 05:47 AM.


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