I just finish a bacteria genome, then compare two genomes of the same genus. I hesitate between the MUMmer and blastp. Please give me some suggestions. And also, is there any software to get a dot plot from the result of blastp?
-
-
yes, you can try BSR, it is based on blastp.
good luck! -
Thank you for bio-x's kind help. I think accurately I'd like to get a dot plot of two genomes's ortholog plot by reciprocal best-hit pair from the result of pair-wise BLASTP. I guess there must be some software which is suitable to do this. Please give me some suggestions.Comment
-
"promer" from the Mummer suite will effectively do want you want. Although it doesn't examine "genes" it matches translated regions, which if there is a good match, are likely to be orthologs (or paralogs) anyway.Originally posted by anyone1985 View PostComment
-
Yes, promer is pretty good. However, I'd like to draw the graph exactly. The promer can not give a clear position to identify the gene. Many literatures have figures as below. I am intrested in which software can draw such figure.
Originally posted by Torst View PostComment
-
does anyone know how to fix the following MUMmer error?
where textlen is the same as the length of build36 (hg18)
./mummer: suffix tree construction failed: textlen=3080436075 larger than maximal textlen=536870908Comment
-
Try:
http://seqanswers.com/wiki/Talk:MUMmer
The corresponding paper can be found here:
http://bioinformatics.oxfordjournals...ull/25/13/1609Comment
-
@robs
thanks, but the link says:
[Thanks for the hint about the recent "MEMs" implementation that "can serve as a drop-in replacement for the MEMs algorithm in MUMmer 3". Seems to be working perfectly where MUMmer crashed :-) --Dan 15:08, 5 January 2010 (UTC)]
what was the original hint? (sorry for any stupidity)Comment
-
Not sure what you mean with "original hint". There is only one hint and if you take a look at the paper you will see that the quote is actually from there.
The MEMs implementation has a mummer module that can be used instead of the MUMmer3 one. By doing this, mummer shouldn't crash anymore. More details are in the paper.Comment
-
@robs
thanks so much for your help.Comment
-
To solve the textlen problem, build a 64-bit version, and make sure
#define SIXTYFOURBITS is set in types.h.
(Grep for the one place this is used, to see where to set it).Comment
-
First, thanks for, as far as I can tell, the only solution to the mummer textlen problem on the web...Originally posted by Will Nelson View Post
Not being an hardcore programmer, but a biologist-come-bioinformatician (more like a user of programs), it took me some time to figure out how to implement it, though. What you meant was to add this line
to the src/kurtz/libbasedir/types.h file before the line where it saysCode:#define SIXTYFOURBITS
and then runCode:#ifdef SIXTYFOURBITS
again.Code:make install
(I had begun to 'uncomment' the ifdef line, causing all kinds of errors - like I said, not a programmer)
Just clarifying for future generations...Comment
-
This thread also supplies another solution:
"
If you recompile the package like so:
> make clean
> make CPPFLAGS="-O3 -DSIXTYFOURBITS"
you should have no more problems with space limitations.
Best,
-Adam
"Comment
-
I want to compare to genomes with 87%similarity, I cant decide should I use Promer or NUCMER?Comment
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...09-01-2026, 05:41 AM
Comment