Blast score-based single-linkage clustering
blastclust [\|-\|] [\|-C\|] [\|-L X\|] [\|-S X\|] [\|-W N\|] [\|-a N\|] [\|-b F\|] [\|-c filename\|] [\|-d filename\|] [\|-e F\|] [\|-i filename\|] [\|-l filename\|] [\|-o filename\|] [\|-p F\|] [\|-r filename\|] [\|-s filename\|] [\|-v [\|filename\|]\|]
blastclust automatically and systematically clusters protein or DNA sequences based on pairwise matches found using the BLAST algorithm in case of proteins or Mega BLAST algorithm for DNA. In the latter case a single Mega BLAST search is performed for all the sequences combined against a database created from the same sequences. blastclust finds pairs of sequences that have statistically significant matches and clusters them using single-linkage clustering.
A summary of options is included below.
Print usage message
Complete unfinished clustering
Length coverage threshold (default = 0.9)
Score coverage threshold (bit score / length if < 3.0, percentage of identities otherwise; default = 1.75)
Use words of size N (length of best perfect match; zero invokes default behavior: 3 for proteins, 32 for nucleotides)
Number of CPU's to use (default = 1)
Do not require coverage on both neighbours
Read advanced options from configuration file filename
Input as a database
Disable id parsing in database formatting
FASTA input file (program will format the database and remove files in the end; default = stdin)
Restrict reclustering to id list in filename
Output file for list of clusters (default = stdout)
Input is nucleotides, not proteins.
Restore neighbors for reclustering from filename
Save all neighbours to filename
Print verbose progress messages (to filename)
The National Center for Biotechnology Information.
blast(1), formatdb(1), /usr/share/doc/blast2/blastclust.html, <http://www.ncbi.nlm.nih.gov/BLAST/>