seqstat(1)

πŸ’… seqstat(1) β€” the Biosquid Manual, but make it cute ✨

okay bestie, buckle up, because this is the little tool that looks at your sequence file and just tells you things. it's from the Biosquid Manual and it's lowkey iconic. πŸ¦‹


🌸 NAME 🌸

seqstat β€” show statistics and format for a sequence file πŸ‘‘


✨ SYNOPSIS ✨

seqstat [options] seqfile

πŸ’– DESCRIPTION πŸ’–

so basically, seqstat reads a sequence file, seqfile, and shows you a bunch of simple stats about it. no cap, that's the whole job. πŸŽ€

the stats it prints out are:

  • the name of the format πŸ“
  • the residue type of the first sequence β€” so like, protein, RNA, or DNA 🧬
  • the number of sequences πŸ”’
  • the total number of residues 🌈
  • and the average and range of the sequence lengths πŸ“

it's giving "vibe check but for your FASTA file," fr fr.


πŸŽ€ OPTIONS πŸŽ€

  • -a β€” show additional verbose info: a table with one line per sequence showing name, length, and description line. these lines are prefixed with a * character so you can easily grep them out and sort them. it's literally built for piping, she understood the assignment. πŸ’…
  • -h β€” print brief help; includes the version number and a summary of all options, including expert options. bestie behaviour honestly 🫢
  • -B β€” (Babelfish 🐠). autodetect and read a sequence file format other than the default (FASTA). almost any common sequence file format is recognized (including Genbank, EMBL, SWISS-PROT, PIR, and GCG unaligned sequence formats, and Stockholm, GCG MSF, and Clustal alignment formats). see the printed documentation for a complete list of supported formats. she speaks so many languages, multilingual queen ✨

πŸ‘‘ EXPERT OPTIONS πŸ‘‘

  • --informat <s> β€” specify that the sequence file is in format <s>, rather than the default FASTA format. common examples include Genbank, EMBL, GCG, PIR, Stockholm, Clustal, MSF, or PHYLIP; see the printed documentation for a complete list of accepted format names. this option overrides the default expected format (FASTA) and the -B Babelfish autodetection option. so like, it's the final boss of format flags, it said "i'll handle this" πŸ’€
  • --quiet β€” suppress the verbose header (program name, release number and date, the parameters and options in effect). for when you want the tea but not the small talk 🀫

🌈 SEE ALSO 🌈

afetch(1), alistat(1), compalign(1), compstruct(1), revcomp(1), seqsplit(1), sfetch(1), shuffle(1), sindex(1), sreformat(1), stranslate(1), weight(1).

the whole friend group!! go say hi to all of them πŸ«ΆπŸ’œ


πŸ¦‹ AUTHOR πŸ¦‹

Biosquid and its documentation are Copyright (C) 1992-2003 HHMI/Washington University School of Medicine. Freely distributed under the GNU General Public License (GPL). See COPYING in the source code distribution for more details, or contact me.

Sean Eddy HHMI/Department of Genetics Washington University School of Medicine 4444 Forest Park Blvd., Box 8510 St Louis, MO 63108 USA Phone: 1-314-362-7666 FAX : 1-314-362-2157 Email: eddy@genetics.wustl.edu


and that's the whole page, bestie!! πŸ’– this one shipped as Biosquid 1.9g, dated January 2003 β€” vintage but still eating, no crumbs left. now go run seqstat on something and feel powerful. periodt. βœ¨πŸŽ€πŸ¦‹