BioJava
BioJava is an open-source software project dedicated to provide Java tools to process biological data.[1][2][3] BioJava is a set of library functions written in the programming language Java for manipulating sequences, protein structures, file parsers, Common Object Request Broker Architecture (CORBA) interoperability, Distributed Annotation System (DAS), access to AceDB, dynamic programming, and simple statistical routines. BioJava supports a huge range of data, starting from DNA and protein sequences to the level of 3D protein structures. The BioJava libraries are useful for automating many daily and mundane bioinformatics tasks such as to parsing a Protein Data Bank (PDB) file, interacting with Jmol and many more.[4] This application programming interface (API) provides various file parsers, data models and algorithms to facilitate working with the standard data formats and enables rapid application development and analysis.
![]() | |
Original author(s) | Andreas Prlić |
---|---|
Developer(s) | Amr ALHOSSARY, Andreas Prlic, Dmytro Guzenko, Hannes Brandstätter-Müller, Jose Manuel Duarte, Thomas Down, Michael L Heuer, Peter Troshin, JianJiong Gao, Aleix Lafita, Peter Rose, Spencer Bliven |
Initial release | 2002 |
Stable release | 6.0.3
/ December 19, 2021 |
Repository | github |
Written in | Java |
Platform | Web browser with Java SE |
Available in | English |
Type | Bioinformatics |
License | Lesser GPL 2.1 |
Website | biojava |
Additional projects from BioJava include rcsb-sequenceviewer, biojava-http, biojava-spark, and rcsb-viewers.
Features
BioJava provides software modules for many of the typical tasks of bioinformatics programming. These include:
- Accessing nucleotide and peptide sequence data from local and remote databases
- Transforming formats of database/ file records
- Protein structure parsing and manipulation
- Manipulating individual sequences
- Searching for similar sequences
- Creating and manipulating sequence alignments
History and publications
The BioJava project grew out of work by Thomas Down and Matthew Pocock to create an API to simplify development of Java-based Bioinformatics tools. BioJava is an active open source project that has been developed over more than 12 years and by more than 60 developers. BioJava is one of a number of Bio* projects designed to reduce code duplication.[5] Examples of such projects that fall under Bio* apart from BioJava are BioPython,[6] BioPerl,[7] BioRuby,[8] EMBOSS[9] etc.
In October 2012, the first paper on BioJava was published.[10] This paper detailed BioJava's modules, functionalities, and purpose.
As of November 2018 Google Scholar counts more than 130 citations.[11]
The most recent paper on BioJava was written in February 2017.[12] This paper detailed a new tool named BioJava-ModFinder. This tool can be used for identification and subsequent mapping of protein modifications to 3D in the Protein Data Bank (PBD). The package was also integrated with the RCSB PDB web application and added protein modification annotations to the sequence diagram and structure display. More than 30,000 structures with protein modifications were identified by using BioJava-ModFinder and can be found on the RCSB PDB website.
In the year 2008, BioJava's first Application note was published.[2] It was migrated from its original CVS repository to GitHub in April 2013.[13] The project has been moved to a separate repository, BioJava-legacy, and is still maintained for minor changes and bug fixes.[14]
Version 3 was released in December 2010. It was a major update to the prior versions. The aim of this release was to rewrite BioJava so that it could be modularized into small, reusable components. This allowed developers to contribute more easily and reduced dependencies. The new approach seen in BioJava 3 was modeled after the Apache Commons.
Version 4 was released in January 2015. This version brought many new features and improvements to the packages biojava-core, biojava-structure, biojava-structure-gui, biojava-phylo, as well as others. BioJava 4.2.0 was the first release to be available using Maven from the Maven Central.
Version 5 was released in March 2018. This represents a major milestone for the project. BioJava 5.0.0 is the first released based on Java 8 which introduces the use of lambda functions and streaming API calls. There were also major changes to biojava-structure module. Also, the previous data models for macro-molecular structures have been adapted to more closely represent the mmCIF data model. This was the first release in over two years. Some of the other improvements include optimizations in the biojava-structure module to improve symmetry detection and added support for MMTF formats. Other general improvements include Javadoc updates, dependency versions, and all tests are now Junit4. The release contains 1,170 commits from 19 contributors.
Modules
During 2014-2015, large parts of the original code base were rewritten. BioJava 3 is a clear departure from the version 1 series. It now consists of several independent modules built using an automation tool called Apache Maven.[15] These modules provide state-of-the-art tools for protein structure comparison, pairwise and multiple sequence alignments, working with DNA and protein sequences, analysis of amino acid properties, detecting protein modifications, predicting disordered regions in proteins, and parsers for common file formats using a biologically meaningful data model. The original code has been moved into a separate BioJava legacy project, which is still available for backward compatibility.[16]
BioJava 5 introduced new features to two modules, biojava-alignment and biojava-structure.
The following sections will describe several of the new modules and highlight some of the new features that are included in the latest version of BioJava.
![](../I/BioJava_5_Module_Layout.png.webp)
Core Module
This module provides Java classes to model amino acid or nucleotide sequences. The classes were designed so that the names are familiar and make sense to biologists and also provide a concrete representation of the steps in going from a gene sequence to a protein sequence for computer scientists and programmers.
A major change between the legacy BioJava project and BioJava3 lies in the way framework has been designed to exploit then-new innovations in Java. A sequence is defined as a generic interface allowing the rest of the modules to create any utility that operates on all sequences. Specific classes for common sequences such as DNA and proteins have been defined in order to improve usability for biologists. The translation engine really leverages this work by allowing conversions between DNA, RNA and amino acid sequences. This engine can handle details such as choosing the codon table, converting start codons to methionine, trimming stop codons, specifying the reading frame and handing ambiguous sequences.
Special attention has been paid to designing the storage of sequences to minimize space needs. Special design patterns such as the Proxy pattern allowed the developers to create the framework such that sequences can be stored in memory, fetched on demand from a web service such as UniProt, or read from a FASTA file as needed. The latter two approaches save memory by not loading sequence data until it is referenced in the application. This concept can be extended to handle very large genomic datasets, such as NCBI GenBank or a proprietary database.
Protein structure modules
The protein structure modules provide tools to represent and manipulate 3D biomolecular structures. They focus on protein structure comparison.
The following algorithms have been implemented and included in BioJava.
- FATCAT algorithm for flexible and rigid body alignment.[17]
- The standard Combinatorial Extension (CE) algorithm.[18]
- A new version of CE that can detect circular permutations in proteins.[19]
These algorithms are used to provide the RCSB Protein Data Bank (PDB)[20] Protein Comparison Tool as well as systematic comparisons of all proteins in the PDB on a weekly basis.[21]
Parsers for PDB[22] and mmCIF[23] file formats allow the loading of structure data into a reusable data model. This feature is used by the SIFTS project to map between UniProt sequences and PDB structures.[24] Information from the RCSB PDB can be dynamically fetched without the need to manually download data. For visualization, an interface to the 3D viewer Jmol is provided.[4]
Genome and Sequencing modules
This module is focused on the creation of gene sequence objects from the core module. This is realized by supporting the parsing of the following popular standard file formats generated by open source gene prediction applications:
- GTF files generated by GeneMark[25]
- GFF2 files generated by GeneID[26]
- GFF3 files generated by Glimmer[27]
Then the gene sequence objects are written out as a GFF3 format and is imported into GMOD.[28] These file formats are well defined but what gets written in the file is very flexible.
For providing input-output support for several common variants of the FASTQ file format from the next generation sequencers,[29] a separate sequencing module is provided. For samples on how to use this module please go to this link.
Alignment module
This module contains several classes and methods that allow users to perform pairwise and multiple sequence alignment. Sequences can be aligned in both a single and multi-threaded fashion. BioJava implements the Needleman-Wunsch[30] algorithm for optimal global alignments and the Smith and Waterman's[31] algorithm for local alignments. The outputs of both local and global alignments are available in standard formats. In addition to these two algorithms, there is an implementation of Guan–Uberbacher algorithm[32] which performs global sequence alignment very efficiently since it only uses linear memory.
For Multiple Sequence Alignment, any of the methods discussed above can be used to progressively perform a multiple sequence alignment.
ModFinder module
![](../I/An_example_application_using_the_ModFinder_module_and_the_protein_structure_module.png.webp)
The ModFinder module provides new methods to identify and classify protein modifications in protein 3D structures. Over 400 different types of protein modifications such as phosphorylation, glycosylation, disulfide bonds metal chelation etc. were collected and curated based on annotations in PSI-MOD,[34] RESID[35] and RCSB PDB.[36] The module also provides an API for detecting pre-, co-, and post-translational protein modifications within protein structures. This module can also identify phosphorylation and print all pre-loaded modifications from a structure.
Amino acid properties module
This module attempts to provide accurate physio-chemical properties of proteins. The properties that can calculated using this module are as follows:
- Molecular mass
- Extinction coefficient
- Instability index
- Aliphatic index
- Grand average of hydropathy
- Isoelectric point
- Amino acid composition
The precise molecular weights for common isotopically labelled amino acids are included in this module. There also exists flexibility to define new amino acid molecules with their molecular weights using simple XML configuration files. This can be useful where the precise mass is of high importance such as mass spectrometry experiments.
Protein disorder module
The goal of this module is to provide users ways to find disorders in protein molecules. BioJava includes a Java implementation of the RONN predictor. The BioJava 3.0.5 makes use of Java's support for multithreading to improve performance by up to 3.2 times,[37] on a modern quad-core machine, as compared to the legacy C implementation.
There are two ways to use this module:
- Using library function calls
- Using command line
Some features of this module include:
- Calculating the probability of disorder for every residue in a sequence
- Calculating the probability of disorder for every residue in the sequence for all proteins from a FASTA input file
- Get the disordered regions of the protein for a single protein sequence or for all the proteins from a FASTA input file
Web service access module
As per the current trends in bioinformatics, web based tools are gaining popularity. The web service module allows bioinformatics services to be accessed using REST protocols. Currently, two services are implemented: NCBI Blast through the Blast URLAPI (previously known as QBlast) and the HMMER web service.[38]
Comparisons with other alternatives
The need for customized software in the field of bioinformatics has been addressed by several groups and individuals. Similar to BioJava, open-source software projects such as BioPerl, BioPython, and BioRuby all provide tool-kits with multiple functionality that make it easier to create customized pipelines or analysis.
As the names suggest, the projects mentioned above use different programming languages. All of these APIs offer similar tools so on what criteria should one base their choice? For programmers who are experienced in only one of these languages, the choice is straightforward. However, for a well-rounded bioinformaticist who knows all of these languages and wants to choose the best language for a job, the choice can be made based on the following guidelines given by a software review done on the Bio* tool-kits.[5]
In general, for small programs (<500 lines) that will be used by only an individual or small group, it is hard to beat Perl and BioPerl. These constraints probably cover the needs of 90 per cent of personal bioinformatics programming.
For beginners, and for writing larger programs in the Bio domain, especially those to be shared and supported by others, Python’s clarity and brevity make it very attractive.
For those who might be leaning towards a career in bioinformatics and who want to learn only one language, Java has the widest general programming support, very good support in the Bio domain with BioJava, and is now the de facto language of business (the new COBOL, for better or worse).
Apart from these Bio* projects there is another project called STRAP which uses Java and aims for similar goals. The STRAP-toolbox, similar to BioJava is also a Java-toolkit for the design of Bioinformatics programs and scripts. The similarities and differences between BioJava and STRAP are as follows:
Similarities
- Both provide comprehensive collections of methods for protein sequences.
- Both are used by Java programmers to code bioinformatics algorithms.
- Both separate implementations and definitions by using java interfaces.
- Both are open source projects.
- Both can read and write many sequence file formats.
Differences
- BioJava is applicable to nucleotide and peptide sequences and can be applied for entire genomes. STRAP cannot cope with single sequences as long as an entire chromosome. Instead STRAP manipulates peptide sequences and 3D- structures of the size of single proteins. Nevertheless, it can hold a high number of sequences and structures in memory. STRAP is designed for protein sequences but can read coding nucleotide files, which are then translated to peptide sequences.
- STRAP is very fast since the graphical user interface must be highly responsive. BioJava is used where speed is less critical.
- BioJava is well designed in terms of type safety, ontology and object design. BioJava uses objects for sequences, annotations and sequence positions. Even single amino acids or nucleotides are object references. To enhance speed, STRAP avoids frequent object instantiations and invocation of non-final object-methods.
- In BioJava peptide sequences and nucleotide sequences are lists of symbols. The symbols can be retrieved one after the other with an iterator or sub-sequences can be obtained. The advantages are that the entire sequence does not necessarily reside in memory and that programs are less susceptible to programming errors. Symbol objects are immutable elements of an alphabet. In STRAP however simple byte arrays are used for sequences and float arrays for coordinates. Besides speed the low memory consumption is an important advantage of basic data types. Classes in Strap expose internal data. Therefore, programmers might commit programming errors like manipulating byte arrays directly instead of using the setter methods. Another disadvantage is that no checks are performed in STRAP whether the characters in sequences are valid with respect to an underlying alphabet.
- In BioJava sequence positions are realized by the class Location. Discontiguous Location objects are composed of several contiguous RangeLocation objects or PointLocation objects. For the class StrapProtein however, single residue positions are indicated by integer numbers between 0 and countResidues()-1. Multiple positions are given by boolean arrays. True at a given index means selected whereas false means not selected.
- BioJava throws exceptions when methods are invoked with invalid parameters. STRAP avoids the time-consuming creation of Throwable objects. Instead, errors in methods are indicated by the return values NaN, -1 or null. From the point of program design however Throwable objects are nicer.
- In BioJava a Sequence object is either a peptide sequence or a nucleotide sequence. A StrapProtein can hold both at the same time if a coding nucleotide sequence was read and translated into protein. Both, the nucleotide sequence and the peptide sequence are contained in the same StrapProtein object. The coding or non-coding regions can be changed and the peptide sequence alters accordingly.
Projects using BioJava
The following projects make use of BioJava.
- Metabolic Pathway Builder: Software suite dedicated to the exploration of connections among genes, proteins, reactions and metabolic pathways
- DengueInfo Archived 2006-12-08 at the Wayback Machine: a Dengue genome information portal that uses BioJava in the middleware and talks to a biosql database.
- Dazzle: A BioJava based DAS server.
- BioSense: A plug-in for the InforSense Suite, an analytics software platform by IDBS that unitizes BioJava.
- Bioclipse: A free, open source, workbench for chemo- and bioinformatics with powerful editing and visualizing abilities for molecules, sequences, proteins, spectra, etc.
- PROMPT: A free, open source framework and application for the comparison and mapping of protein sets. Uses BioJava for handling most input data formats.
- Cytoscape: An open source bioinformatics software platform to visualize molecular interaction networks.
- BioWeka: An open source biological data mining application.
- Geneious: A molecular biology toolkit.
- MassSieve: An open source application to analyze mass spec proteomics data.
- STRAP: A tool for multiple sequence alignment and sequence-based structure alignment.
- Jstacs: A Java framework for statistical analysis and classification of biological sequences
- jLSTM: "Long Short-Term Memory" for protein classification
- LaJolla: An open source structural alignment tool for RNA and proteins using an index structure for fast alignment of thousands of structures; includes an easy-to-use command line interface.
- GenBeans: A rich client platform for bioinformatics primarily focused on molecular biology and sequence analysis.
- JEnsembl: A version-aware Java API to Ensembl data systems.[39]
- MUSI: An integrated system to identify multiple specificity from very large peptide or nucleic acid data sets.[40]
- Bioshell: A utility library for structural bioinformatics[41]
See also
References
- Prlić A, Yates A, Bliven SE, et al. (October 2012). "BioJava: an open-source framework for bioinformatics in 2012". Bioinformatics. 28 (20): 2693–5. doi:10.1093/bioinformatics/bts494. PMC 3467744. PMID 22877863.
- Holland RC, Down TA, Pocock M, Prlić A, Huen D, James K, et al. (2008). "BioJava: an open-source framework for bioinformatics". Bioinformatics. 24 (18): 2096–7. doi:10.1093/bioinformatics/btn397. PMC 2530884. PMID 18689808.
- VS Matha and P Kangueane, 2009, Bioinformatics: a concept-based introduction, 2009. p26
- Hanson, R.M. (2010) Jmol a paradigm shift in crystallographic visualization.
- Mangalam H (2002). "The Bio* toolkits--a brief overview". Briefings in Bioinformatics. 3 (3): 296–302. doi:10.1093/bib/3.3.296. PMID 12230038.
- Cock PJ, Antao T, Chang JT, et al. (June 2009). "Biopython: freely available Python tools for computational molecular biology and bioinformatics". Bioinformatics. 25 (11): 1422–3. doi:10.1093/bioinformatics/btp163. PMC 2682512. PMID 19304878.
- Stajich JE, Block D, Boulez K, et al. (October 2002). "The Bioperl toolkit: Perl modules for the life sciences". Genome Res. 12 (10): 1611–8. doi:10.1101/gr.361602. PMC 187536. PMID 12368254.
- Goto N, Prins P, Nakao M, Bonnal R, Aerts J, Katayama T (October 2010). "BioRuby: bioinformatics software for the Ruby programming language". Bioinformatics. 26 (20): 2617–9. doi:10.1093/bioinformatics/btq475. PMC 2951089. PMID 20739307.
- Rice P, Longden I, Bleasby A (June 2000). "EMBOSS: the European Molecular Biology Open Software Suite". Trends Genet. 16 (6): 276–7. doi:10.1016/S0168-9525(00)02024-2. PMID 10827456.
- Prlić A, Yates A, Bliven SE, et al. (October 2012). "BioJava: an open-source framework for bioinformatics in 2012". Bioinformatics. 28 (20): 2693–5. doi:10.1093/bioinformatics/bts494. PMC 3467744. PMID 22877863.
- "Google Scholar". scholar.google.com. Retrieved 2018-11-22.
- Gao, Jianjiong; Prlić, Andreas; Bi, Chunxiao; Bluhm, Wolfgang F.; Dimitropoulos, Dimitris; Xu, Dong; Bourne, Philip E.; Rose, Peter W. (2017-02-17). "BioJava-ModFinder: identification of protein modifications in 3D structures from the Protein Data Bank". Bioinformatics. 33 (13): 2047–2049. doi:10.1093/bioinformatics/btx101. ISSN 1367-4803. PMC 5870676. PMID 28334105.
- "History". Retrieved 30 Jan 2015.
- BioJava-legacy Archived 2013-01-09 at the Wayback Machine
- Maven, Apache. "Maven". Apache.
- BioJava legacy project Archived 2013-01-09 at the Wayback Machine
- Ye Y, Godzik A (October 2003). "Flexible structure alignment by chaining aligned fragment pairs allowing twists". Bioinformatics. 19 (Suppl 2): ii246–55. doi:10.1093/bioinformatics/btg1086. PMID 14534198.
- Shindyalov IN, Bourne PE (September 1998). "Protein structure alignment by incremental combinatorial extension (CE) of the optimal path". Protein Eng. 11 (9): 739–47. doi:10.1093/protein/11.9.739. PMID 9796821.
- Bliven S, Prlić A (2012). "Circular permutation in proteins". PLOS Comput. Biol. 8 (3): e1002445. Bibcode:2012PLSCB...8E2445B. doi:10.1371/journal.pcbi.1002445. PMC 3320104. PMID 22496628.
- Rose PW, Beran B, Bi C, et al. (January 2011). "The RCSB Protein Data Bank: redesigned web site and web services". Nucleic Acids Res. 39 (Database issue): D392–401. doi:10.1093/nar/gkq1021. PMC 3013649. PMID 21036868.
- Prlić A, Bliven S, Rose PW, et al. (December 2010). "Pre-calculated protein structure alignments at the RCSB PDB website". Bioinformatics. 26 (23): 2983–5. doi:10.1093/bioinformatics/btq572. PMC 3003546. PMID 20937596.
- Bernstein FC, Koetzle TF, Williams GJ, et al. (May 1977). "The Protein Data Bank: a computer-based archival file for macromolecular structures". J. Mol. Biol. 112 (3): 535–42. doi:10.1016/s0022-2836(77)80200-3. PMID 875032.
- Fitzgerald, P.M.D. et al. (2006) Macromolecular dictionary (mmCIF). In Hall, S.R.
- Velankar S, McNeil P, Mittard-Runte V, et al. (January 2005). "E-MSD: an integrated data resource for bioinformatics". Nucleic Acids Res. 33 (Database issue): D262–5. doi:10.1093/nar/gki058. PMC 540012. PMID 15608192.
- Besemer J, Borodovsky M (July 2005). "GeneMark: web software for gene finding in prokaryotes, eukaryotes and viruses". Nucleic Acids Res. 33 (Web Server issue): W451–4. doi:10.1093/nar/gki487. PMC 1160247. PMID 15980510.
- Blanco E, Abril JF (2009). "Computational gene annotation in new genome assemblies using GeneID". Bioinformatics for DNA Sequence Analysis. Methods in Molecular Biology. Vol. 537. pp. 243–61. doi:10.1007/978-1-59745-251-9_12. ISBN 978-1-58829-910-9. PMID 19378148.
- Kelley DR, Liu B, Delcher AL, Pop M, Salzberg SL (January 2012). "Gene prediction with Glimmer for metagenomic sequences augmented by classification and clustering". Nucleic Acids Res. 40 (1): e9. doi:10.1093/nar/gkr1067. PMC 3245904. PMID 22102569.
- Stein LD, Mungall C, Shu S, et al. (October 2002). "The generic genome browser: a building block for a model organism system database". Genome Res. 12 (10): 1599–610. doi:10.1101/gr.403602. PMC 187535. PMID 12368253.
- Cock PJ, Fields CJ, Goto N, Heuer ML, Rice PM (April 2010). "The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants". Nucleic Acids Res. 38 (6): 1767–71. doi:10.1093/nar/gkp1137. PMC 2847217. PMID 20015970.
- Needleman SB, Wunsch CD (March 1970). "A general method applicable to the search for similarities in the amino acid sequence of two proteins". J. Mol. Biol. 48 (3): 443–53. doi:10.1016/0022-2836(70)90057-4. PMID 5420325.
- Smith TF, Waterman MS (March 1981). "Identification of common molecular subsequences". J. Mol. Biol. 147 (1): 195–7. CiteSeerX 10.1.1.63.2897. doi:10.1016/0022-2836(81)90087-5. PMID 7265238.
- Guan X, Uberbacher EC (February 1996). "Alignments of DNA and protein sequences containing frameshift errors". Comput. Appl. Biosci. 12 (1): 31–40. doi:10.1093/bioinformatics/12.1.31. PMID 8670617.
- Chen K, Jung YS, Bonagura CA, et al. (February 2002). "Azotobacter vinelandii ferredoxin I: a sequence and structure comparison approach to alteration of [4Fe-4S]2+/+ reduction potential". J. Biol. Chem. 277 (7): 5603–10. doi:10.1074/jbc.M108916200. PMID 11704670.
- Montecchi-Palazzi L, Beavis R, Binz PA, et al. (August 2008). "The PSI-MOD community standard for representation of protein modification data". Nat. Biotechnol. 26 (8): 864–6. doi:10.1038/nbt0808-864. PMID 18688235. S2CID 205270043.
- Garavelli JS (June 2004). "The RESID Database of Protein Modifications as a resource and annotation tool". Proteomics. 4 (6): 1527–33. doi:10.1002/pmic.200300777. PMID 15174122. S2CID 25712150.
- Berman HM, Westbrook J, Feng Z, et al. (January 2000). "The Protein Data Bank". Nucleic Acids Res. 28 (1): 235–42. doi:10.1093/nar/28.1.235. PMC 102472. PMID 10592235.
- Yang ZR, Thomson R, McNeil P, Esnouf RM (August 2005). "RONN: the bio-basis function neural network technique applied to the detection of natively disordered regions in proteins". Bioinformatics. 21 (16): 3369–76. doi:10.1093/bioinformatics/bti534. PMID 15947016.
- Finn RD, Clements J, Eddy SR (July 2011). "HMMER web server: interactive sequence similarity searching". Nucleic Acids Res. 39 (Web Server issue): W29–37. doi:10.1093/nar/gkr367. PMC 3125773. PMID 21593126.
- Paterson T, Law A (November 2012). "JEnsembl: a version-aware Java API to Ensembl data systems". Bioinformatics. 28 (21): 2724–31. doi:10.1093/bioinformatics/bts525. PMC 3476335. PMID 22945789.
- Kim T, Tyndel MS, Huang H, et al. (March 2012). "MUSI: an integrated system for identifying multiple specificity from very large peptide or nucleic acid data sets". Nucleic Acids Res. 40 (6): e47. doi:10.1093/nar/gkr1294. PMC 3315295. PMID 22210894.
- Gront D, Kolinski A (February 2008). "Utility library for structural bioinformatics". Bioinformatics. 24 (4): 584–5. doi:10.1093/bioinformatics/btm627. PMID 18227118.