extracting sequences from a Blast output seqsblast seqsblast seqlab seqsblast perl "seqsblast" 0 infile Blast output File perl " $value" 10 blast_output perl 1 howmany How many sequences? (-n) perl ($value) ? " -n $value" : "" 1 You must either enter a number of sequences or a significance threshold. perl ! ($signif || $howmany) You must either enter a number of sequences or a significance threshold (see: Other options). other_options Other options 1 outformat Output Sequence Format (readseq formats) (-f) perl " -f $value" 8 1 1 1. IG/Stanford 2 2. GenBank/GB 3 3. NBRF 4 4. EMBL 6 6. DNAStrider 7 7. Fitch 8 8. Pearson/Fasta 9 9. Zuker 10 10. Olsen (in-only) 14 14. PIR/CODATA 16 16. ASN.1 signif Signifiance threshold (-s) perl (defined $value)? " -s $value" : "" 1 print_hsp Get HSP (instead of entry) (-H) perl ($value)? " -H" : "" 1 1 extend_hsp Extend HSP to query (-e) perl ($value)? " -e" : "" 1 extend_left Extend to # additional position on the left (-l) perl ($value)? " -l $value" : "" 1 extend_right Extend to # additional position on the right (-r) perl ($value)? " -r $value" : "" 1 seqsfile *.seqs seqsfile perl 1