Schema for TransMap ESTs - TransMap EST Mappings Version 5
  Database: hg38    Primary Table: transMapEstV5 Data last updated: 2019-06-11
Big Bed File: /gbdb/hg38/transMap/V5/hg38.est.transMapV5.bigPsl
Item Count: 4,722,982
Format description: bigPsl derived pairwise alignment with additional information
fieldexampledescription
chromchr1Reference sequence chromosome or scaffold
chromStart165937016Start position in chromosome
chromEnd166062379End position in chromosome
namerheMac8:CB230316.1-1.1alignment Id
score928Score (0-1000), faction identity * 1000
strand++ or - indicates whether the query aligns to the + or - strand on the reference
thickStart166062379Start of where display should be thick (start codon)
thickEnd166062379End of where display should be thick (stop codon)
reserved0RGB value (use R,G,B string in input file)
blockCount13Number of blocks
blockSizes374,99,11,64,28,18,5,10,5,17,5,12,10,Comma separated list of block sizes
chromStarts0,5668,5768,125177,125242,125270,125288,125294,125309,125314,125331,125341,125353,Start positions relative to chromStart
oChromStart3Start position in other chromosome
oChromEnd681End position in other chromosome
oStrand++ or -, - means that psl was reversed into BED-compatible coordinates
oChromSize789Size of other chromosome.
oChromStarts3,377,476,487,551,580,599,606,623,629,647,658,671,Start positions relative to oChromStart or from oChromStart+oChromSize depending on strand
oSequencegggactagatgccagtagatctgccttccagttgtaacaaccaaaactgtctacagacattgcctaatgtccgcctgtgggtcaaaatcgtccctgttttgagaaccaaggatccagacccgtgcttctcaaaatacagtccgtggtccagtaagcgtctgcctcatctagaactttgttagaagtgaagactcttgaggcccactccaaatctactgactcagaaactgcattttaacaaggtcctcaggtgattcatatgcagttagaattgacaagcgctgacctaaaccaacttcttcagtttacaagtgaggagatggctttccagagggcagatcacctactcaaggtcatgtctcctgaatcctgcacatgaaagtttcaacacagtggaagagacagaccgatcaaaactgcaatacaaagcgaaatgtaaaacggagtaagcgaggcacccccaaagtgctacgggattgagaagaagacggaaaagcgaagatccaagagatggcctgcatcaaacccccatcaagccaaaatcaggcctcggggtctcaagtctcctgggcccggacacccaccctgtgccaatgaaacactggggtttgatacacaaaaaccgggctttggacagatctaattcctactggaaaaaaaaaagacaggaatcagggcaaggctggcaaaaaccttcgggaaatccaactggagttttaaaagcctggcttaattccttttccacaaaaaaatggcgccccaatttttcccctttctggacaaaaaaaaSequence on other chrom (or edit list, or empty)
oCDSCDS in NCBI format
chromSize248956422Size of target chromosome
match256Number of bases matched.
misMatch47 Number of bases that don't match
repMatch355 Number of bases that match but are part of repeats
nCount0 Number of 'N' bases
seqType10=empty, 1=nucleotide, 2=amino_acid
srcDbrheMac8source database
srcTransIdCB230316.1source transcript id
srcChromchr1source chromosome
srcChromStart201533114start position in source chromosome
srcChromEnd201667101end position in source chromosome
srcIdent986source score (fraction identity * 1000)
srcAligned833fraction of source transcript aligned (fraction aligned * 1000)
geneNamegene name
geneIdgene id
geneTypegene type
transcriptTypetranscript type
chainTypesyntype of chains used for mapping
commonNameRhesuscommon name
scientificNameMacaca mulattascientific name
orgAbbrevMacacaorganism abbreviation

Sample Rows
 
chromchromStartchromEndnamescorestrandthickStartthickEndreservedblockCountblockSizeschromStartsoChromStartoChromEndoStrandoChromSizeoChromStartsoSequenceoCDSchromSizematchmisMatchrepMatchnCountseqTypesrcDbsrcTransIdsrcChromsrcChromStartsrcChromEndsrcIdentsrcAlignedgeneNamegeneIdgeneTypetranscriptTypechainTypecommonNamescientificNameorgAbbrev
chr1165937016166062379rheMac8:CB230316.1-1.1928+166062379166062379013374,99,11,64,28,18,5,10,5,17,5,12,10,0,5668,5768,125177,125242,125270,125288,125294,125309,125314,125331,125341,125353,3681+7893,377,476,487,551,580,599,606,623,629,647,658,671,gggactagatgccagtagatctgccttccagttgtaacaaccaaaactgtctacagacattgcctaatgtccgcctgtgggtcaaaatcgtccctgttttgagaaccaaggatccagacccgtgcttc ...2489564222564735501rheMac8CB230316.1chr1201533114201667101986833synRhesusMacaca mulattaMacaca
chr1166060009166069884mm10:AV381729.2-1.1760-16606988416606988402057,31,14,21,24,24,30,22,32,30,45,39,32,11,65,56,18,20,27,52,0,57,93,126,147,171,201,232,254,288,320,366,405,455,468,533,589,607,627,9823,1678+6780,58,89,103,137,164,188,218,242,274,304,349,389,421,432,498,555,574,598,625,cgatgttgcttcccaaccaattagaagcgagagtgaaaacatccacaccaggtgtctcctagtaggacaaggcaagcccacccagcaagcagaaccacatggcaagacgctgagaaggagggatatga ...2489564224551563901mm10AV381729.2chr1167079415167088602997998synMouseMus musculusMus
chr1166060319166070238mm10:CA327390.1-1.1844-1660702381660702380148,45,39,32,11,65,56,18,20,27,304,83,9,7,0,10,56,95,145,158,223,279,297,317,9515,9819,9902,9912,0735+7350,8,53,93,125,136,202,259,278,302,329,634,719,728,gcagactgcccaactgnntgtgtggaacacatgccactgagacccgcagtggggtgaccaggggtggcaggagccaaactgagtttctgagccaaagcagaccctcttggtgtgccagcctttgcagc ...2489564225701123931mm10CA327390.1chr1167079061167088306998993synMouseMus musculusMus
chr1166070081166166146rn6:CX569547.1-1.1945-16616614616616614601315,17,5,40,41,7,532,86,5,5,1,3,58,0,36,53,59,99,142,150,95904,95991,95996,96001,96002,96007,0848+92274,115,133,140,181,222,229,761,847,853,859,861,864,gaggagacctcgcccatcgtcctgcgctacaagaccccttacttcaaagcctcggcccgtgtggtcatgccccccatcccccgccacgagacctgggtggtgggctggatccagtcttgcaaccagat ...248956422761440101rn6CX569547.1chr138511112185231314998874synRatRattus norvegicusRattus
chr1166070218166166182mm10:CB518481.1-1.1955-1661661821661661820312,532,197,0,13,95767,0741+7443,15,547,gtgtgcgccaccatcgaccagtgccccacgcgcatcgaggagacctcgcccatcgtcctgcgctacaagaccccttacttcaaagcctcggcccgtgtggtcatgccccccatcccccgccacgagac ...24895642270533031mm10CB518481.1chr1167001631167079080997991synMouseMus musculusMus
chr1166070238166070763xenTro9:EL832716.1-1.1809-16607076316607076301525,0,58583+700117,aaaatggtggctcaagatgaaagaatgcaatatctacttatatatatactccagtttggtcaagctgggagcttcctgatttgagggaaggaagagtaaaagcaatcagcgattctgatggtgtaagt ...248956422425100001xenTro9EL832716.1chr410896877110900485910001000rbestX. tropicalisXenopus tropicalisXenopu
chr1166070361166166234mm10:BB647672.1-1.1943-1661662341661662340614,7,6,11,348,249,0,22,37,43,54,95624,1660+70545,59,88,95,107,455,taagcatcacctgcaaggcgcggatccggcgcgagaacatcgtggtgtacgatgtgtgcgccaccatcgaccagtgccccacgcgcatcgaggagacctcgcccatcgtcctgcgctacaagacccct ...24895642259936001mm10BB647672.1chr1167001579167078938981900synMouseMus musculusMus
chr1166070690166166366rn6:FM038873.1-1.1975-16616636616616636601673,86,5,5,1,3,82,53,24,31,18,10,1,47,3,6,0,95295,95382,95387,95392,95393,95398,95481,95534,95558,95589,95607,95618,95619,95666,95670,0455+4550,73,159,165,171,173,176,259,312,336,367,386,396,398,446,449,gctgagctggaccgggagcgagcggcgcgtcgccatgccggcgtgacgggcgcccgcggctgcccgcgccgggcccccgcgctgccccacgccgcgccgcgccggcgccgggcggcaggatgggctgt ...24895642243711001rn6FM038873.1chr138510997685230708995993synRatRattus norvegicusRattus
chr1166165901166185318mm10:BY091203.1-1.1764+166185318166185318065,68,79,19,29,93,0,6,74,19262,19281,19324,1348+3481,6,79,205,226,255,ggaggggttggggagggggttcagaatggagagtggagaaaggggccaggaccccgctcccgcgccgcgtggcgtggccggtcgcttacatgcccaggtcgctgtaggtgttgaagaactccatctgg ...24895642222469001mm10BY091203.1chr11669745901670019161000997synMouseMus musculusMus
chr1166165985166166313xenTro9:CX796837.1-1.1802-16616631316616631304263,9,10,39,0,269,278,289,453775+86489,352,362,372,atgagagagagagagcggctgaatccgagcgcagggaaggattcaccatccgacaaggacatcagcccctggccacctgggacatgacaagccctggcagccccggcagcagggcacaggagcagcgg ...24895642225663021xenTro9CX796837.1chr4108969323109020500995997rbestX. tropicalisXenopus tropicalisXenopu

TransMap ESTs (transMapEstV5) Track Description
 

Description

This track contains GenBank spliced EST alignments produced by the TransMap cross-species alignment algorithm from other vertebrate species in the UCSC Genome Browser. For closer evolutionary distances, the alignments are created using syntenically filtered BLASTZ alignment chains, resulting in a prediction of the orthologous genes in human.

Display Conventions and Configuration

This track follows the display conventions for PSL alignment tracks.

This track may also be configured to display codon coloring, a feature that allows the user to quickly compare cDNAs against the genomic sequence. For more information about this option, click here. Several types of alignment gap may also be colored; for more information, click here.

Methods

  1. Source transcript alignments were obtained from vertebrate organisms in the UCSC Genome Browser Database. BLAT alignments of RefSeq Genes, GenBank mRNAs, and GenBank Spliced ESTs to the cognate genome, along with UCSC Genes, were used as available.
  2. For all vertebrate assemblies that had BLASTZ alignment chains and nets to the human (hg38) genome, a subset of the alignment chains were selected as follows:
    • For organisms whose branch distance was no more than 0.5 (as computed by phyloFit, see Conservation track description for details), syntenic filtering was used. Reciprocal best nets were used if available; otherwise, nets were selected with the netfilter -syn command. The chains corresponding to the selected nets were used for mapping.
    • For more distant species, where the determination of synteny is difficult, the full set of chains was used for mapping. This allows for more genes to map at the expense of some mapping to paralogous regions. The post-alignment filtering step removes some of the duplications.
  3. The pslMap program was used to do a base-level projection of the source transcript alignments via the selected chains to the human genome, resulting in pairwise alignments of the source transcripts to the genome.
  4. The resulting alignments were filtered with pslCDnaFilter with a global near-best criteria of 0.5% in finished genomes (human and mouse) and 1.0% in other genomes. Alignments where less than 20% of the transcript mapped were discarded.

To ensure unique identifiers for each alignment, cDNA and gene accessions were made unique by appending a suffix for each location in the source genome and again for each mapped location in the destination genome. The format is:

   accession.version-srcUniq.destUniq
Where srcUniq is a number added to make each source alignment unique, and destUniq is added to give the subsequent TransMap alignments unique identifiers.

For example, in the cow genome, there are two alignments of mRNA BC149621.1. These are assigned the identifiers BC149621.1-1 and BC149621.1-2. When these are mapped to the human genome, BC149621.1-1 maps to a single location and is given the identifier BC149621.1-1.1. However, BC149621.1-2 maps to two locations, resulting in BC149621.1-2.1 and BC149621.1-2.2. Note that multiple TransMap mappings are usually the result of tandem duplications, where both chains are identified as syntenic.

Data Access

The raw data for these tracks can be accessed interactively through the Table Browser or the Data Integrator. For automated analysis, the annotations are stored in bigPsl files (containing a number of extra columns) and can be downloaded from our download server, or queried using our API. For more information on accessing track data see our Track Data Access FAQ. The files are associated with these tracks in the following way:

  • TransMap Ensembl - hg38.ensembl.transMapV4.bigPsl
  • TransMap RefGene - hg38.refseq.transMapV4.bigPsl
  • TransMap RNA - hg38.rna.transMapV4.bigPsl
  • TransMap ESTs - hg38.est.transMapV4.bigPsl
Individual regions or the whole genome annotation can be obtained using our tool bigBedToBed which can be compiled from the source code or downloaded as a precompiled binary for your system. Instructions for downloading source code and binaries can be found here. The tool can also be used to obtain only features within a given range, for example:

bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/transMap/V4/hg38.refseq.transMapV4.bigPsl -chrom=chr6 -start=0 -end=1000000 stdout

Credits

This track was produced by Mark Diekhans at UCSC from cDNA and EST sequence data submitted to the international public sequence databases by scientists worldwide and annotations produced by the RefSeq, Ensembl, and GENCODE annotations projects.

References

Siepel A, Diekhans M, Brejová B, Langton L, Stevens M, Comstock CL, Davis C, Ewing B, Oommen S, Lau C et al. Targeted discovery of novel human exons by comparative genomics. Genome Res. 2007 Dec;17(12):1763-73. PMID: 17989246; PMC: PMC2099585

Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008 Mar 1;24(5):637-44. PMID: 18218656

Zhu J, Sanborn JZ, Diekhans M, Lowe CB, Pringle TH, Haussler D. Comparative genomics search for losses of long-established genes on the human lineage. PLoS Comput Biol. 2007 Dec;3(12):e247. PMID: 18085818; PMC: PMC2134963