SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Mira assembly jmpi Bioinformatics 12 02-16-2015 08:47 PM
MIRA assembly with MID barcodes for 454? raw937 Bioinformatics 4 06-15-2011 12:54 AM
MIRA assembly with MID barcodes for 454? raw937 454 Pyrosequencing 1 06-13-2011 12:52 PM
Discussion about MIRA hybrid assembly of 454 reads with Illumina unpaired data edge De novo discovery 5 11-16-2009 02:17 AM
A qusetion about denovo assembly 454 sequence using MIRA kentnf Bioinformatics 7 04-24-2009 06:36 AM

Reply
 
Thread Tools
Old 12-20-2010, 09:34 AM   #1
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,317
Question Extract scaffold info from Mira PE 454 assembly?

Does Mira create scaffolds during paired-end read sequence assemblies? gsAssembler does, so I just assumed Mira would. But maybe I need to run downstream software (Bambus?) to extract this information?

--
Phillip
pmiguel is offline   Reply With Quote
Old 12-20-2010, 11:56 AM   #2
sklages
Senior Member
 
Location: Berlin, DE

Join Date: May 2008
Posts: 628
Default

No, MIRA doesn't create scaffolds; you should use an external tool like BAMBUS.

cheers,
Sven
sklages is offline   Reply With Quote
Old 12-21-2010, 06:19 AM   #3
themerlin
Member
 
Location: Flagstaff, AZ

Join Date: Feb 2010
Posts: 51
Default

You can use the MIRA output for scaffolding with Bambus (see tutorial)

http://mira-assembler.sourceforge.ne...IRA_BAMBUS.pdf

You can also use the MIRA output with SSPACE with only minimal manipulation:

http://www.baseclear.com/sequencing/...-tools/sspace/

J
themerlin is offline   Reply With Quote
Old 12-21-2010, 12:24 PM   #4
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,317
Question Bambus Obsolete?

I am getting the mysterious "script failed" error mentioned in this thread:

http://seqanswers.com/forums/showthread.php?p=21874

Therein it is posited that there is some maximum number of contigs that can be processed. Anyone know what that maximum is or have other insights?

Also, given that no modifications to Bambus have been checked in to its SourceForge site since 2005 -- isn't this really a pre-next gen program?

--
Phillip
pmiguel is offline   Reply With Quote
Old 12-22-2010, 04:53 AM   #5
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,317
Default Bambus works for gsAssembler .ace file

On the other hand, for smaller 454 jobs it does work for Mira assemblies. And with a few modifications of the normal protocol, works directly on a gsAssembler .ace file.

Not sure if gsAssembler creates a .ace file by default. If not, use the appropriate option to have one created. Then create the .contig file from the .ace file as normal using the amostools utility:

ace2contig -i 454Contigs.ace 454Contigs.contig

The ace file created by gsAssembler uses a different paradigm to denote forward and reverse reads. So I used the following 454Contigs.mates file:

Code:
library all     1000    6000    (.......).*
pair    (.*)_left       (.*)_right
As always for bambus, the field delimiter (shown above as whitespace) for a .mates must be tabs. So a copy/paste from this post will not work without modification of the resultant file.

Code:
perl -i -pe 's/  */\t/g' 454Contigs.mates
should do the trick. Note that is two spaces prior to the "*", not one. If you only put one, that would allow perl to place a tab in between every character.

One other catch for me. I had been skipping the .conf file, allowing goBambus to substitute its own default one. But I got the following this time:

Quote:
Step 400: Running scaffolder
Error: Priority not specified: at least one library must be assigned a priority in module grommit
Grommit(/usr/local/bambus-2.33/bin/grommit -i 454Contigs_bambus.inp -o 454Contigs_bambus.out.xml -C 454Contigs_bambus.grommit.conf --append --logfile goBambus.log --debug 1) script failed
So I created a .conf file, 454Contigs_bambus.conf with the following in it:

Code:
# Priorities
priority ALL 1

# Redundancies
redundancy 2

# min group size
mingroupsize 0

Then the following command:

Code:
goBambus -c 454Contigs.contig -m 454Contigs.mates -o 454Contigs_bambus -C 454Contigs_bambus.conf
appears to work.

Of course, gsAssembler should create scaffolds itself. But, for this project it did a very poor job, creating a few tiny scaffolds. So it is nice to have Bambus as a backup.

--
Phillip
pmiguel 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:13 PM.


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