SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Visualizing Mauve Alignments sphil Bioinformatics 6 04-22-2014 11:24 PM
PubMed: Mauve assembly metrics. Newsbot! Literature Watch 0 01-04-2012 03:10 AM
Are there any alignment programs that take SNPs into account during alignment? sdarko Bioinformatics 2 06-04-2011 06:09 AM
PubMed: Exact and complete short read alignment to microbial genomes using GPU progra Newsbot! Literature Watch 0 04-01-2011 03:00 AM
.bam alignment format from 1000 Genomes Project CellsDividing Bioinformatics 2 02-06-2009 09:21 AM

Reply
 
Thread Tools
Old 03-31-2012, 10:58 AM   #1
gmer
Junior Member
 
Location: UK

Join Date: Mar 2012
Posts: 3
Default Alignment of >30 genomes with mauve?

Hi all,

Relatively new to all this, I tried progressiveMauve to align more than >25 bacterial genomes and it repeatedly crashed after a few hours. When I align ~10 it works perfectly. I was told by other members of my lab that mauve tends to not work so much with a larger number of genomes.

The question may seem trivial to most of you, sorry in advance, but which genome alignment tool would you use to do this? Still, I like mauve; so would there be a way for it to work with a larger number of genomes?

Thanks for your time and help!

Best,

G
gmer is offline   Reply With Quote
Old 03-31-2012, 11:33 AM   #2
bckirkup
Member
 
Location: Washington DC

Join Date: Jan 2011
Posts: 17
Default Whole Bacterial Genome Alignment

Of course, Mummer is the classic 'alignment' tool; and Progressive Mauve is built to purpose; but I may have something which can help you - a colleague has built a tool that uses progressive ortholog identification to construct phylogenies; but by virtue of the method, it can provide the core information for a progressive alignment as well - and currently it is robust to a relatively large number of genomes. It would need some extension to provide assembled whole genome or contig alignment.
bckirkup is offline   Reply With Quote
Old 03-31-2012, 11:52 AM   #3
nickloman
Senior Member
 
Location: Birmingham, UK

Join Date: Jul 2009
Posts: 356
Default

I've successfully aligned 20 genomes with Mauve but it took some time, over a day and I needed to assign extra memory to the task (do this via the Java options a la http://www.mail-archive.com/mauve-us.../msg00137.html)

You might also want to consider looking at Mugsy from the Salzberg group (http://bioinformatics.oxfordjournals...tq665.abstract) which can apparently align 31 Streptococcus genomes in 2 hours.

And of course bckirkup's suggestion sounds interesting too!
nickloman is offline   Reply With Quote
Old 04-01-2012, 07:01 AM   #4
gmer
Junior Member
 
Location: UK

Join Date: Mar 2012
Posts: 3
Default

I will try Mugsy and come back here if any other problem. Thanks!

@bckirkup: what's the name of the tool made by your colleague?

G
gmer is offline   Reply With Quote
Old 04-07-2012, 08:57 PM   #5
srividya22
Junior Member
 
Location: New Jersey

Join Date: Mar 2012
Posts: 4
Default Reg Mugsy Usage

Hi,

I am trying to align human chr1 and mouse chr 1 using Mugsy. I am facing difficulties in getting it to work.

These are the steps I followed. Can someone help.

I have the 2 files Ch1_human, chr1_mouse in FASTA format.

I could install Mugsy
It asks me to set the the MUGSY_INSTALL path. I set it using
bash
export MUGSY_INSTALL= <working directory>

source <path of mugsyenv.sh>

I dont have root permissions on the server I use.

I get the following error in mugsy.
20120407|225222| 4601| ERROR: The following critical files could not be used
20120407|225222| 4601| ./mugsy_output/tmp.chr1.delta
20120407|225222| 4601| ./mugsy_output/tmp.chr1.cluster
20120407|225222| 4601| ./mugsy_output/tmp.chr1.mgaps
20120407|225222| 4601| ./mugsy_output/tmp.chr1.ntref
20120407|225222| 4601| /afs/cad.njit.edu/research/ccs/usman/12/srividya/Thesis/mugsy_x86-64-v1r2.2/mugsy_output/mugsy_output/tmp.chr1.queries.fsa
20120407|225222| 4601| /afs/cad.njit.edu/research/ccs/usman/12/srividya/Thesis/mugsy_x86-64-v1r2.2/mugsy_output/mugsy_output/chr1
20120407|225222| 4601| Check your paths and file permissions and try again

I checked the file permissions. it is proper.

can anyone please help ?
srividya22 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 06:11 PM.


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