SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Cuffmerge: merging reference and transcript GTFs? mra Bioinformatics 0 06-27-2013 11:08 AM
Create new Reference by merging SNV list with reference genome rdoan Bioinformatics 0 10-12-2012 07:17 AM
Cufflinks merging more than one transcript on bacterial genomes Noa Bioinformatics 2 01-24-2012 12:33 PM
Cufflinks merging more than one transcript on bacterial genomes Noa RNA Sequencing 0 01-24-2012 05:41 AM
Merging 2 genomes oiiio Bioinformatics 0 01-25-2011 03:16 PM

Reply
 
Thread Tools
Old 01-20-2015, 05:58 AM   #1
g.galazzo
Junior Member
 
Location: Groningen

Join Date: Dec 2014
Posts: 2
Default Merging different reference genomes

Hi to everyone,
i'm trying to replicate a pipeline from an article in which they says that they had merged different information from different reference genomes to build an unique and complete reference genome.
The article is by Billy Li and here is the link

http://www.nature.com/nmeth/journal/...nmeth.1982.pdf

and what they said is : " we mapped RNA-seq
reads not only to the reference genome or to the transcriptome
but to a combination of the hg19 reference genome plus exonic
sequences surrounding all currently known splicing junctions
from gene models available in annotation from Gencode, RefSeq,
Ensembl and UCSC Genes"

So here is my problem, i don't understand how they performed this task.
I tried in different ways to merge the reference genomes,
i tried merging the genome.fa file from Wholegenome folder from my reference files,
i tried to merge the chromosome files together and after fuse them in a unique file,
i also tried to merge the genome.fa from the BWAIndex folder that is already present in to the reference folders, but nothing works.
I index the merged files but the program doesn't recognize the index.
does anyone have and idea about how to perform this merging?

thanks all
g.galazzo is offline   Reply With Quote
Old 01-20-2015, 06:27 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Honestly, it would make more sense to just use a split-read mapper like STAR, tophat, or segemehl. You can then just map to the genome and be done with that step instead of needing to play around with adding splice junction sequences in.
dpryan is offline   Reply With Quote
Reply

Tags
bwa, index, reference genome

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:30 PM.


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