SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Has anyone used the Genovo assembler? aroth8 De novo discovery 10 07-30-2012 11:23 AM
k-range in Meta-IDBA loba17 Metagenomics 3 04-05-2012 01:09 PM
Velvet assembler bioinf Bioinformatics 31 08-24-2011 10:19 AM
PE Assembler ewalt98 Bioinformatics 2 04-01-2011 02:21 PM
SOAPdenovo assembler shailesh Bioinformatics 6 10-04-2010 07:35 PM

Reply
 
Thread Tools
Old 11-11-2012, 01:10 PM   #1
Mark
Member
 
Location: Raleigh, NC

Join Date: Nov 2008
Posts: 48
Default IDBA-UD assembler

Hi All

Has anyone given the idba-ud assembler a try? My first go at it produced a bacterial genome assembly quite inferior to what I had achieved with Velvetoptimiser, but proper parameterization of idba is not exactly clear to me. What has your experience been?

Mark
Mark is offline   Reply With Quote
Old 02-11-2013, 10:05 AM   #2
lanzen
Junior Member
 
Location: Bilbao, Basque Country

Join Date: Mar 2012
Posts: 2
Default

Hi Mark,

I am about to try it now with some mixed (by mistake) metagenome and metatranscriptome data. I found the algorithm very appealing for this purpose, since the main problem will be extremely variable coverage. I will let you know how it turns out.

Did you try it for metagenome data? I don't think there is much point in using for a normal single genome assembly..
lanzen is offline   Reply With Quote
Old 03-11-2013, 09:16 AM   #3
luperci
Junior Member
 
Location: Nashville, TN

Join Date: Mar 2009
Posts: 3
Default

lanzen,

How did it go? I'm currently trying to use idba-ud to do a metagenomic assembly of a goat gut. Do you know if there is a way to limit the memory usage? We run all our analysis using a cluster, and my job got killed because it went over the requested memory.
luperci is offline   Reply With Quote
Old 03-11-2013, 09:22 AM   #4
lanzen
Junior Member
 
Location: Bilbao, Basque Country

Join Date: Mar 2012
Posts: 2
Default

For me it worked quite fine, but most contigs were short, so it is unclear to me what was the (contaminating) metagenomic data that I tried to remove or what was from shorter transcripts. This may not be the fault of the assembly algorithm though. I will try a more normal dataset and let you know.

Otherwise, it did not use that much memory or calculation time. A 32Gb machine seemed with 8 cores was enough and the job finished in 30 min.

Let me know if there is something else I can do to help you.
lanzen is offline   Reply With Quote
Old 03-11-2013, 09:29 AM   #5
luperci
Junior Member
 
Location: Nashville, TN

Join Date: Mar 2009
Posts: 3
Default

I allocated 30 cores and 300 GB and it ran for around 7 hours before being stopped. I have a feeling the sample was really over sequenced. Right now I looking into doing some preprocessing to clean up the reads a little. Hopefully that will help limit the resource usage. Thanks.
luperci is offline   Reply With Quote
Old 05-29-2013, 06:28 AM   #6
mstagliamonte
Member
 
Location: Florida

Join Date: Feb 2013
Posts: 33
Default

Hi,


sorry to revive an old thread. I did two idba runs, one using only one fasta file, the other using two files (two separate illumina runs for the same sample).

In the first idba run, the file contig.fa had no read counts, while they were present in the scaffold.fa.
At the second attempt (using two fasta files), it was the opposite!

Any idea what's happening?
Max
mstagliamonte is offline   Reply With Quote
Reply

Tags
assembly, idba.velvet

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


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