alignment-thin - Remove sequences or columns from an alignment.
Contents
Column Filtering Options:
-Karg, –keep-columnsarg
Keep columns from this sequence
-marg, –min-lettersarg
Remove columns with fewer than arg letters.
-uarg, –remove-uniquearg
Remove insertions in a single sequence if longer than arg letters
-e, –erase-empty-columns
Remove columns with no characters (all gaps).
Description
Remove sequences or columns from an alignment.
Examples:
Remove columns without a minimum number of letters:
% alignment-thin --min-letters=5 file.fasta > file-thinned.fasta
Remove sequences by name:
% alignment-thin --remove=seq1,seq2 file.fasta > file2.fasta
% alignment-thin --keep=seq1,seq2 file.fasta > file2.fasta
Remove short sequences:
% alignment-thin --longer-than=250 file.fasta > file-long.fasta
Remove similar sequences with <= 5 differences from the closest other sequence:
% alignment-thin --cutoff=5 file.fasta > more-than-5-differences.fasta
Remove similar sequences until we have the right number of sequences:
% alignment-thin --down-to=30 file.fasta > file-30taxa.fasta
Remove dissimilar sequences that are missing conserved columns:
% alignment-thin --remove-crazy=10 file.fasta > file2.fasta
Protect some sequences from being removed:
% alignment-thin --down-to=30 file.fasta --protect=seq1,seq2 > file2.fasta
% alignment-thin --down-to=30 file.fasta --protect=@filename > file2.fasta
General Options:
-h, –help
Print usage information.
-V, –verbose
Output more log messages on stderr.
Name
alignment-thin - Remove sequences or columns from an alignment.
Output Options:
-S, –sort
Sort partially ordered columns to group similar gaps.
-L, –show-lengths
Just print out sequence lengths.
-N, –show-names
Just print out sequence lengths.
-Farg, –find-dupsarg
For each sequence, find the closest other sequence.
Reporting Bugs:
BAli-Phy online help: <http://www.bali-phy.org/docs.php>.
Please send bug reports to <bali-phy-users@googlegroups.com>.
Sequence Filtering Options:
-parg, –protectarg
Sequences that cannot be removed (comma-separated).
-karg, –keeparg
Remove sequences not in comma-separated list arg.
-rarg, –removearg
Remove sequences in comma-separated list arg.
-larg, –longer-thanarg
Remove sequences not longer than arg.
-sarg, –shorter-thanarg
Remove sequences not shorter than arg.
-carg, –cutoffarg
Remove similar sequences with #mismatches < cutoff.
-darg, –down-toarg
Remove similar sequences down to arg sequences.
–remove-crazyarg
Remove arg outlier sequences – defined as sequences that are missing too many conserved sites.
–conservedarg (=0.75)
Fraction of sequences that must contain a letter for it to be considered conserved.
Synopsis
alignment-thinalignment-file [OPTIONS]
