4898794edd81be5285ea6e544acbedeaeb31bf78
max
Tue Nov 23 08:10:57 2021 -0800
Fixing pointers to README file for license in all source code files. refs #27614
diff --git src/hg/gisaid/gisaidTable/getSeq.c src/hg/gisaid/gisaidTable/getSeq.c
index 0df055f..2c0f223 100644
--- src/hg/gisaid/gisaidTable/getSeq.c
+++ src/hg/gisaid/gisaidTable/getSeq.c
@@ -1,100 +1,100 @@
/* getSeq - pages to get protein and nucleic acid sequence. */
/* Copyright (C) 2013 The Regents of the University of California
- * See README in this or parent directory for licensing information. */
+ * See kent/LICENSE or http://genome.ucsc.edu/license/ for licensing information. */
#include "common.h"
#include "jksql.h"
#include "cart.h"
#include "dnautil.h"
#include "hdb.h"
#include "cheapcgi.h"
#include "hPrint.h"
#include "gisaidTable.h"
static void getSeqFromBlob(struct sqlConnection *conn,
struct subjInfo *siList, char *tableName, char *xrefField)
/* Get sequence from blob field in table and print it as fasta. */
{
struct sqlResult *sr;
char **row;
char query[256];
struct subjInfo *si;
int seqCnt = 0;
hPrintf("");
for (si = siList; si != NULL; si = si->next)
{
char *subjId = si->fields[1];
/* currently just 3 Thailand or 4 US */
sqlSafef(query, sizeof(query),
"select id, seq from %s s, gisaidXref g where g.subjId='%s' and g.%s=s.id",
tableName, subjId, xrefField);
sr = sqlGetResult(conn, query);
while ((row = sqlNextRow(sr)) != NULL)
{
char *id = row[0];
char *seq = row[1];
hPrintf(">%s", id);
hPrintf(":%s\n", subjId);
writeSeqWithBreaks(stdout, seq, strlen(seq), 60);
hPrintf("\n");
seqCnt++;
}
sqlFreeResult(&sr);
}
if (seqCnt == 0) hPrintf("No sequence data available for subject(s) selected.");
hPrintf("
");
}
static void getProtein( struct sqlConnection *conn, struct subjInfo *siList)
/* Print out proteins. */
{
getSeqFromBlob(conn, siList, "aaSeq", "aaSeqId");
}
static void getGenomic( struct sqlConnection *conn, struct subjInfo *siList)
/* Print out dna. */
{
getSeqFromBlob(conn, siList, "dnaSeq", "dnaSeqId");
}
void doGetSeq(struct sqlConnection *conn,
struct subjInfo *siList, char *how)
/* Put up the get sequence page. */
{
if (sameString(how, "protein"))
getProtein(conn, siList);
else if (sameString(how, "genomic"))
getGenomic(conn, siList);
else
errAbort("Unrecognized %s value %s", getSeqHowVarName, how);
}
static void howRadioButton(char *how)
/* Put up a getSeqHow radio button. */
{
char *howName = getSeqHowVarName;
char *oldVal = cartUsualString(cart, howName, "protein");
cgiMakeRadioButton(howName, how, sameString(how, oldVal));
}
void doGetSeqPage(struct sqlConnection *conn, struct column *colList)
/* Put up the get sequence page asking how to get sequence. */
{
hPrintf("