feeb7ee7651d3a85ebc34c4fb32243a841eacdf2
braney
  Mon Apr 16 16:23:04 2012 -0700
on-going work on annoGrator (#6152).  Started using Sequence Ontology numbers from EBI, and added a generic variant structure as input to gpFx library
diff --git src/hg/lib/variant.c src/hg/lib/variant.c
index 8719d12..5e7ebb9 100644
--- src/hg/lib/variant.c
+++ src/hg/lib/variant.c
@@ -1,32 +1,36 @@
+/* variant.c -- routines to convert other variant formats to a generic
+ *              variant structure */
+
 #include "common.h"
 #include "variant.h"
 
 struct variant *variantFromPgSnp(struct pgSnp *pgSnp)
 /* convert pgSnp record to variant record */
 {
 struct variant *variant;
 
 // this is probably the wrong way to do this.  Alleles in
 // variant should be their size in query bases
 int alleleLength = pgSnp->chromEnd - pgSnp->chromStart;
 
 // We have a new variant!
 AllocVar(variant);
 variant->chrom = cloneString(pgSnp->chrom);
 variant->chromStart = pgSnp->chromStart;
 variant->chromEnd = pgSnp->chromEnd;
+variant->numAlleles = pgSnp->alleleCount;
 
 // get the alleles.
 char *nextAlleleString = pgSnp->name;
 int alleleNumber = 0;
 for( ; alleleNumber < pgSnp->alleleCount; alleleNumber++)
     {
     if (nextAlleleString == NULL)
 	errAbort("number of alleles in pgSnp doesn't match number in name");
     
     char *thisAlleleString = nextAlleleString;
 
     // advance pointer to next variant string
     // probably there's some kent routine to do this behind the curtain
     nextAlleleString = strchr(thisAlleleString, '/');
     if (nextAlleleString)	 // null out '/' and move to next char