bc21bd3d27fe3d29971231955b3fc544fa1c3d1e angie Wed Oct 16 11:51:39 2013 -0700 Two new tracks for Locus Reference Genomic (LRG) (#11863) with customhandlers: LRG Regions and LRG Transcripts. LRGs are frozen reference sequences for a particular gene plus some upstream and downstream sequence. They are intended to provide a stable coordinate system for gene annotations that won't change with every new genome assembly, but can be mapped to each genome assembly. Since there is a lot of metadata associated with each region, I made LRG Regions a bigBed 12 + with fields describing mismatches and indels, so that PSL can be derived from the bigBed and the original LRG sequence can be reconstructed using genome assembly sequence and the mismatch/indel info. hgTracks shows differences and LRG insertions into the reference assembly using the cds.c baseColor code. (LRG deletions from the reference appear as gaps, which we get for free with bed12 info). For LRG Transcripts, I found the genePred codon-coloring code inadequate for showing an insertion into hg19 (or even mismatches), so instead of genePred I ended up using PSL + sequence, more like the mRNA track representation and display. diff --git src/hg/hgTrackUi/hgTrackUi.c src/hg/hgTrackUi/hgTrackUi.c index cec6d00..8157a80 100644 --- src/hg/hgTrackUi/hgTrackUi.c +++ src/hg/hgTrackUi/hgTrackUi.c @@ -2734,30 +2734,34 @@ else if (sameString(track, "switchDbTss")) switchDbScoreUi(tdb); else if (sameString(track, "dgv") || (startsWith("dgvV", track) && isdigit(track[4]))) dgvUi(tdb); else if (sameString(track, "all_mrna") || sameString(track, "mrna") || sameString(track, "all_est") || sameString(track, "est") || sameString(track, "tightMrna") || sameString(track, "tightEst") || sameString(track, "intronEst") || sameString(track, "xenoMrna") || sameString(track, "xenoEst")) mrnaCfgUi(cart, tdb, tdb->track, NULL, boxed); +else if (sameString(track, "lrg")) + lrgCfgUi(cart, tdb, tdb->track, NULL, boxed); +else if (sameString(track, "lrgTranscriptAli")) + lrgTranscriptAliCfgUi(cart, tdb, tdb->track, NULL, boxed); else if (tdb->type != NULL) { // NOTE for developers: please avoid special cases and use cfgTypeFromTdb//cfgByCfgType() // When you do, then multi-view cfg and subtrack cfg will work. eCfgType cType = cfgTypeFromTdb(tdb,FALSE); if (cType != cfgNone) { cfgByCfgType(cType,database, cart, tdb,tdb->track, NULL, boxed); } // NOTE: these cases that fall through the cracks should probably get folded into cfgByCfgType() else if (startsWithWord("expRatio", tdb->type)) expRatioUi(tdb); else if (startsWith("chromGraph", tdb->type)) chromGraphUi(tdb); else if (startsWith("sample", tdb->type)) genericWiggleUi(tdb,7);