]> Scans a sparse protein signature against swissprot (EMBOSS) protein:3d structure http://bioweb.pasteur.fr/docs/EMBOSS/sigscan.html sigscan &emboss_init; input input Section sigin Name of signature file for input (-sigin) perl " -sigin=$value" 1 database Name of sequence database to search (-database) perl " -database=$value -sformat=fasta" 2 8 seqsfile perl 1 targetf Name of (optionally grouped) scop families file for input (-targetf) perl " -targetf=$value" 3 required required Section thresh Minimum length (residues) of overlap required for two hits with the same code to be counted as the same hit. (-thresh) 20 perl " -thresh=$value" 4 sub Residue substitution matrix (-sub) EPAM60 EPAM60 EPAM290 EPAM290 EPAM470 EPAM470 EPAM110 EPAM110 EBLOSUM50 EBLOSUM50 EPAM220 EPAM220 EBLOSUM62-12 EBLOSUM62-12 EPAM400 EPAM400 EPAM150 EPAM150 EPAM330 EPAM330 EBLOSUM55 EBLOSUM55 EPAM30 EPAM30 EPAM260 EPAM260 EBLOSUM90 EBLOSUM90 EPAM440 EPAM440 EPAM190 EPAM190 EPAM370 EPAM370 EPAM70 EPAM70 EPAM480 EPAM480 EPAM120 EPAM120 EDNAMAT EDNAMAT EPAM300 EPAM300 EBLOSUM60 EBLOSUM60 EPAM230 EPAM230 EBLOSUM62 EBLOSUM62 EPAM410 EPAM410 EPAM160 EPAM160 EPAM340 EPAM340 EBLOSUM65 EBLOSUM65 EPAM40 EPAM40 EPAM270 EPAM270 EPAM450 EPAM450 EPAM380 EPAM380 EPAM80 EPAM80 EPAM490 EPAM490 EBLOSUM30 EBLOSUM30 EBLOSUMN EBLOSUMN EPAM200 EPAM200 EPAM130 EPAM130 EBLOSUM35 EBLOSUM35 EPAM310 EPAM310 EBLOSUM70 EBLOSUM70 EPAM10 EPAM10 EPAM240 EPAM240 EPAM420 EPAM420 EPAM170 EPAM170 EBLOSUM75 EBLOSUM75 EPAM350 EPAM350 EPAM280 EPAM280 EPAM50 EPAM50 EPAM460 EPAM460 EPAM390 EPAM390 EPAM90 EPAM90 EPAM100 EPAM100 EBLOSUM40 EBLOSUM40 EPAM210 EPAM210 EPAM140 EPAM140 EBLOSUM45 EBLOSUM45 EPAM320 EPAM320 EBLOSUM80 EBLOSUM80 EPAM500 EPAM500 EPAM20 EPAM20 EPAM250 EPAM250 EPAM430 EPAM430 EPAM180 EPAM180 EBLOSUM85 EBLOSUM85 EPAM360 EPAM360 ./EBLOSUM62 perl " -sub=$value" 5 gapo Gap insertion penalty (-gapo) 10 perl " -gapo=$value" 6 The gap insertion 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 EDNAMAT matrix for nucleotide sequences. Allowed values: Floating point number from 1.0 to 100.0 1. 100. gape Gap extension penalty (-gape) 0.5 perl " -gape=$value" 7 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. Allowed values: Floating point number from 0.0 to 10.0 0.0 10. nterm Select number -- N-terminal matching options (-nterm) 1 Align anywhere and allow only complete signature-sequence fit 2 Align anywhere and allow partial signature-sequence fit 3 Use empirical gaps only 1 perl " -nterm=$value" 8 output output Section nhits Number of hits to output (-nhits) 100 perl " -nhits=$value" 9 hitsf Name of signature hits file for output (-hitsf) test.hits perl " -hitsf=$value" 10 alignf Name of signature alignments file for output (-alignf) test.align perl " -alignf=$value" 11 auto perl " -auto -stdout" 12