GETCDS(1)
πΈ NAME πΈ
okay so bestie, meet getcds.pl (a Perl situation, obviously) β it does CDS extraction from various annotation formats π
like it just. goes and gets your coding sequences out of whatever annotation chaos you handed it. iconic behavior, honestly.
β¨ SYNOPSIS β¨
the vibe check, literally:
getcds.pl [options] input file [output file] [input file [outputfile] ...]
so yeah, it takes options, then an input file, then optionally an output file, and then it can just. keep going. pairs of files, over and over, as many as you want π«Ά
π OPTIONS π
the options, as promised, include:
redirect=> turns on a redirect at runtime. and yes, it's explained below, keep reading babe, don't scroll past me πforce=> lets you write empty.cdsfiles. that's it, that's the whole thing. sometimes empty is a valid lifestyle β¨
π tiny but important: if you don't specify an [output file], it just uses [input file].cds. so it AutoPicks for you. no cap.
π¦ the redirect lore π¦
there's a lot of dark magic you can do with .redirect files π. redirects exist to modify the ordinary behavior of the program, and you can apply them:
- globally β a
.redirectfile sitting in the same directory as [input file] π - per-input β a
[input file].redirectfile, just for that one file π
and the order matters: they get read in that order, so per-input effects get applied on top of the global effects. global is the base, per-input is the glitter on top π
also, a cute little convenience: if you left the filename exactly as-is from an Ensembl download, the right Ensembl DB source gets guessed automatically for you. no extra work. we love a self-sufficient file π
π§ possible redirects π§
here's the whole menu, bestie, every single one:
CDSβ changes what tag type gets searched for when running on bioperl files. (the default isCDS, in case you were wondering.)offsetβ adds some value to all coordinates. shift the whole thing, why not β¨+bioperlβ on top of any other redirect-based data, this also extracts annotation via bioperl (using a direct disables bioperl's parser by default). yes, the+is doing something. it's giving additive π₯ensembl_buildβ grab data from the ensembl (in the specified schema).ucsc_buildβ grab data from ucsc (in the specified schema).murasaki_synthβ synthesize annotation data for each anchor from an alignment. and if this is a directory, then each input file (eginput.lav) is checked for a corresponding.anchorsfile (eg.input.anchors). very "we'll find it ourselves" of it πprimaryβ what type of tags to create from data gathered from outside data (eg:ensembl,ucsc, ormurasaki). the default is the same as theCDSredirect. so like, it inherits the vibe πchromosomeβ forces a specific chromosome rather than deriving it from the filename. assertive. we stan πgtfβ extract annotation from a.gtffile. that's the whole redirect, lowkey perfect π¦
π warning-ish energy, but in a nice way: the redirect files are genuinely dark magic territory. you can stack them, and per-input beats global, so if your output looks cursed, check your .redirect files first bestie. not me speaking from experience π
and that's the whole page, fr fr π«Ά
sign-off moment, because you deserve it: this lovely thing is documented as part of the User Contributed Perl Documentation β¨, running under perl v5.40.1 π, and this version of the page is dated 2026-08-12. go extract some CDS, you absolute main character ππΈ
π¦ where this page comes from
- ubuntu stonking: murasaki-common 1.68.6-18 Β· GPL-3.0-or-later
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms