gp_getseq [-n name] [-q] [-v] [-d] [-h] start end [inputfile] [outputfile]
gp_getseq -r [-n name] [-q] [-v] [-d] [-h] [inputfile] [outputfile]
gp_getseq cuts out a specific fragment from a sequence or sequences. If the start position is greater then the end position, i.e. the sequence is on the minus strang, gp_getseq accordingly produces a reverse complement.
start denotes the number of the first base to be contained in the new sequence, and end the number of the last base which will be contained in the new sequence.
If the -r option is used, gp_getseq produces the reverse complement of the sequence(s) read.
Usually, gp_getseq takes the name of the original sequence and ads a comment containing the position of the newly produced fragment. If you supply it with the option -n name, name will be assigned to the sequence.
gp_getseq 1 100 ust.seq
All Genpak programs complain in situations you would also complain, like when they cannot find a sequence you gave them or the sequence is not valid.
The Genpak programs do not write over existing files. I have found this feature very useful :-)
I'm sure there are plenty left, so please mail me if you find them. I tried to clean up every bug I could find.
January Weiner III <january@bioinformatics.org>