static void cytoBandDrawAt(struct track *tg, void *item,
struct hvGfx *hvg, int xOff, int y, double scale,
MgFont *font, Color color, enum trackVisibility vis)
/* Draw cytoBand items. */
{
struct cytoBand *band = item;
int heightPer = tg->heightPer;
int x1,x2,w;
x1 = round((double)((int)tg->itemStart(tg,band) - winStart)*scale) + xOff;
x2 = round((double)((int)tg->itemEnd(tg,band) - winStart)*scale) + xOff;
/* Clip here so that text will tend to be more visible... */
if (x1 < xOff)
x1 = xOff;
if (x2 > insideX + insideWidth)
x2 = insideX + insideWidth;
w = x2-x1;
if (w < 1)
w = 1;
hCytoBandDrawAt(band, hvg, x1, y, w, heightPer, hCytoBandDbIsDmel(database), font,
mgFontPixelHeight(font), tg->ixColor, tg->ixAltColor,
shadesOfGray, maxShade);
if(tg->mapsSelf)
tg->mapItem(tg, hvg, band, band->name, band->name, band->chromStart, band->chromEnd,
x1, y, w, heightPer);
else
mapBoxHc(hvg, band->chromStart, band->chromEnd, x1,y,w,heightPer, tg->track,
band->name, band->name);
}
void genoLayDrawBandedChroms(struct genoLay *gl, struct hvGfx *hvg, char *db,
struct sqlConnection *conn, Color *shadesOfGray, int maxShade,
int defaultColor)
/* Draw chromosomes with centromere and band glyphs.
* Get the band data from the database. If the data isn't
* there then draw simple chroms in default color instead */
{
char *bandTable = "cytoBandIdeo";
int yOffset = gl->chromOffsetY;
genoLayDrawSimpleChroms(gl, hvg, defaultColor);
if (sqlTableExists(conn, bandTable) && !gl->allOneLine)
{
int centromereColor = hCytoBandCentromereColor(hvg);
double pixelsPerBase = 1.0/gl->basesPerPixel;
int height = gl->chromHeight;
int innerHeight = gl->chromHeight-2;
struct genoLayChrom *chrom;
boolean isDmel = hCytoBandDbIsDmel(db);
boolean bColor = hvGfxFindColorIx(hvg, 200, 150, 150);
int fontPixelHeight = mgFontPixelHeight(gl->font);
for (chrom = gl->chromList; chrom != NULL; chrom = chrom->next)
{
boolean gotAny = FALSE;
struct sqlResult *sr;
char **row;
char query[256];
int cenX1=BIGNUM, cenX2=0;
int y = chrom->y + yOffset;
/* Fetch bands from database and draw them. */
safef(query, sizeof(query), "select * from %s where chrom='%s'",
bandTable, chrom->fullName);
sr = sqlGetResult(conn, query);
while ((row = sqlNextRow(sr)) != NULL)
{
struct cytoBand band;
int x1, x2;
cytoBandStaticLoad(row, &band);
x1 = pixelsPerBase*band.chromStart;
x2 = pixelsPerBase*band.chromEnd;
if (sameString(band.gieStain, "acen"))
{
/* Centromere is represented as two adjacent bands.
* We'll just record the extents of it here, and draw it
* in one piece later. */
if (x1 < cenX1)
cenX1 = x1;
if (x2 > cenX2)
cenX2 = x2;
}
else
{
/* Draw band */
hCytoBandDrawAt(&band, hvg, x1+chrom->x, y+1, x2-x1, innerHeight,
isDmel, gl->font, fontPixelHeight, MG_BLACK, bColor,
shadesOfGray, maxShade);
gotAny = TRUE;
}
}
sqlFreeResult(&sr);
/* Draw box around chromosome */
hvGfxBox(hvg, chrom->x, y, chrom->width, 1, MG_BLACK);
hvGfxBox(hvg, chrom->x, y+height-1, chrom->width, 1, MG_BLACK);
hvGfxBox(hvg, chrom->x, y, 1, height, MG_BLACK);
hvGfxBox(hvg, chrom->x+chrom->width-1, y, 1, height, MG_BLACK);
/* Draw centromere if we found one. */
if (cenX2 > cenX1)
{
hCytoBandDrawCentromere(hvg, cenX1+chrom->x, y, cenX2-cenX1, height,
MG_WHITE, centromereColor);
}
}
}
}