int main(int argc, char **argv)
{
AjPSeqset seqset = NULL;
AjPStr cl = NULL;
AjPSeqout seqout = NULL;
AjBool full = ajFalse;
AjPStr fn = NULL;
AjPStr stmp = NULL;
AjPStr outfname = NULL;
embInitPV("echlorop", argc, argv, "CBSTOOLS", VERSION);
seqset = ajAcdGetSeqset("sequence");
outfname= ajAcdGetOutfileName("outfile");
full = ajAcdGetBoolean("full");
cl = ajStrNewS(ajAcdGetpathC("chlorop"));
fn = ajStrNew();
stmp = ajStrNew();
ajFilenameSetTempname(&fn);
seqout = ajSeqoutNew();
if(!ajSeqoutOpenFilename(seqout, fn))
ajFatal("Cannot open temporary file %S",fn);
ajSeqoutSetFormatC(seqout, "fasta");
ajSeqoutWriteSet(seqout,seqset);
ajSeqoutClose(seqout);
if(full)
ajStrAppendC(&cl," -F");
ajFmtPrintS(&stmp," %S",fn);
ajStrAppendS(&cl,stmp);
#if 0
ajFmtPrint("%S\n",cl);
#endif
#if 1
ajSysExecOutnameAppendS(cl, outfname);
#endif
ajSysFileUnlinkS(fn);
ajStrDel(&cl);
ajStrDel(&fn);
ajStrDel(&stmp);
ajSeqoutDel(&seqout);
ajSeqsetDel(&seqset);
ajStrDel(&outfname);
embExit();
return 0;
}
int main(int argc, char **argv)
{
/* ACD data item variables */
AjPFile hmmfile = NULL;
AjPSeqall seqfile = NULL;
AjBool nuc = 0;
ajint A = 0;
float E = 0.0;
float T = 0.0;
ajint Z = 0;
AjBool acc = ajFalse;
AjBool compat = ajFalse;
ajint cpu = 0;
AjBool cutga = ajFalse;
AjBool cuttc = ajFalse;
AjBool cutnc = ajFalse;
float dome = 0.0;
float domt = 0.0;
AjBool forward = ajFalse;
AjBool nulltwo = ajFalse;
AjBool pvm = ajFalse;
AjBool xnu = ajFalse;
AjPStr outname = NULL;
/* Housekeeping variables */
AjPStr cmd = NULL;
AjPStr rnd = NULL;
AjPSeqout rndo = NULL;
AjPSeq seq = NULL;
/* ACD file processing */
embInitPV("ehmmpfam",argc,argv,"HMMERNEW",VERSION);
hmmfile = ajAcdGetInfile("hmmfile");
seqfile = ajAcdGetSeqall("seqfile");
nuc = ajAcdGetBoolean("nuc");
A = ajAcdGetInt("A");
E = ajAcdGetFloat("E");
T = ajAcdGetFloat("T");
Z = ajAcdGetInt("Z");
acc = ajAcdGetBoolean("acc");
compat = ajAcdGetBoolean("compat");
cpu = ajAcdGetInt("cpu");
cutga = ajAcdGetBoolean("cutga");
cuttc = ajAcdGetBoolean("cuttc");
cutnc = ajAcdGetBoolean("cutnc");
dome = ajAcdGetFloat("dome");
domt = ajAcdGetFloat("domt");
forward = ajAcdGetBoolean("forward");
nulltwo = ajAcdGetBoolean("nulltwo");
pvm = ajAcdGetBoolean("pvm");
xnu = ajAcdGetBoolean("xnu");
outname = ajAcdGetOutfileName("outfile");
/* MAIN APPLICATION CODE */
/* 1. Housekeeping */
cmd = ajStrNew();
rnd = ajStrNew();
/* 2. Re-write seqfile to a temporary file in a format (fasta) HMMER can understand.
We cannot just pass the name of seqfile to HMMER as the name provided might be a
USA which HMMER would not understand. */
ajFilenameSetTempname(&rnd);
rndo = ajSeqoutNew();
if(!ajSeqoutOpenFilename(rndo, rnd))
ajFatal("Terminal ajSeqFileNewOut failure. Email EMBOSS helpdesk!\n");
ajSeqoutSetFormatC(rndo, "fasta");
while(ajSeqallNext(seqfile, &seq))
ajSeqoutWriteSeq(rndo, seq);
ajSeqoutClose(rndo);
ajSeqoutDel(&rndo);
/* 2. Build hmmpfam command line */
/* Command line is built in this order:
i. Application name.
ii. HMMER 'options' (in order they appear in ACD file)
iii.HMMER 'options' (that don't appear in ACD file)
iv. HMMER & new parameters.
*/
ajStrAssignS(&cmd, ajAcdGetpathC("hmmpfam"));
if(nuc)
ajStrAppendC(&cmd, " -n ");
ajFmtPrintAppS(&cmd, " -A %d -E %f -T %f -Z %d", A, E, T, Z);
if(acc)
ajStrAppendC(&cmd, " --acc ");
if(compat)
ajStrAppendC(&cmd, " --compat ");
if(cpu)
ajFmtPrintAppS(&cmd, " --cpu %d ", cpu);
if(cutga)
ajStrAppendC(&cmd, " --cutga ");
if(cuttc)
ajStrAppendC(&cmd, " --cuttc ");
if(cutnc)
ajStrAppendC(&cmd, " --cutnc ");
ajFmtPrintAppS(&cmd, " --domE %f --domT %f ", dome, domt);
if(forward)
ajStrAppendC(&cmd, " --forward ");
if(nulltwo)
ajStrAppendC(&cmd, " --null2 ");
if(pvm)
ajStrAppendC(&cmd, " --pvm ");
if(xnu)
ajStrAppendC(&cmd, " --xnu ");
/* Note output redirected to outname.
rnd is the name of the rewritten seqfile.
MUST specify FASTA format explicitly. */
ajFmtPrintAppS(&cmd, " --informat FASTA %s %S",
ajFileGetNameC(hmmfile),
rnd);
/* 3. Close ACD files. */
ajFileClose(&hmmfile);
ajSeqallDel(&seqfile);
/* 4. Call hmmpfam. Use C system call instead of ajSystem
so that redirect in cmd works ok. */
ajFmtPrint("\n%S\n\n", cmd);
ajSysExecOutnameAppendS(cmd,outname);
/* 5. Exit cleanly */
ajSysFileUnlinkS(rnd);
ajStrDel(&cmd);
ajStrDel(&rnd);
ajStrDel(&outname);
embExit();
return 0;
}
int main(int argc, char **argv)
{
AjPSeqset seqset = NULL;
AjPStr cl = NULL;
AjPSeqout seqout = NULL;
float cutoff = 0.;
AjBool best = ajFalse;
AjBool gff = ajFalse;
AjBool two = ajFalse;
AjBool kinase = ajFalse;
AjPStr rsd = NULL;
AjBool addseq = ajFalse;
AjBool plot = ajFalse;
const AjPStr ofn = NULL;
AjPStr fn = NULL;
AjPStr stmp = NULL;
AjPFile outf = NULL;
embInitPV("enetphos", argc, argv, "CBSTOOLS", VERSION);
seqset = ajAcdGetSeqset("sequence");
outf = ajAcdGetOutfile("outfile");
plot = ajAcdGetBoolean("plot");
best = ajAcdGetBoolean("best");
gff = ajAcdGetBoolean("gff");
two = ajAcdGetBoolean("two");
kinase = ajAcdGetBoolean("kinase");
addseq = ajAcdGetBoolean("addseq");
cutoff = ajAcdGetFloat("cutoff");
rsd = ajAcdGetListSingle("residue");
cl = ajStrNewC("netphos ");
fn = ajStrNew();
stmp = ajStrNew();
ajFilenameSetTempname(&fn);
seqout = ajSeqoutNew();
if(!ajSeqoutOpenFilename(seqout, fn))
ajFatal("Cannot open temporary file %S",fn);
ajSeqoutSetFormatC(seqout, "fasta");
ajSeqoutWriteSet(seqout,seqset);
ajSeqoutClose(seqout);
if(ajStrMatchC(rsd,"serine"))
ajStrAppendC(&cl," -s");
else if(ajStrMatchC(rsd,"threonine"))
ajStrAppendC(&cl," -t");
else if(ajStrMatchC(rsd,"tyrosine"))
ajStrAppendC(&cl," -y");
if(plot)
ajStrAppendC(&cl," -g");
if(two)
ajStrAppendC(&cl," -2");
if(best)
ajStrAppendC(&cl," -b");
if(gff)
ajStrAppendC(&cl," -f gff");
if(kinase)
ajStrAppendC(&cl," -k");
if(addseq)
ajStrAppendC(&cl," -S");
ajFmtPrintS(&stmp," -c %f",cutoff);
ajStrAppendS(&cl,stmp);
ajFmtPrintS(&stmp," %S",fn);
ajStrAppendS(&cl,stmp);
ofn = ajFileGetNameS(outf);
ajFmtPrintS(&stmp," > %S",ofn);
ajStrAppendS(&cl,stmp);
ajFileClose(&outf);
#if 0
ajFmtPrint("%S\n",cl);
#endif
#if 1
if(system(ajStrGetPtr(cl)) == -1)
ajFatal("Command %S failed",cl);
#endif
ajSysFileUnlink(fn);
ajStrDel(&cl);
ajStrDel(&fn);
ajStrDel(&stmp);
ajStrDel(&rsd);
ajSeqoutDel(&seqout);
ajSeqsetDel(&seqset);
embExit();
return 0;
}
int main(int argc, char **argv)
{
/* ACD data item variables */
AjPSeqset dataset = NULL;
AjPFile bfile = NULL;
AjPFile plib = NULL;
AjPStr mod = NULL;
ajint nmotifs = 0;
AjBool text = ajFalse;
AjPStr prior = NULL;
float evt = 0.0;
ajint nsites = 0;
ajint minsites = 0;
ajint maxsites = 0;
float wnsites = 0.0;
ajint w = 0;
ajint minw = 0;
ajint maxw = 0;
AjBool nomatrim = ajFalse;
ajint wg = 0;
ajint ws = 0;
AjBool noendgaps = ajFalse;
AjBool revcomp = ajFalse;
AjBool pal = ajFalse;
AjBool nostatus = ajFalse;
ajint maxiter = 0;
float distance = 0.0;
float b = 0.0;
float spfuzz = 0.0;
AjPStr spmap = NULL;
AjPStr cons = NULL;
ajint maxsize = 0;
ajint p = 0;
ajint time = 0;
AjPStr sf = NULL;
ajint heapsize = 64;
AjBool xbranch = ajFalse;
AjBool wbranch = ajFalse;
ajint bfactor = 0;
AjPFile outtext = NULL;
/* Housekeeping variables */
AjPStr cmd = NULL;
AjPStr ssname = NULL;
AjPSeqout outseq = NULL;
AjPStr tmp = NULL;
char option;
/* ACD file processing */
embInitPV("ememetext",argc,argv,"MEME",VERSION);
dataset = ajAcdGetSeqset("dataset");
bfile = ajAcdGetInfile("bfile");
plib = ajAcdGetInfile("plibfile");
mod = ajAcdGetSelectSingle("mod");
nmotifs = ajAcdGetInt("nmotifs");
text = ajAcdGetBoolean("text");
prior = ajAcdGetSelectSingle("prior");
evt = ajAcdGetFloat("evt");
nsites = ajAcdGetInt("nsites");
minsites = ajAcdGetInt("minsites");
maxsites = ajAcdGetInt("maxsites");
wnsites = ajAcdGetFloat("wnsites");
w = ajAcdGetInt("w");
minw = ajAcdGetInt("minw");
maxw = ajAcdGetInt("maxw");
nomatrim = ajAcdGetBoolean("nomatrim");
wg = ajAcdGetInt("wg");
ws = ajAcdGetInt("ws");
noendgaps = ajAcdGetBoolean("noendgaps");
revcomp = ajAcdGetBoolean("revcomp");
pal = ajAcdGetBoolean("pal");
nostatus = ajAcdGetBoolean("nostatus");
maxiter = ajAcdGetInt("maxiter");
distance = ajAcdGetFloat("distance");
b = ajAcdGetFloat("b");
spfuzz = ajAcdGetFloat("spfuzz");
spmap = ajAcdGetSelectSingle("spmap");
cons = ajAcdGetString("cons");
maxsize = ajAcdGetInt("maxsize");
p = ajAcdGetInt("p");
time = ajAcdGetInt("time");
sf = ajAcdGetString("sf");
heapsize = ajAcdGetInt("heapsize");
xbranch = ajAcdGetBoolean("xbranch");
wbranch = ajAcdGetBoolean("wbranch");
bfactor = ajAcdGetInt("bfactor");
outtext = ajAcdGetOutfile("outtext");
outseq = ajAcdGetSeqoutset("outseq");
/* MAIN APPLICATION CODE */
/* 1. Housekeeping */
cmd = ajStrNew();
tmp = ajStrNew();
/* 2. Re-write dataset to a temporary file in a format (fasta) MEME
** can understand.
** Can't just pass the name of dataset to MEME as the name provided
** might be a USA which MEME would not understand.
*/
ssname = ajStrNewS(ajFileGetNameS(outseq->File));
ajSeqoutSetFormatC(outseq, "fasta");
ajSeqoutWriteSet(outseq, dataset);
ajSeqoutClose(outseq);
ajSeqoutDel(&outseq);
/* 3. Build ememe command line */
/* Command line is built in this order:
i. Application name.
ii. Original MEME options (in order they appear in ACD file)
iii.Original MEME options (that don't appear in ACD file)
iv. EMBASSY MEME new qualifiers and parameters.
*/
ajStrAssignS(&cmd, ajAcdGetpathC("meme"));
ajFmtPrintAppS(&cmd, " %S", ssname);
if(bfile)
ajFmtPrintAppS(&cmd, " -bfile %s ", ajFileGetNameC(bfile));
if(plib)
ajFmtPrintAppS(&cmd, " -plib %s ", ajFileGetNameC(plib));
option = ajStrGetCharFirst(mod);
if(option == 'o')
ajStrAppendC(&cmd, " -mod oops ");
else if(option == 'z')
ajStrAppendC(&cmd, " -mod zoops ");
else if(option == 'a')
ajStrAppendC(&cmd, " -mod anr ");
if(nmotifs != 1)
ajFmtPrintAppS(&cmd, " -nmotifs %d ", nmotifs);
if(text)
ajFmtPrintAppS(&cmd, " -text ");
ajFmtPrintAppS(&cmd, " -prior %S ", prior);
if(evt != -1)
ajFmtPrintAppS(&cmd, " -evt %f ", evt);
if(nsites != -1)
ajFmtPrintAppS(&cmd, " -nsites %d ", nsites);
else
{
if(minsites != -1)
ajFmtPrintAppS(&cmd, " -minsites %d ", minsites);
if(maxsites != -1)
ajFmtPrintAppS(&cmd, " -maxsites %d ", maxsites);
}
if(wnsites < 0.7999 || wnsites > .8001)
ajFmtPrintAppS(&cmd, " -wnsites %f ", wnsites);
if(w != -1)
ajFmtPrintAppS(&cmd, " -w %d ", w);
if(minw != 8)
ajFmtPrintAppS(&cmd, " -minw %d ", minw);
if(maxw != 50)
ajFmtPrintAppS(&cmd, " -maxw %d ", maxw);
if(nomatrim)
ajFmtPrintAppS(&cmd, " -nomatrim ");
if(wg != 11)
ajFmtPrintAppS(&cmd, " -wg %d ", wg);
if(ws != 1)
ajFmtPrintAppS(&cmd, " -ws %d ", ws);
if(noendgaps)
ajFmtPrintAppS(&cmd, " -noendgaps ");
if(revcomp)
ajFmtPrintAppS(&cmd, " -revcomp ");
if(pal && ajSeqsetIsNuc(dataset))
ajFmtPrintAppS(&cmd, " -pal ");
if(nostatus)
ajFmtPrintAppS(&cmd, " -nostatus ");
if(maxiter != 50)
ajFmtPrintAppS(&cmd, " -maxiter %d ", maxiter);
if(distance < 0.00099 || distance > 0.00101)
ajFmtPrintAppS(&cmd, " -distance %f ", distance);
if(b != -1)
ajFmtPrintAppS(&cmd, " -b %f ", b);
if(spfuzz != -1)
ajFmtPrintAppS(&cmd, " -spfuzz %f ", spfuzz);
if(!ajStrMatchC(spmap,"default"))
ajFmtPrintAppS(&cmd, " -spmap %S ", spmap);
if(MAJSTRGETLEN(cons))
ajFmtPrintAppS(&cmd, "-cons %S", cons);
if(maxsize != -1)
ajFmtPrintAppS(&cmd, " -maxsize %d ", maxsize);
if(p > 0)
ajFmtPrintAppS(&cmd, " -p %d ", p);
if(time > 0)
ajFmtPrintAppS(&cmd, " -time %d ", time);
if(MAJSTRGETLEN(sf))
ajFmtPrintAppS(&cmd, " -sf %S", sf);
if(heapsize != 64)
ajFmtPrintAppS(&cmd, " -heapsize %d ", heapsize);
if(xbranch)
ajFmtPrintAppS(&cmd, " -x_branch");
if(wbranch)
ajFmtPrintAppS(&cmd, " -w_branch");
if(bfactor != 3)
ajFmtPrintAppS(&cmd, " -bfactor %d ", bfactor);
if(ajSeqsetIsProt(dataset))
ajFmtPrintAppS(&cmd, "-protein ");
else
ajFmtPrintAppS(&cmd, "-dna ");
ajFmtPrintAppS(&cmd, " -text");
ajFmtPrintAppS(&cmd, " > %S ", ajFileGetNameS(outtext));
/* 4. Close files from ACD before calling meme */
ajFileClose(&bfile);
ajFileClose(&plib);
/* 5. Call meme */
/* ajFmtPrint("\n%S\n", cmd); */
system(ajStrGetPtr(cmd));
/* 6. Exit cleanly */
ajSeqsetDel(&dataset);
ajStrDel(&cons);
ajStrDel(&sf);
ajStrDel(&mod);
ajStrDel(&prior);
ajStrDel(&spmap);
ajStrDel(&cmd);
ajStrDel(&ssname);
ajStrDel(&tmp);
ajFileClose(&bfile);
ajFileClose(&plib);
ajFileClose(&outtext);
ajSeqoutDel(&outseq);
embExit();
return 0;
}
int main(int argc, char **argv)
{
/* ACD data item variables */
AjPFile hmmfile = NULL;
AjPSeqset seqfile = NULL;
AjPFile mapali = NULL;
AjPFile withali = NULL;
AjPAlign o = NULL;
AjBool m = ajFalse;
AjBool q = ajFalse;
/* Housekeeping variables */
AjPStr cmd = NULL;
AjPStr fmt = NULL;
AjBool fmtok = ajFalse;
AjPStr rnd = NULL;
AjPSeqout rndo = NULL;
/* ACD file processing */
embInitPV("ehmmalign",argc,argv,"HMMERNEW",VERSION);
hmmfile = ajAcdGetInfile("hmmfile");
seqfile = ajAcdGetSeqset("seqfile");
mapali = ajAcdGetInfile("mapali");
withali = ajAcdGetInfile("withali");
o = ajAcdGetAlign("o");
m = ajAcdGetBoolean("m");
q = ajAcdGetBoolean("q");
/* MAIN APPLICATION CODE */
/* 1. Housekeeping */
cmd = ajStrNew();
fmt = ajStrNew();
rnd = ajStrNew();
/* 2. Re-write seqfile to a temporary file in a format (fasta) HMMER can understand.
We cannot just pass the name of seqfile to HMMER as the name provided might be a
USA which HMMER would not understand. */
ajFilenameSetTempname(&rnd);
rndo = ajSeqoutNew();
if(!ajSeqoutOpenFilename(rndo, rnd))
ajFatal("Terminal ajSeqFileNewOut failure. Email EMBOSS helpdesk!\n");
ajSeqoutSetFormatC(rndo, "fasta");
ajSeqoutWriteSet(rndo, seqfile);
ajSeqoutClose(rndo);
ajSeqoutDel(&rndo);
/* 3. Build hmmalign command line */
/* Command line is built in this order:
i. Application name.
ii. HMMER 'options' (in order they appear in ACD file)
iii.HMMER 'options' (that don't appear in ACD file)
iv. HMMER & new parameters.
*/
ajFmtPrintS(&cmd, "%S ", ajAcdGetpathC("hmmalign"));
if(mapali)
ajFmtPrintAppS(&cmd, " --mapali %s ", ajFileGetNameC(mapali));
if(withali)
ajFmtPrintAppS(&cmd, " --withali %s ", ajFileGetNameC(withali));
if(m)
ajStrAppendC(&cmd, " -m ");
if(q)
ajStrAppendC(&cmd, " -q ");
/* Ensure output alignment is in user-specified format. */
fmtok=ajTrue;
ajStrAssignS(&fmt, ajAlignGetFormat(o));
/* fasta and a2m are identical formats. */
if(ajStrMatchC(fmt, "fasta"))
ajStrAssignC(&fmt, "A2M");
else if(ajStrMatchC(fmt, "a2m"))
ajStrAssignC(&fmt, "A2M");
else if(ajStrMatchC(fmt, "msf"))
ajStrAssignC(&fmt, "MSF");
else if(ajStrMatchC(fmt, "phylip"))
ajStrAssignC(&fmt, "PHYLIP");
/* hmmer also supports stockholm, SELEX & Clustal output, EMBOSS does not.
EMBOSS supports unknown/multiple/simple and srs output, hmmer does not. */
else
fmtok = ajFalse;
if(!fmtok)
{
/* This could be replaced with code to reformat the file. */
ajWarn("Specified output alignment format ('o' ACD option) is "
"not understood by HMMER. Using stockholm format instead.");
ajStrAssignC(&fmt, "Stockholm");
}
/* rnd is the name of the rewritten seqfile. MUST specify FASTA format explicitly. */
ajFmtPrintAppS(&cmd, " --informat FASTA --outformat %S -o %s %s %S",
fmt,
ajAlignGetFilename(o),
ajFileGetNameC(hmmfile),
rnd);
/* 4. Close ACD files */
ajFileClose(&hmmfile);
ajSeqsetDel(&seqfile);
ajFileClose(&mapali);
ajFileClose(&withali);
ajAlignClose(o);
ajAlignDel(&o);
/* 5. Call hmmalign */
ajFmtPrint("\n%S\n\n", cmd);
ajSysExecS(cmd);
/* 6. Exit cleanly */
ajSysFileUnlinkS(rnd);
ajStrDel(&cmd);
ajStrDel(&fmt);
ajStrDel(&rnd);
embExit();
return 0;
}
int main(int argc, char **argv)
{
AjPSeqset seqset = NULL;
AjPStr cl = NULL;
AjPSeqout seqout = NULL;
AjBool pc = ajFalse;
AjBool sigp = ajFalse;
AjBool plot = ajFalse;
const AjPStr ofn = NULL;
AjPStr fn = NULL;
AjPStr stmp = NULL;
AjPFile outf = NULL;
embInitPV("eprop", argc, argv, "CBSTOOLS", VERSION);
seqset = ajAcdGetSeqset("sequence");
outf = ajAcdGetOutfile("outfile");
plot = ajAcdGetBoolean("plot");
pc = ajAcdGetBoolean("pcprediction");
sigp = ajAcdGetBoolean("signalp");
cl = ajStrNewC("prop ");
fn = ajStrNew();
stmp = ajStrNew();
ajFilenameSetTempname(&fn);
seqout = ajSeqoutNew();
if(!ajSeqoutOpenFilename(seqout, fn))
ajFatal("Cannot open temporary file %S",fn);
ajSeqoutSetFormatC(seqout, "fasta");
ajSeqoutWriteSet(seqout,seqset);
ajSeqoutClose(seqout);
if(plot)
ajStrAppendC(&cl," -g");
if(pc)
ajStrAppendC(&cl," -p");
if(sigp)
ajStrAppendC(&cl," -s");
ajFmtPrintS(&stmp," %S",fn);
ajStrAppendS(&cl,stmp);
ofn = ajFileGetNameS(outf);
ajFmtPrintS(&stmp," > %S",ofn);
ajStrAppendS(&cl,stmp);
ajFileClose(&outf);
#if 0
ajFmtPrint("%S\n",cl);
#endif
#if 1
if(system(ajStrGetPtr(cl)) == -1)
ajFatal("Command %S failed",cl);
#endif
ajSysFileUnlink(fn);
ajStrDel(&cl);
ajStrDel(&fn);
ajStrDel(&stmp);
ajSeqoutDel(&seqout);
ajSeqsetDel(&seqset);
embExit();
return 0;
}
int main(int argc, char **argv)
{
AjPStr cl = NULL;
AjPSeqset seqset = NULL;
AjPSeqout seqout = NULL;
AjPStr stmp = NULL;
AjPStr email = NULL;
AjPStr fmt = NULL;
AjPStr trtab = NULL;
const AjPStr ofn = NULL;
AjPStr *applist = NULL;
AjBool crc = ajFalse;
AjBool alt = ajFalse;
AjBool iprlook = ajFalse;
AjBool goterms = ajFalse;
ajint trlen = 0;
ajuint i = 0;
ajuint lcnt = 0;
AjPStr fn = NULL;
AjPFile outf = NULL;
embInitPV("eiprscan", argc, argv, "IPRSCAN", VERSION);
seqset = ajAcdGetSeqset("sequence");
email = ajAcdGetString("email");
crc = ajAcdGetBoolean("crc");
applist = ajAcdGetList("appl");
fmt = ajAcdGetListSingle("format");
trtab = ajAcdGetListSingle("trtable");
trlen = ajAcdGetInt("trlen");
alt = ajAcdGetBoolean("altjobs");
iprlook = ajAcdGetBoolean("iprlookup");
goterms = ajAcdGetBoolean("goterms");
outf = ajAcdGetOutfile("outfile");
stmp = ajStrNew();
cl = ajStrNewC("iprscan -cli");
fn = ajStrNew();
ajFilenameSetTempname(&fn);
seqout = ajSeqoutNew();
if(!ajSeqoutOpenFilename(seqout, fn))
ajFatal("Cannot open temporary file %S",fn);
ajSeqoutSetFormatC(seqout, "fasta");
ajSeqoutWriteSet(seqout,seqset);
ajSeqoutClose(seqout);
ajFmtPrintS(&stmp," -i %S",fn);
ajStrAppendS(&cl,stmp);
if(!ajSeqsetIsProt(seqset))
ajStrAppendC(&cl," -seqtype n");
if(!crc)
ajStrAppendC(&cl," -nocrc");
if(ajStrGetLen(email))
{
ajFmtPrintS(&stmp," -email %S",email);
ajStrAppendS(&cl,stmp);
}
i = 0;
lcnt = 0;
while(applist[i])
{
++lcnt;
++i;
}
i = 0;
while(applist[i])
{
if(ajStrMatchC(applist[i],"all"))
{
if(lcnt != 1)
ajFatal("Cannot specify 'all' with multiple "
"applications");
++i;
continue;
}
ajFmtPrintS(&stmp," -appl %S",applist[i]);
ajStrAppendS(&cl,stmp);
++i;
}
ajFmtPrintS(&stmp," -format %S",fmt);
ajStrAppendS(&cl,stmp);
ajFmtPrintS(&stmp," -trtable %S",trtab);
ajStrAppendS(&cl,stmp);
ajFmtPrintS(&stmp," -trlen %d",trlen);
ajStrAppendS(&cl,stmp);
if(alt)
ajStrAppendC(&cl," -altjobs");
if(iprlook)
ajStrAppendC(&cl," -iprlookup");
if(goterms)
ajStrAppendC(&cl," -goterms");
ofn = ajFileGetNameS(outf);
ajFmtPrintS(&stmp," -o %S",ofn);
ajFileClose(&outf);
ajStrAppendS(&cl,stmp);
#if 0
ajFmtPrint("%S\n",cl);
#endif
#if 1
system(ajStrGetPtr(cl));
#endif
ajSysFileUnlink(fn);
ajStrDel(&cl);
ajStrDel(&fn);
ajStrDel(&stmp);
ajStrDel(&email);
ajStrDel(&fmt);
ajStrDel(&trtab);
ajStrDelarray(&applist);
ajSeqoutDel(&seqout);
ajSeqsetDel(&seqset);
embExit();
return 0;
}
int main(int argc, char **argv)
{
AjPSeqset seqset = NULL;
AjPStr cl = NULL;
AjPSeqout seqout = NULL;
float thresh = 0.;
AjBool netphos = ajFalse;
AjPStr format = NULL;
AjBool plot = ajFalse;
AjPStr fn = NULL;
AjPStr stmp = NULL;
AjPStr outfname = NULL;
embInitPV("eyinoyang", argc, argv, "CBSTOOLS",VERSION);
seqset = ajAcdGetSeqset("sequence");
outfname= ajAcdGetOutfileName("outfile");
plot = ajAcdGetBoolean("plot");
netphos = ajAcdGetBoolean("netphos");
thresh = ajAcdGetFloat("threshold");
format = ajAcdGetListSingle("format");
cl = ajStrNewS(ajAcdGetpathC("yinOyang"));
fn = ajStrNew();
stmp = ajStrNew();
ajFilenameSetTempname(&fn);
seqout = ajSeqoutNew();
if(!ajSeqoutOpenFilename(seqout, fn))
ajFatal("Cannot open temporary file %S",fn);
ajSeqoutSetFormatC(seqout, "fasta");
ajSeqoutWriteSet(seqout,seqset);
ajSeqoutClose(seqout);
if(ajStrMatchC(format,"short"))
ajStrAppendC(&cl," -f s");
else if(ajStrMatchC(format,"long"))
ajStrAppendC(&cl," -f l");
if(plot)
ajStrAppendC(&cl," -g");
if(netphos)
ajStrAppendC(&cl," -y");
if(netphos)
{
ajFmtPrintS(&stmp," -t %f",thresh);
ajStrAppendS(&cl,stmp);
}
ajFmtPrintS(&stmp," %S",fn);
ajStrAppendS(&cl,stmp);
#if 0
ajFmtPrint("%`S\n",cl);
#endif
#if 1
ajSysExecOutnameAppendS(cl, outfname);
#endif
ajSysFileUnlinkS(fn);
ajStrDel(&cl);
ajStrDel(&fn);
ajStrDel(&stmp);
ajStrDel(&format);
ajSeqoutDel(&seqout);
ajSeqsetDel(&seqset);
ajStrDel(&outfname);
embExit();
return 0;
}