logo
Free, unlimited AI code reviews that run on commit
git-lrc git-lrc GitHub Install Now We'd appreciate a star git-lrc - Free, unlimited AI code reviews that run on commit | Product Hunt git-lrc - Free, unlimited AI code reviews that run on commit | Product Hunt

skipredundant - Remove redundant sequences from an input set

Author

DebianMedPackagingTeam <debian-med-packaging@lists.alioth.debian.org>
           Wrote the script used to autogenerate this manual page.

Bugs

       Bugs can be reported to the Debian Bug Tracking system (http://bugs.debian.org/emboss), or directly to
       the EMBOSS developers (http://sourceforge.net/tracker/?group_id=93650&atid=605031).

Description

skipredundant is a command line program from EMBOSS (“the European Molecular Biology Open Software
       Suite”). It is part of the "Edit" command group(s).

Name

       skipredundant - Remove redundant sequences from an input set

Options

Inputsection-featuretoggle
           Sequence feature information will be retained if this option is set.

       -sequencesseqset-datafilematrixf
           This is the scoring matrix file used when comparing sequences. By default it is the file 'EBLOSUM62'
           (for proteins) or the file 'EDNAFULL' (for nucleic sequences). These files are found in the 'data'
           directory of the EMBOSS installation.

   Requiredsection-modelist
           This option specifies whether to remove redundancy at a single threshold percentage sequence
           similarity or remove redundancy outside a range of acceptable threshold percentage similarity. All
           permutations of pair-wise sequence alignments are calculated for each set of input sequences in turn
           using the EMBOSS implementation of the Needleman and Wunsch global alignment algorithm. Redundant
           sequences are removed in one of two modes as follows: (i) If a pair of proteins achieve greater than
           a threshold percentage sequence similarity (specified by the user) the shortest sequence is
           discarded. (ii) If a pair of proteins have a percentage sequence similarity that lies outside an
           acceptable range (specified by the user) the shortest sequence is discarded. Default value: 1

       -thresholdfloat
           This option specifies the percentage sequence identity redundancy threshold. The percentage sequence
           identity redundancy threshold determines the redundancy calculation. If a pair of proteins achieve
           greater than this threshold the shortest sequence is discarded. Default value: 95.0

       -minthresholdfloat
           This option specifies the percentage sequence identity redundancy threshold (lower limit). The
           percentage sequence identity redundancy threshold determines the redundancy calculation. If a pair of
           proteins have a percentage sequence similarity that lies outside an acceptable range the shortest
           sequence is discarded. Default value: 30.0

       -maxthresholdfloat
           This option specifies the percentage sequence identity redundancy threshold (upper limit). The
           percentage sequence identity redundancy threshold determines the redundancy calculation. If a pair of
           proteins have a percentage sequence similarity that lies outside an acceptable range the shortest
           sequence is discarded. Default value: 90.0

       -gapopenfloat
           The gap open penalty is the score taken away when a gap is created. The best value depends on the
           choice of comparison matrix. The default value assumes you are using the EBLOSUM62 matrix for protein
           sequences, and the EDNAFULL matrix for nucleotide sequences. Default value: @($(acdprotein)? 10.0 :
           10.0 )

       -gapextendfloat
           The gap extension, penalty is added to the standard gap penalty for each base or residue in the gap.
           This is how long gaps are penalized. Usually you will expect a few long gaps rather than many short
           gaps, so the gap extension penalty should be lower than the gap penalty. An exception is where one or
           both sequences are single reads with possible sequencing errors in which case you would expect many
           single base gaps. You can get this result by setting the gap open penalty to zero (or very low) and
           using the gap extension penalty to control gap scoring. Default value: @($(acdprotein)? 0.5 : 0.5 )

   AdvancedsectionOutputsection-outseqseqoutall-redundantoutseqseqoutall

See Also

       skipredundant is fully documented via the tfm(1) system.

Synopsis

skipredundant-featuretoggle-sequencesseqset [-datafilematrixf] -modelist-thresholdfloat-minthresholdfloat-maxthresholdfloat-gapopenfloat-gapextendfloat-outseqseqoutall-redundantoutseqseqoutallskipredundant-help

See Also