Question: How to create chain file to change genome coordinates from one genome version to another
3
gravatar for iraun
6.1 years ago by
iraun3.8k
Norway
iraun3.8k wrote:

Hi there,

I'm looking for a 'protocol' or I would like to know the best way to create a chain file in order to get the new coordinates of a new assembly of a given genome. I have an old genome, with a gtf file, and a new genome (only the fasta). I would like to create a new GTF, and I'm thinking about Crossmap, but first, I need this chain file... I have read something about blat if the specie is the same (my case)..: http://genomewiki.ucsc.edu/index.php/Minimal_Steps_For_LiftOver , but my genome is larger than 4Mb!

Anyone knows a pipeline o protocol for doing this?

Thanks in advance,

liftover remap chain format • 5.2k views
ADD COMMENTlink modified 5.9 years ago by Biostar ♦♦ 20 • written 6.1 years ago by iraun3.8k
1

Hi,

I run into the same problem and Crossmap looks promising.. However, it looks like the chain file is critical.. 

Have you come up with an alternative strategy to resolve this issue? Any help/guide is more than appreciated.. Thanks..

ADD REPLYlink written 5.9 years ago by Jada10

I just found this: https://github.com/wurmlab/flo haven't tried yet, but it seems to do the job...

ADD REPLYlink written 4.1 years ago by fransua390

what about http://genomewiki.ucsc.edu/index.php/Same_species_lift_over_construction ?

ADD REPLYlink written 6.1 years ago by Pierre Lindenbaum131k

Thanks Pierre. Yes, I have already read that link but... "expecting small genomes with chromosome sizes less than 10,000,000 bases each ". What's up if my genome is larger?

ADD REPLYlink written 6.1 years ago by iraun3.8k

Maybe I can change the parameter export targetChunkSize="10000000" (inside the script), up to my genome size but... maybe the rest of variables in the script are set to that size and not a bigger size..

ADD REPLYlink modified 6.1 years ago • written 6.1 years ago by iraun3.8k

the 'targetChunkSize' is not the genomeSize. The UCSC use that script to generate the chain files for the  human genome.

ADD REPLYlink written 6.1 years ago by Pierre Lindenbaum131k

Look this: https://groups.google.com/a/soe.ucsc.edu/forum/#!searchin/genome/chain/genome/8GlealJxiRw/l8b7u6PeKgoJ , looking the last comment... he says that those scripts are only for small genomes...

ADD REPLYlink written 6.1 years ago by iraun3.8k
Please log in to add an answer.

Help
Access

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