This default behaviour can be changed by using the qualifier '-ossingle' which writes many sequences into many files, each containing one sequence.
The program seqretsplit will take a file containing many sequences and will output many files, each containing one sequence. However you have no choice over the naming of the files - they are named after the ID name fo the sequence they contain.
If, however you have the situation where you have a file containing multiple sequences and you wish to extract one of them, then this application may be useful.
nthseq allows you to specify the name of the output file, so you may find that it is useful to include this program in scripts where you need to be able to specify the name of the resulting sequence files you create.
This application extracts the indicated sequence from a multiple set of sequences and writes it out.
|
In the example, the second sequence from the input file will be written out to the specified output file.
For example:
#!/usr/local/bin/perl -w if ($#ARGV !=1) { die "Usage: scriptname in out\n"; } $count=1; @list = `infoseq $ARGV[0] -auto -only -name`; while ($count <= $#list+1) { system("nthseq -auto $ARGV[0] -n $count $ARGV[1]-$count.seq"); $count++; }
The program seqretsplit will take a file containing many sequences and will output many files, each containing one sequence. However you have no choice over the naming of the files - they are named after the ID name fo the sequence they contain.