esl-alimask
Quick Start
- Command:
esl-alimask [options] <msafile> ... - Local executable:
/home/vimalinx/miniforge3/envs/bio/bin/esl-alimask - Detailed reference: See
references/help.mdfor full option list and usage details
When To Use This Tool
- Use
esl-alimaskwhen you need to remove or retain alignment columns based on an explicit mask, coordinate range, gap fraction, posterior probability, or RF annotation. - It is appropriate for MSA cleanup before HMM building, consensus trimming, or downstream comparative analyses.
- Use
-tfor coordinate-based trimming,-gfor gap-based masking, and-pfor posterior-probability-based masking. - Reach for
--rf-is-maskwhen the alignment already carries a trusted#=GC RFannotation defining the columns to keep.
Common Patterns
# Keep only columns 23..100 of an alignment
esl-alimask -t alignment.sto 23..100 > trimmed.sto
# Mask columns with too many gaps
esl-alimask -g --gapthresh 0.3 alignment.sto > gapmasked.sto
# Mask using posterior probabilities
esl-alimask -p --pavg 0.9 alignment.sto > ppmasked.sto
# Apply an explicit 0/1 mask file
esl-alimask alignment.sto mask.txt > masked.sto
Recommended Workflow
- Prepare your MSA file in a supported format (e.g., Stockholm, aligned FASTA)
- Choose a masking strategy: external mask file, gap-based, posterior probability-based, or RF annotation
- Run
esl-alimaskwith the appropriate mode flag (-t,-g,-p, or--rf-is-mask) or provide a mask file - Verify the output alignment has the expected columns masked or truncated
Guardrails
- Use
-hfor help;--helpand--versionare not supported options - Ensure the MSA file format matches the expected input (check with Easel tools if needed)
- When using
-t, provide coordinates in the correct format and range for the alignment -gand-pcan be combined, but the other major usage modes are mutually exclusive