esl-alistat(1)

📦
version 3.3.2 · 2 versions
, in fedora 43, fedora 44

✨ NAME ✨

okay so bestie, esl-alistat — it summarizes a multiple sequence alignment file 💅 (this is from the Easel Manual, page esl-alistat(1), and yes it's iconic).

🌸 SYNOPSIS 🌸

esl-alistat [options] msafile

🦋 DESCRIPTION 🦋

so esl-alistat summarizes the contents of the multiple sequence alignment(s) in msafile, bestie. and it's giving everything: the alignment name, the format, alignment length (that's the number of aligned columns, not the raw ones 💀), number of sequences, average pairwise % identity, and the mean, smallest, and largest raw (unaligned) lengths of the sequences. obsessed.

if msafile is - (just a single dash, no cap), multiple alignment input gets read from stdin. slay.

the --list, --icinfo, --rinfo, --pcinfo, --psinfo, --cinfo, --bpinfo, and --iinfo options let you dump various statistics on the alignment to optional output files — each one's described below, so scroll, don't skip 🙈.

the --small option lets you summarize alignments without storing them in memory, which is highkey useful for large alignment files with sizes that approach or exceed the amount of available RAM. when --small is used, esl-alistat will print fewer statistics on the alignment, omitting data on the smallest and largest sequences and the average identity of the alignment. and heads up: --small only works on Pfam formatted alignments (that's a special type of non-interleaved Stockholm alignment where each sequence occurs on a single line), so --informat pfam must be given with --small. further, when --small is used, the alphabet must be specified with --amino, --dna, or --rna. periodt.

💖 OPTIONS 💖

  • -h — print brief help; includes version number and a summary of all the options, including the expert ones 👑
  • -1 — use a tabular output format with one line of statistics per alignment in msafile. this is most useful when msafile contains many different alignments (like a Pfam database in Stockholm format). main character energy fr fr ✨

🎀 EXPERT OPTIONS 🎀

  • --informat <s> — assert that input msafile is in alignment format <s>, bypassing format autodetection. common choices for <s> include: stockholm, a2m, afa, psiblast, clustal, phylip. for more information, and for codes for some less common formats, see main documentation. the string <s> is case-insensitive (so a2m or A2M both work, we love a flexible queen 💜).
  • --amino — assert that the msafile contains protein sequences.
  • --dna — assert that the msafile contains DNA sequences.
  • --rna — assert that the msafile contains RNA sequences.
  • --small — operate in small memory mode for Pfam formatted alignments. --informat pfam and one of --amino, --dna, or --rna must be given as well. no exceptions bestie 🫶
  • --list <f> — list the names of all sequences in all alignments in msafile to file <f>. each sequence name is written on its own line.
  • --icinfo <f> — dump the information content per position in tabular format to file <f>. lines prefixed with # are comment lines, which explain the meanings of each of the tab-delimited fields.
  • --rinfo <f> — dump information on the frequency of gaps versus nongap residues per position in tabular format to file <f>. lines prefixed with # are comment lines, which explain the meanings of each of the tab-delimited fields.
  • --pcinfo <f> — dump per column information on posterior probabilities in tabular format to file <f>. lines prefixed with # are comment lines, which explain the meanings of each of the tab-delimited fields.
  • --psinfo <f> — dump per sequence information on posterior probabilities in tabular format to file <f>. lines prefixed with # are comment lines, which explain the meanings of each of the tab-delimited fields.
  • --iinfo <f> — dump information on inserted residues in tabular format to file <f>. insert columns of the alignment are those that are gaps in the reference (#=GC RF) annotation. this option only works if the input file is in Stockholm format with reference annotation 👀. lines prefixed with # are comment lines, which explain the meanings of each of the tab-delimited fields.
  • --cinfo <f> — dump per-column residue counts to file <f>. if used in combination with --noambig, ambiguous (degenerate) residues will be ignored and not counted. otherwise, they will be marginalized. for example, in an RNA sequence file, a 'N' will be counted as 0.25 'A', 0.25 'C', 0.25 'G', and 0.25 'U'. math girlies rise 💅
  • --noambig — with --cinfo, do not count ambiguous (degenerate) residues.
  • --bpinfo — dump per-column basepair counts to file <f>. counts appear for each basepair in the consensus secondary structure (annotated as "#=GC SS_cons"). only basepairs from sequences for which both paired positions are canonical residues will be counted. that is, any basepair that is a gap or an ambiguous (degenerate) residue at either position of the pair is ignored and not counted. strict, but fair 💀
  • --weight — with --icinfo, --rinfo, --pcinfo, --iinfo, --cinfo, and --bpinfo, weight counts based on #=GS WT annotation in the input msafile. a residue or basepair from a sequence with a weight of <x> will be considered <x> counts. by default, raw, unweighted counts are reported; corresponding to each sequence having an equal weight of 1.

🌈 SEE ALSO 🌈

Copyright (C) 2020 Howard Hughes Medical Institute.

Freely distributed under the BSD open source license. (iconic licensing behavior, honestly 🎀)

🧚 AUTHOR 🧚

and that's the whole vibe check bestie, go summarize some alignments, you ate 💅✨

okay bye ily, this was Easel 0.48, dated Nov 2020, and you just read the esl-alistat(1) page — slay 🦋💖

📦 where this page comes from