1 /* @source seqretset application
2 **
3 ** Read and write sequences as a set
4 **
5 ** @author Copyright (C) Peter Rice
6 ** @@
7 **
8 ** This program is free software; you can redistribute it and/or
9 ** modify it under the terms of the GNU General Public License
10 ** as published by the Free Software Foundation; either version 2
11 ** of the License, or (at your option) any later version.
12 **
13 ** This program is distributed in the hope that it will be useful,
14 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
15 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 ** GNU General Public License for more details.
17 **
18 ** You should have received a copy of the GNU General Public License
19 ** along with this program; if not, write to the Free Software
20 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
21 ******************************************************************************/
22
23 #include "emboss.h"
24
25
26
27
28 /* @prog nthseqset *********************************************************
29 **
30 ** Reads and writes (returns) one set from input containing multiple sets
31 ** of sequences
32 **
33 ******************************************************************************/
34
main(int argc,char ** argv)35 int main(int argc, char **argv)
36 {
37 AjPSeqset* seqsets;
38 AjPSeqout seqout;
39 ajint n;
40 ajint count = 0;
41 ajuint i;
42
43 embInit("nthseqset", argc, argv);
44
45 seqsets = ajAcdGetSeqsetall("sequence");
46 n = ajAcdGetInt("number");
47 seqout = ajAcdGetSeqoutall("outseq");
48
49 while(seqsets[count])
50 count++;
51
52 if(count < n)
53 {
54 ajSeqoutClose(seqout);
55 ajFatal("No such sequence set - only %d sequence sets were input.",
56 count);
57 embExitBad();
58 }
59
60 ajSeqoutWriteSet(seqout, seqsets[n-1]);
61 ajSeqoutClose(seqout);
62
63 for(i=0;seqsets[i];i++)
64 ajSeqsetDel(&seqsets[i]);
65 AJFREE(seqsets);
66 ajSeqoutDel(&seqout);
67
68 embExit();
69
70 return 0;
71 }
72