Biostar Beta. Not for public use.
ANI calculation: alternatives to Jspecies
0
Entering edit mode
12 months ago
fhsantanna • 440
Brazil

I am trying to calculate ANI values (ANIm and ANIb) using Jspecies in Ubuntu 14.04.

For the majority of comparisons, Jspecies is not able to attribute a ANI value (instead it gives "nan" result). Funny is that sometimes it works well.

Here is the log when errors occur:

[12:57]ERROR509909[Thread-11] - api.calculation.anim.RunANIm.calculateEMA(RunANIm.java:75) - .cluster file doesn't exist # seems to be related to nucmer  
[12:58]ERROR580395[Thread-11] - api.calculation.blast.FormatDB.runFormatDB(FormatDB.java:61) - [formatdb] FATAL ERROR: Database base name must be provided with multiple input files

Do you know an alternative software to jspecies?

Thanks.

jspecies ani • 4.8k views
ADD COMMENTlink
1
Entering edit mode

I discovered this http://jspecies.ribohost.com/jspeciesws/

Apparently they moved on the web. Maybe it can help you.

ADD REPLYlink
0
Entering edit mode

Great! Thank you very much. I am going to give it a try.

ADD REPLYlink
0
Entering edit mode

Let me understand, you submit different file in Jspecies and sometimes it works and sometimes not. Is that right?

ADD REPLYlink
0
Entering edit mode

Yes. I was not able to figure out what is the problem. Maybe it is a incompatibility with the new version of blast, because when I utilize an script for ani calculation, blast crashes sometimes...

ADD REPLYlink
0
Entering edit mode

I had a similar problem. Are you in Linux system? Did you import the files from windows? Which blast version are you using?

ADD REPLYlink
0
Entering edit mode

Yes, I am using Linux (Ubuntu 14.04). No, I didn' t import files from windows. Blastall version 2.2.26.

ADD REPLYlink
1
Entering edit mode

If I am not wrong ANIm uses mummer, does mummer work in your system? for the formatdb error, maybe it does not like the Ids of you seqs. You could give it a tray using one string, no space, only alphanumeric characters e.g. > MYSEQUENCE

ADD REPLYlink
5
Entering edit mode
4.2 years ago
aleimba • 130
Germany

An alternative that I can highly recommend is _pyani_ from Leighton Pritchard. No GUI just command-line, but that was one of the major drawbacks of the Java JSpecies anyway. Therefore, _pyani_ can handle lots of microbial genomes at once.

https://github.com/widdowquinn/pyani

Additionally, _pyani_ can also create heatmap diagrams corresponding to the ANI or Tetra values with a hierarchical clustering of the results (with additional QC graphics).

pyani relies heavily on some very common Python libraries and R or matplotlib for graphical output, but most bioinformaticians will have these installed anyway.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1