Same species lift over construction: Difference between revisions

From genomewiki
Jump to navigationJump to search
(add description of this step 1 procedure)
(add category tags)
 
(3 intermediate revisions by the same user not shown)
Line 1: Line 1:
==2018 UPDATE NOTE==
This page is an interesting historical discussion and well worth the read.
<b><span style="color:#bb0000">
HOWEVER, please note</span></b>,
the UCSC tool chain command: [[DoSameSpeciesLiftOver.pl]] can now perform this entire sequence of events
in your environment with your selected genome sequences.
=Same Species Lift Over=
=Same Species Lift Over=
Same species lift over procedure.  There is a different procedure for different species lift over construction: [[Whole_genome_alignment_howto]]  
Same species lift over procedure.  There is a different procedure for different species lift over construction: [[Whole_genome_alignment_howto]]  
Line 12: Line 21:
uses [[File:BlatJob.csh.txt]]
uses [[File:BlatJob.csh.txt]]


The 'sameSpeciesBlatSetup.sh' script partitions each genome sequence into single parts, expecting small genomes
The ''sameSpeciesBlatSetup.sh'' script partitions each genome sequence into single parts, expecting small genomes
with chromosome sizes less than 10,000,000 bases each so that no chromosome is broken into separate
with chromosome sizes less than 10,000,000 bases each so that no chromosome is broken into separate
parts.  The two resulting parts lists are matched together so that each part from each genome is matched
parts.  The two resulting parts lists are matched together so that each part from each genome is matched
with every other part of the other genome.  This generates a job-list shell script for each of those parings.
with every other part of the other genome.  This generates a job-list shell script for each of those parings.
The BlatJob.csh script runs the blat for each pair to generate PSL file results which are used in the chain/net step.
The ''BlatJob.csh'' script runs the blat for each pair to generate PSL file results which are used in the chain/net step.


=Step 2 - chain/net=
=Step 2 - chain/net=
The second step is encapsulated in the script: [[File:SameSpeciesChainNet.sh.txt]]
The second step is encapsulated in the script: [[File:SameSpeciesChainNet.sh.txt]]
This script uses the PSL files from the blat runs, chaining and netting to obtain a single
coverage lift over file.
Similar setup as for Step 1, working in the run.chain directory, constructs the result file.
Similar setup as for Step 1, working in the run.chain directory, constructs the result file.
Read the script for information on how to use it and where the result is.
Read the script for information on how to use it and where the result is.
[[Category:Cluster FAQ]]
[[Category:Technical FAQ]]

Latest revision as of 15:06, 26 April 2018

2018 UPDATE NOTE

This page is an interesting historical discussion and well worth the read.

HOWEVER, please note, the UCSC tool chain command: DoSameSpeciesLiftOver.pl can now perform this entire sequence of events in your environment with your selected genome sequences.

Same Species Lift Over

Same species lift over procedure. There is a different procedure for different species lift over construction: Whole_genome_alignment_howto See also: Minimal_Steps_For_LiftOver

Prerequisites

These procedures are unix command line shell scripted operations. Users should be familiar with unix command line shell programming and have the kent source tree available to build the necessary kent programs.

Step 1 - blat

Decide on a workDirectory where all this is to be performed. Note the script: File:SameSpeciesBlatSetup.sh.txt and the comments in it to setup the workDirectory and your two genome sequences. That script also uses File:BlatJob.csh.txt

The sameSpeciesBlatSetup.sh script partitions each genome sequence into single parts, expecting small genomes with chromosome sizes less than 10,000,000 bases each so that no chromosome is broken into separate parts. The two resulting parts lists are matched together so that each part from each genome is matched with every other part of the other genome. This generates a job-list shell script for each of those parings. The BlatJob.csh script runs the blat for each pair to generate PSL file results which are used in the chain/net step.

Step 2 - chain/net

The second step is encapsulated in the script: File:SameSpeciesChainNet.sh.txt This script uses the PSL files from the blat runs, chaining and netting to obtain a single coverage lift over file.

Similar setup as for Step 1, working in the run.chain directory, constructs the result file. Read the script for information on how to use it and where the result is.