SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
hg19/hg19 liftover Oscar.K Bioinformatics 1 03-14-2014 08:25 AM
How to get sequence from hg19.fa lewewoo Bioinformatics 3 09-03-2012 10:48 PM
dbSNP129 for hg19 Jane M Bioinformatics 5 07-02-2012 02:41 AM
How to get hg19.fa? ninad Bioinformatics 26 08-24-2011 07:40 AM
hg18 or hg19 Layla Bioinformatics 1 11-09-2009 11:38 PM

Reply
 
Thread Tools
Old 10-01-2014, 08:03 AM   #1
spaceace507
Junior Member
 
Location: Ann Arbor

Join Date: Jun 2013
Posts: 4
Default Cufflinks hg19 vs hg19_Broad_variant

I am trying to use Cufflinks (Cuffquant and Cuffnorm) to generate gene/transcript expression levels from the CCLE RNAseq data on cghub. The RNAseq data is provided as bam files aligned to genome assembly hg19_Broad_variant. The transcriptome assembly I would like to use is based on hg19/GRCh37. I have not found much information on hg19_Broad_variant to understand how different it is from hg19/GRCh37.

Do I need to convert the CCLE bam files to hg19/GRCh37 before putting them through a Cufflinks workflow? Or is hg19_Broad_variant close enough that it will work without conversion?

If I have to convert to hg19/GRCh37, is it possible to liftOver bam files? Or do I have to do bam->fastq using hg19_Broad_variant, then start over with alignment to hg19/GRCh37?

Thanks in advance for any advice.
spaceace507 is offline   Reply With Quote
Old 10-02-2014, 12:19 PM   #2
ske
Junior Member
 
Location: Santa Cruz

Join Date: Oct 2014
Posts: 1
Default HG19_Broad_variant

This is the explanation I received from the Broad:

The file is essentially what UCSC uses for the HG19 primary assembly plus EBV.

All the files we submit have a sequence dictionary that clearly gives
the accession numbers for all of the sequences. In addition it
contains a URL to the broad's public FTP site for the fasta file.
ske is offline   Reply With Quote
Old 01-26-2016, 09:29 AM   #3
icortes
Junior Member
 
Location: MA

Join Date: Jan 2016
Posts: 1
Default

So I guess that the gene coordinates do not change between UCSC hg19 and hg19_broad. Am I right?
Thanks a lot!
icortes is offline   Reply With Quote
Reply

Tags
broad_variant cufflinks

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:57 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO