ff9313091f342ffe62a124427b5c9d6a516faaea hiram Thu Oct 15 15:03:07 2015 -0700 fixup gcc warnings for -Wunused-but-set-variable refs #16121 diff --git src/hg/hgc/lowelab.c src/hg/hgc/lowelab.c index f865e17..a367911 100644 --- src/hg/hgc/lowelab.c +++ src/hg/hgc/lowelab.c @@ -989,162 +989,151 @@ hFreeConn(&conn); } /*Code to display Sargasso Sea Information*/ void doSargassoSea(struct trackDb *tdb, char *trnaName) { struct bed *cb=NULL; struct sargassoSeaXra *cbs=NULL, *cbs2, *list=NULL; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; char query[512]; char *dupe, *words[16]; char **row; char tempstring[255]=""; -int flag, z, dashes, wordCount, rowOffset; -int start = cartInt(cart, "o"), num = 0, flag2=0; +int z, dashes, wordCount; +int start = cartInt(cart, "o"), num = 0; float sequenceLength, dashlength=60; genericHeader(tdb,trnaName); dupe = cloneString(tdb->type); wordCount = chopLine(dupe, words); if (wordCount > 1) num = atoi(words[1]); if (num < 3) num = 3; genericBedClick(conn, tdb, trnaName, start, num); -rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, trnaName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) cb=bedLoadN(row+1, 6); sequenceLength=(cb->chromEnd - cb->chromStart); if(sequenceLength<0){ sequenceLength=sequenceLength*-1;} sequenceLength=sequenceLength/3; dashlength=sequenceLength/60; /*Query the database for the extrainfo file for sargassoSea*/ conn=hAllocConn(database);/*sqlConnect(dupe);*/ sqlSafef(tempstring, sizeof(tempstring),"select * from sargassoSeaXra where qname = '%s'", trnaName); sr = sqlGetResult(conn, tempstring); /*Load the required data from the database*/ while ((row = sqlNextRow(sr)) != NULL) { cbs=sargassoSeaXraLoad(row); slAddHead(&list, cbs); } slReverse(&list); -flag=0; -flag2=0; - /*Print out table with Blast information*/ printf(" \n\n

"); printf(""); printf(" \n \n"); printf(" \n \n \n \n \n \n \n \n \n "); -flag=0; for(cbs2=list;cbs2!=NULL;cbs2=cbs2->next) { printf("\n \n \n \n \n \n \n \n \n \n \n"); - flag=0; } /* printf("

\n
Organism"); printf("
\n
"); printf("Nucleotides \naligned begin"); printf("
\n
"); printf("Nucelotides \naligned end"); printf("
\n
NCBI Link"); printf("
\n
"); printf("Evalue"); printf("
\n
"); printf("Percent Identity"); printf("
\n
"); printf("Alignment Length"); printf("
\n
"); printf("Gap openings"); printf("
"); printf("",cbs2->GI); printf("%s",cbs2->GI,cbs2->species); printf("
\n
"); printf("%i",cbs2->thisseqstart); printf("
\n
"); printf("%i",cbs2->thisseqend); printf("
\n
"); printf("NCBI Link",cbs2->GI); printf("
\n
"); printf("%s",cbs2->evalue); printf("
\n
"); printf("%.2f",cbs2->PI); printf("
\n
"); printf("%i",cbs2->length); printf("
\n
"); printf("%i",cbs2->gap); printf("
\n"); */ /*Print out the table for the alignments*/ printf(""); printf(" \n \n"); printf(" \n \n \n "); printf(" \n \n \n "); -flag=0; -flag2=0; for(cbs2=list;cbs2!=NULL;cbs2=cbs2->next) { printf("\n \n \n \n \n"); - flag=0; - - } /*Free the data*/ hFreeConn(&conn); sargassoSeaXraFree(&cbs); printTrackHtml(tdb); } /*Function to print out full code name from code letter*/ void printCode(char code) { switch(code) @@ -1197,40 +1186,36 @@ case 'r': printf("Facultative Aerobe"); break; default: break; } } void doTrnaGenes(struct trackDb *tdb, char *trnaName) { struct tRNAs *trna; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; -char *dupe, *words[16]; char **row; -int wordCount; int rowOffset; char* chrom = cartString(cart, "c"); genericHeader(tdb,trnaName); -dupe = cloneString(tdb->type); -wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where chrom = '%s' and name = '%s'", tdb->table, chrom, trnaName); sr = sqlGetResult(conn, query); printf("
Organism"); printf("
\n
Alignment"); printf("
%s where each plus(+) is approx. %.1f amino acids", trnaName, dashlength); printf("
\n
"); printf("\n"); for(z=0; z<60; z++) { printf("+"); } printf(""); printf("
"); dashes=cbs2->queryseqstart-cbs2->queryseqend; if(dashes<0) dashes=dashes*-1; printf("",cbs2->GI); printf("%s, %s",cbs2->GI, cbs2->species,cbs2->qName); printf("\n "); printf("\n"); dashes=dashes/dashlength; if(cbs2->queryseqstart>cbs2->queryseqend) for(z=0; z<((cbs2->queryseqend)/dashlength); z++) { printf(" "); } else for(z=0; z<((cbs2->queryseqstart)/dashlength); z++) { printf(" "); } if(dashes<1) printf("+"); for(z=0; z"); printf("
\n"); while ((row = sqlNextRow(sr)) != NULL) { printf("\n"); printf(""); } printf("
\n"); trna = tRNAsLoad(row+rowOffset); printf("tRNA name: %s
\n",trna->name); printf("tRNA Isotype: %s
\n",trna->aa); printf("tRNA anticodon: %s
\n",trna->ac); printf("tRNAscan-SE score: %.2f
\n",trna->trnaScore); @@ -1265,38 +1250,34 @@ printf("
"); sqlFreeResult(&sr); hFreeConn(&conn); printTrackHtml(tdb); tRNAsFree(&trna); } void doSnornaGenes(struct trackDb *tdb, char *snornaName) { struct snoRNAs *snorna; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; -char *dupe, *words[16]; char **row; -int wordCount; int rowOffset; genericHeader(tdb,snornaName); -dupe = cloneString(tdb->type); -wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, snornaName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { snorna = snoRNAsLoad(row+rowOffset); printf("sRNA name: %s
\n",snorna->name); printf("Snoscan score: %.2f
\n",snorna->snoScore); printf("HMM snoRNA score: %.2f
\n",snorna->hmmScore); printf("Predicted targets: %s
\n",snorna->targetList); printf("Predicted guide interactions:
%s
\n",snorna->guideStr); printf("Possible sRNA homolog(s): %s
\n",snorna->orthologs); @@ -1310,39 +1291,34 @@ if (snorna->next != NULL) printf("
\n"); } sqlFreeResult(&sr); hFreeConn(&conn); printTrackHtml(tdb); snoRNAsFree(&snorna); } void doGbRnaGenes(struct trackDb *tdb, char *gbRnaName) { struct gbRNAs *gbRna; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; -char *dupe, *words[16]; char **row; -int wordCount; int rowOffset; genericHeader(tdb,gbRnaName); -dupe = cloneString(tdb->type); -wordCount = chopLine(dupe, words); - rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, gbRnaName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { gbRna = gbRNAsLoad(row+rowOffset); printf("Genbank ncRNA name: %s
\n",gbRna->name); printf("Product Description/Note: %s
\n",gbRna->product); printf ("Intron(s): %s
\n",gbRna->intron); @@ -1408,42 +1384,41 @@ printTrackHtml(tdb); easyGeneFreeList(&egList); } void doCodeBlast(struct trackDb *tdb, char *trnaName) { struct pepPred *pp=NULL; struct codeBlast *cb=NULL; struct codeBlastScore *cbs=NULL, *cbs2, *list=NULL; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; char query[512]; char *dupe, *words[16]; char **row; char tempstring[255]=""; -int flag, z, dashes, wordCount, rowOffset, currentGI=0; -int start = cartInt(cart, "o"), num = 0, flag2=0; +int flag, z, dashes, wordCount, currentGI=0; +int start = cartInt(cart, "o"), num = 0; float sequenceLength, dashlength=60; genericHeader(tdb,trnaName); dupe = cloneString(tdb->type); wordCount = chopLine(dupe, words); if (wordCount > 1) num = atoi(words[1]); if (num < 3) num = 3; genericBedClick(conn, tdb, trnaName, start, num); -rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, trnaName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { cb=codeBlastLoad(row); } sequenceLength=(cb->chromEnd - cb->chromStart); if(sequenceLength<0){ sequenceLength=sequenceLength*-1;} sequenceLength=sequenceLength/3; dashlength=sequenceLength/60; conn=hAllocConn(database);/*sqlConnect(dupe);*/ sqlSafef(query, sizeof query, "select * from gbProtCodePep where name = '%s'", trnaName); @@ -1669,31 +1644,30 @@ printf(" \n \n"); printf(" Organism"); printf("
\n \n Alignment"); printf("\n \n "); printf(" %s where each plus(+) is approx. %f amino acids", trnaName, dashlength); printf("
\n \n "); printf("\n"); for(z=0; z<60; z++) { printf("+"); } printf(""); printf("\n \n "); flag=0; -flag2=0; for(cbs2=list;cbs2!=NULL;cbs2=cbs2->next) { if(sameString(trnaName, cbs2->qName)) { if(flag==0) { currentGI=cbs2->GI; printf("\n \n "); flag=1; } if((cbs2->next!=NULL) && (currentGI== cbs2->GI) && (currentGI== cbs2->next->GI) ) { } else @@ -1883,41 +1857,34 @@ printTrackHtml(tdb);*/ /* clean up */ /*sqlFreeResult(&sr); hFreeConn(&conn); tigrOperonFree(&op); } */ void doTigrCmrGene(struct trackDb *tdb, char *tigrName) /* Handle the TIRG CMR gene track. */ { struct tigrCmrGene *tigr; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; - char *dupe, *words[16]; char **row; - int wordCount; - int rowOffset; - /* int start = cartInt(cart, "o"), num = 0; */ genericHeader(tdb,tigrName); - dupe = cloneString(tdb->type); - wordCount = chopLine(dupe, words); - rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, tigrName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { tigr = tigrCmrGeneLoad(row); if (tigr != NULL) { printf("TIGR locus name: %s
\n",tigrName); printf("TIGR gene description: %s
\n",tigr->tigrCommon); printf("Alternate TIGR gene name: "); if (strlen(tigr->tigrGene) >0) { printf("%s
\n",tigr->tigrGene); } else { printf("None
"); @@ -1948,38 +1915,34 @@ } } sqlFreeResult(&sr); hFreeConn(&conn); printTrackHtml(tdb); tigrCmrGeneFree(&tigr); } void doJgiGene(struct trackDb *tdb, char *jgiName) /* Handle the JGI gene track. */ { struct jgiGene *jgi; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; - char *dupe, *words[16]; char **row; - int wordCount; int rowOffset; genericHeader(tdb,jgiName); - dupe = cloneString(tdb->type); - wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s'", tdb->table, jgiName); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { jgi = jgiGeneLoad(row+rowOffset); printf("JGI locus name: %s
\n",jgiName); printf("JGI gene symbol: %s
\n",jgi->jgiSymbol); printf("JGI gene description: %s
\n",jgi->jgiDescription); printf("JGI gene id: " "", jgi->jgiGeneId); printf("%s
\n", jgi->jgiGeneId); printf("GC content: %.0f %%
\n",jgi->jgiGc); @@ -1995,41 +1958,37 @@ } sqlFreeResult(&sr); hFreeConn(&conn); printTrackHtml(tdb); jgiGeneFree(&jgi); } void doPfamHit(struct trackDb *tdb, char *hitName) /* Handle the Pfam hits track. */ { struct lowelabPfamHits *pfamHit; char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlConnection *spConn = NULL; struct sqlResult *sr; - char *dupe, *words[16]; char **row; - int wordCount; int rowOffset; char *description; int start = cartInt(cart, "o"); spConn = sqlConnect(UNIPROT_DB_NAME); genericHeader(tdb,hitName); - dupe = cloneString(tdb->type); - wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database, seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s' and chrom = '%s' and chromStart = %d", tdb->table, hitName,seqName,start); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { pfamHit = lowelabPfamHitsLoad(row+rowOffset); sqlSafef(query, sizeof(query), "select description from proteome.pfamDesc where pfamAC='%s'", pfamHit->pfamAC); if (!sqlTableExists(spConn,"proteome.pfamDesc")) { sqlSafef(query, sizeof(query), "select extDbRef.extAcc1 from extDbRef,extDb " @@ -3279,39 +3238,30 @@ /* Close table */ printf("\n"); printf("\n"); printf("\n"); printf("\n"); bedFree(&primer); hFreeConn(&conn); printTrackHtml(tdb); } void doWiki(char *track, struct trackDb *tdb, char *itemName) { char strand[2]; - char wikiea[] = "wikiea"; - char wikibme[] = "microbewiki"; - char *wiki; - - if(sameWord(track, "wiki")) - wiki = wikiea; - else - wiki = wikibme; - printf(""); if(startsWith("Make", itemName)) { strand[0] = itemName[strlen(itemName)-1]; strand[1] = 0; printf("", database, seqName, winStart, winEnd, strand); } else { printf("", database, seqName, itemName); } printf(""); @@ -3406,57 +3356,53 @@ freeMem(rnaHyb); } sqlFreeResult(&sr); } hFreeConn(&conn); } void doarCOGs(struct trackDb *tdb, char *itemName) { char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; - char *dupe, *words[16]; char **row; - int wordCount; int rowOffset; struct arCOGs *infoload; int start = cartInt(cart, "o"); int end = cartInt(cart, "t"); sqlSafef(query, sizeof query, "select * from mgCommonDb.arcogdesc where name = '%s'", itemName); //sqlSafef(query, sizeof query, "select * from %s where query = '%s'", blastpHitsTable, queryName); struct sqlResult *srarcogdesc = sqlGetResult(conn, query); struct arcogdesc *description = NULL; while ((row = sqlNextRow(srarcogdesc)) != NULL) { struct arcogdesc *element; element = arcogdescLoad(row); slAddTail(&description, element); } sqlFreeResult(&srarcogdesc); genericHeader(tdb,itemName); - dupe = cloneString(tdb->type); - wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database,seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s' and chrom = '%s' and chromStart = %d and chromEnd = '%d';", tdb->table, itemName,seqName,start, end); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { infoload = arCOGsLoad(row+rowOffset); printf("Name: %s
\n", infoload->name); printf("Description: %s
\n", description->description); printf("Code: %s
\n", description->code); printf("Gene: %s
\n", infoload->gene); printf("Position: " "", @@ -3582,42 +3528,38 @@ printf("Genomic size: %d nt
\n", (infoload->chromEnd - infoload->chromStart)); if (infoload->next != NULL) printf("
\n"); } sqlFreeResult(&sr); hFreeConn(&conn); cddInfoFree(&infoload); printTrackHtml(tdb); } void domegablastInfo(struct trackDb *tdb, char *itemName) { char query[512]; struct sqlConnection *conn = hAllocConn(database); struct sqlResult *sr; - char *dupe, *words[16]; char **row; - int wordCount; int rowOffset; struct megablastInfo *infoload; int start = cartInt(cart, "o"); int end = cartInt(cart, "t"); genericHeader(tdb,itemName); - dupe = cloneString(tdb->type); - wordCount = chopLine(dupe, words); rowOffset = hOffsetPastBin(database,seqName, tdb->table); sqlSafef(query, sizeof query, "select * from %s where name = '%s' and chrom = '%s' and chromStart = %d and chromEnd = '%d';", tdb->table, itemName,seqName,start, end); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { infoload = megablastInfoLoad(row+rowOffset); printf("Name: %s
\n", infoload->name); printf("Accession: %s
", infoload->name, infoload->name); printf(" Link to NCBI Site
\n"); printf("Description: %s
\n", infoload->fullname); printf("E-value: %0.0e", infoload->evalue); #ifdef LISTUI