Question: Mapping using UCSC and counting using RefSeq
gravatar for Fawzi Yassine
21 months ago by
Fawzi Yassine10 wrote:

Hi, I am analyzing public RNA-seq data to find genes differentially expressed between diseased patients and controls. The paper for which the paper belong says: ” The reads were mapped to human genome (hg19) by TopHat. The number of fragments in each known gene from RefSeq database (downloaded from UCSC Genome Browser) was enumerated….”

Q1: Are the 2 files compatible? (one file used for Mapping using UCSC hg19 Genome but another file used for annotation was downloaded from TableBrowser using track Refseq ) Q2: which track shall I use? ‘NCBI Refseq’ or ‘Other Refseq'. Thank you,

ADD COMMENTlink modified 21 months ago by David Langenberger9.5k • written 21 months ago by Fawzi Yassine10
gravatar for David Langenberger
21 months ago by
David Langenberger9.5k wrote:

The two files are compatible, if for the gene annotation file the 'hg19' assembly (same as genome assembly) was selected. The genome file holds the sequence of the human genome and the gene annotation file holds the predicted genes from the RefSeq database for this genome. First file is a fasta file and the second file probably a bed file with genomic locations of genes. Use the NCBI RefSeq file and not the 'Others RefSeq'.

ADD COMMENTlink written 21 months ago by David Langenberger9.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1201 users visited in the last hour