SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
[Velvet,assembly] core dumped occured by runnning velvet matador3000 De novo discovery 0 12-17-2011 08:31 AM
Velvet assembly katussa10 De novo discovery 2 12-12-2010 11:22 AM
Scaffolding output from a velvet assembly of SOLiD data xavier De novo discovery 6 02-24-2010 04:00 AM
Scaffolding problem with Velvet Melissa Bioinformatics 2 01-05-2010 07:44 PM
velvet assembly vani s kulkarni Bioinformatics 25 10-07-2009 12:27 PM

Reply
 
Thread Tools
Old 03-28-2011, 02:30 PM   #1
Canadian_philosophy
Junior Member
 
Location: Canada

Join Date: Apr 2010
Posts: 9
Question Scaffolding Translational mapping (STM) - velvet assembly summary?

Hello,

Has anyone used STM [http://genome.cshlp.org/content/early/2010/08/30/gr.103846.109.abstract[/URL] for velvet assembly summaries?

I am trying to use it and have issues (St9bad_alloc) with multiple-k method - cd-est. Has anyone seen this before?

Thanks!
Canadian_philosophy is offline   Reply With Quote
Old 05-10-2011, 08:58 PM   #2
wreid
Junior Member
 
Location: Auburn, Alabama

Join Date: Dec 2010
Posts: 1
Default

Hi Canadian.

We have tried this method and have had good results with the multiple-K (still not able to get the actual STM part to work, may be a problem with biopython?). For cd-hit-est, you have to install it separately. If you google "cd-hit-v4.3-2010-10-25.tgz" it should pull up the address for the tarball download. One thing that we also did was to run Velvet-Oases prior to the multiple K. I think the Oases step is better to include if you're doing a transcriptome (I'm assuming this if you're interested in the STM method). Are you going to also try the STM method?

Cheers!
wreid is offline   Reply With Quote
Old 05-12-2011, 01:24 PM   #3
Canadian_philosophy
Junior Member
 
Location: Canada

Join Date: Apr 2010
Posts: 9
Thumbs down

Hi there,

Thanks for your reply. You are right, I am trying to use STM for transcriptome assembly. I have assemblies from Velvet/oases and I am trying to use that to run the cd-hit-est and multiple-k method. My files are like this :
/home/deepthi/All_kmers_fasta --> has contigs_25.fa contigs_29.fa contigs_31.fa Lane5_velvet_assembled_27.fa Lane678_velvet_seq_27.fa #contigs obtained from different k mers from velvet/oases#
My command-line is:
/home/deepthi/cd-hit-est -i /home/deepthi/All_kmers_fasta -o /home/deepthi/All_kmers_fasta_clus -c 0.9 -n 5 -M 800

Error looks like this :
\"Apprixmated maximum memory consumption: 593M
writing new database

Fatal Error:
file opening failed
Program halted !!

Traceback (most recent call last):
File \"/usr/local/software/STM/multiple-K.py\", line 84, in
<module>
os.remove(\'%s.bak.clstr\'%outfile)
OSError: [Errno 2] No such file or directory:
\'/home/deepthi/All_kmers_fasta.bak.clstr\'


Multiple-k method with this line doesn't seem to work either : sudo /usr/local/software/STM/multiple-K.py -d /home/deepthi/All_kmers_fasta -o /home/deepthi/All_kmers_fasta -v

The first part using the reference genome worked -- was that your question?
Canadian_philosophy 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 08:46 PM.


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