Mercurial > repos > ktnyt > gembassy
diff GEMBASSY-1.0.3/src/gcodoncompiler.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/gcodoncompiler.c Fri Jun 26 05:19:29 2015 -0400 @@ -0,0 +1,163 @@ +/****************************************************************************** +** @source gcodoncompiler +** +** Calculate various kinds of amino acid and codon usage data +** +** @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 gcodoncompiler ****************************************************** +** +** Calculate various kinds of amino acid and codon usage data +** +******************************************************************************/ + +int main(int argc, char *argv[]) +{ + embInitPV("gcodoncompiler", argc, argv, "GEMBASSY", "1.0.3"); + + AjPSeqall seqall; + AjPSeq seq; + AjPStr inseq = NULL; + + AjBool translate = ajFalse; + AjBool startcodon = ajFalse; + AjBool stopcodon = ajFalse; + AjPStr delkey = NULL; + AjPStr data = NULL; + + AjBool accid = ajFalse; + AjPStr restid = NULL; + AjPStr seqid = NULL; + + AjPStr base = NULL; + AjPStr url = NULL; + + AjPFile tmpfile = NULL; + AjPStr tmpname = NULL; + + AjPFile outf = NULL; + + seqall = ajAcdGetSeqall("sequence"); + translate = ajAcdGetBoolean("translate"); + startcodon = ajAcdGetBoolean("startcodon"); + stopcodon = ajAcdGetBoolean("stopcodon"); + delkey = ajAcdGetString("delkey"); + data = ajAcdGetListSingle("data"); + accid = ajAcdGetBoolean("accid"); + outf = ajAcdGetOutfile("outfile"); + + base = ajStrNewC("rest.g-language.org"); + + gAssignUniqueName(&tmpname); + + while(ajSeqallNext(seqall, &seq)) + { + inseq = NULL; + + if(!accid) + { + if(gFormatGenbank(seq, &inseq)) + { + tmpfile = ajFileNewOutNameS(tmpname); + if(!tmpfile) + { + ajDie("Output file (%S) open error\n", tmpname); + } + ajFmtPrintF(tmpfile, "%S", inseq); + ajFileClose(&tmpfile); + ajFmtPrintS(&url, "http://%S/upload/upl.pl", base); + gFilePostSS(url, tmpname, &restid); + ajStrDel(&url); + ajSysFileUnlinkS(tmpname); + } + else + { + ajWarn("Sequence does not have features\n" + "Proceeding with sequence accession ID\n"); + accid = ajTrue; + } + } + + ajStrAssignS(&seqid, ajSeqGetAccS(seq)); + + if(ajStrGetLen(seqid) == 0) + { + ajStrAssignS(&seqid, ajSeqGetNameS(seq)); + } + + if(ajStrGetLen(seqid) == 0) + { + ajWarn("No valid header information\n"); + } + + if(accid) + { + ajStrAssignS(&restid, seqid); + if(ajStrGetLen(seqid) == 0) + { + ajDie("Cannot proceed without header with -accid\n"); + } + + if(!gValID(seqid)) + { + ajDie("Invalid accession ID:%S, exiting\n", seqid); + } + } + + url = ajStrNew(); + + ajFmtPrintS(&url, "http://%S/%S/codon_compiler//translate=%d/" + "startcodon=%d/stopcodon=%d/delkey=%S/data=%S/output=f/", + base, restid, translate, startcodon, stopcodon, delkey, data); + + 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); + ajStrDel(&inseq); + } + + ajFileClose(&outf); + + ajSeqallDel(&seqall); + ajSeqDel(&seq); + ajStrDel(&base); + + ajStrDel(&delkey); + ajStrDel(&data); + + embExit(); + + return 0; +}