CD-HIT-454(1)

๐Ÿ“ฆ
version 4.8.1 ยท 2 versions
, in debian bookworm, debian trixie, ubuntu jammy, ubuntu noble

๐Ÿ’– NAME ๐Ÿ’–

cd-hit-454 โ€” quickly group sequences, optimised for 454 data ๐Ÿฆ‹


โœจ SYNOPSIS โœจ

cd-hit-454 [Options]

๐ŸŽ€ DESCRIPTION ๐ŸŽ€

okay bestie, first thing you see when you run this is the little banner flexing its whole build identity, and honestly? iconic:

====== CD-HIT version 4.8.1 (built on Aug 20 2021) ======

๐Ÿ’… Options ๐Ÿ’…

  • -i โ€” input filename in fasta format, and yes it's required. Can be in .gz format too, we love a compressed queen ๐Ÿซถ
  • -o โ€” output filename, also required, no you cannot skip it
  • -c โ€” sequence identity threshold, default 0.98. This is a "global sequence identity" and it's calculated as: number of identical amino acids or bases in the alignment divided by the full length of the shorter sequence + gaps โœจ
  • -b โ€” band_width of the alignment, default 10
  • -M โ€” memory limit in MB for the program, default 800; set it to 0 for unlimited ๐Ÿ‘‘
  • -T โ€” number of threads, default 1; if you set it to 0, all CPUs will be used ๐Ÿ”ฅ
  • -n โ€” word_length, default 10, and bestie if you're wondering how to pick it, go read the user's guide, it's giving required reading ๐Ÿ“–
  • -aL โ€” alignment coverage for the longer sequence, default 0.0. Set it to 0.9 and the alignment must cover 90% of the sequence ๐Ÿ’…
  • -AL โ€” alignment coverage control for the longer sequence, default 99999999. So if you set it to 60 and the length of the sequence is 400, then the alignment must be >= 340 (400-60) residues, math ate that ๐Ÿงฎ
  • -aS โ€” alignment coverage for the shorter sequence, default 0.0. Set it to 0.9 and the alignment must cover 90% of the sequence ๐Ÿ’…
  • -AS โ€” alignment coverage control for the shorter sequence, default 99999999. If you set it to 60 and the length of the sequence is 400, then the alignment must be >= 340 (400-60) residues, same energy as its longer-sequence sibling ๐Ÿงฎ
  • -B โ€” 1 or 0, default 0. By default sequences are stored in RAM; set it to 1 and sequences are stored on the hard drive.
๐Ÿ˜ญ !! No longer supported !! do not come for me, I'm just the messenger
  • -g โ€” 1 or 0, default 0. With cd-hit's default algorithm, a sequence gets clustered to the first cluster that meets the threshold (fast cluster). Set it to 1 and the program clusters it into the most similar cluster that meets the threshold (accurate but slow mode). But lowkey, either 1 or 0 won't change the representatives of the final clusters, so don't stress ๐ŸŒธ
  • -D โ€” max size per indel, default 1
  • -match โ€” matching score, default 2
  • -mismatch โ€” mismatching score, default -1
  • -gap โ€” gap opening score, default -3
  • -gap-ext โ€” gap extension score, default -1
  • -bak โ€” write backup cluster file (1 or 0, default 0)
  • -h โ€” print this help, the classic ๐ŸŽ€

๐Ÿซถ Questions, bugs, contact ๐Ÿซถ

Got questions? Got bugs? Contact Weizhong Li at liwz@sdsc.edu ๐Ÿ“ง

And if you find cd-hit useful, please kindly cite these, it's the polite thing to do fr fr:

  1. "CD-HIT: a fast program for clustering and comparing large sets of protein or nucleotide sequences", Weizhong Li & Adam Godzik. Bioinformatics, (2006) 22:1658-1659
  2. "CD-HIT: accelerated for clustering the next generation sequencing data", Limin Fu, Beifang Niu, Zhengwei Zhu, Sitao Wu & Weizhong Li. Bioinformatics, (2012) 28:3150-3152
  3. "Beifang Niu, Limin Fu, Shulei Sun and Weizhong Li. Artificial and natural duplicates in pyrosequencing reads of metagenomic data. BMC Bioinformatics (2010) 11:187"

and that's the whole page bestie, cd-hit-454 version 4.8.1, August 2021, page CD-HIT-454(1) โ€” go cluster something gorgeous ๐Ÿ’œ๐ŸŒˆโญ

๐Ÿ“ฆ where this page comes from