Mercurial > repos > ktnyt > gembassy
diff GEMBASSY-1.0.3/src/gaminoinfo.c @ 0:8300eb051bea draft
Initial upload
author | ktnyt |
---|---|
date | Fri, 26 Jun 2015 05:19:29 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/GEMBASSY-1.0.3/src/gaminoinfo.c Fri Jun 26 05:19:29 2015 -0400 @@ -0,0 +1,121 @@ +/****************************************************************************** +** @source gaminoinfo +** +** Prints out basic amino acid sequence statistics +** +** @author Copyright (C) 2012 Hidetoshi Itaya +** @version 1.0.3 +** @modified 2012/1/20 Hidetoshi Itaya Created! +** @modified 2013/6/16 Revision 1 +** @modified 2015/2/7 Refactor +** @@ +** +** This program is free software; you can redistribute it and/or +** modify it under the terms of the GNU General Public License +** as published by the Free Software Foundation; either version 2 +** of the License, or (at your option) any later version. +** +** This program is distributed in the hope that it will be useful, +** but WITHOUT ANY WARRANTY; without even the implied warranty of +** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +** GNU General Public License for more details. +** +** You should have received a copy of the GNU General Public License +** along with this program; if not, write to the Free Software +** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. +******************************************************************************/ + +#include "emboss.h" +#include "glibs.h" + + + + +/* @prog gaminoinfo *********************************************************** +** +** Prints out basic amino acid sequence statisctics +** +******************************************************************************/ + +int main(int argc, char *argv[]) +{ + embInitPV("gaminoinfo", argc, argv, "GEMBASSY", "1.0.3"); + + AjPSeqall seqall; + AjPSeq seq; + + AjPStr restid = NULL; + AjPStr seqid = NULL; + + AjPStr base = NULL; + AjPStr url = NULL; + + AjPStr tmpname = NULL; + AjPSeqout tmpout = NULL; + + AjPFile outf = NULL; + + seqall = ajAcdGetSeqall("sequence"); + outf = ajAcdGetOutfile("outfile"); + + base = ajStrNewC("rest.g-language.org"); + + gAssignUniqueName(&tmpname); + ajStrAppendC(&tmpname, ".fasta"); + + while(ajSeqallNext(seqall, &seq)) + { + tmpout = ajSeqoutNew(); + + if(!ajSeqoutOpenFilename(tmpout, tmpname)) + { + embExitBad(); + } + + ajSeqoutSetFormatS(tmpout,ajStrNewC("fasta")); + ajSeqoutWriteSeq(tmpout, seq); + ajSeqoutClose(tmpout); + ajSeqoutDel(&tmpout); + + ajFmtPrintS(&url, "http://%S/upload/upl.pl", base); + gFilePostSS(url, tmpname, &restid); + ajStrDel(&url); + ajSysFileUnlinkS(tmpname); + + ajStrAssignS(&seqid, ajSeqGetAccS(seq)); + + if(ajStrGetLen(seqid) == 0) + { + ajStrAssignS(&seqid, ajSeqGetNameS(seq)); + } + + if(ajStrGetLen(seqid) == 0) + { + ajWarn("No valid header information\n"); + } + + url = ajStrNew(); + + ajFmtPrintS(&url, "http://%S/%S/amino_info", base, restid); + + ajFmtPrintF(outf, "Sequence: %S\n", seqid); + if(!gFileOutURLS(url, &outf)) + { + ajDie("Failed to download result from:\n%S\n", url); + } + + ajStrDel(&url); + ajStrDel(&restid); + ajStrDel(&seqid); + } + + ajFileClose(&outf); + + ajSeqallDel(&seqall); + ajSeqDel(&seq); + ajStrDel(&base); + + embExit(); + + return 0; +}