This directory contain human/rat alignments made
using the July 2003 human assembly (also known as build 34)
vs. the June 2003 rat assembly (also known as rn3).
The subdirectory axtTight contains relatively stringent
human/rat alignments filtered so that only the best
alignment for any given region of the human genome
is used.
The alignments are in 'axt' format. Each alignment
contains three lines and is separated from the next
alignment by a space:
Line 1 - summarizes the alignment.
Line 2 - contains the human sequence with inserts.
Line 3 - contains the rat sequence with inserts.
The summary line contains 9 blank separated fields with the
following meanings:
1 - Alignment number. The first alignment in a file
is numbered 0, the next 1, and so forth.
2 - Human chromosome.
3 - Start in human chromosome. The first base is
numbered 1.
4 - End in human chromosome. The end base is included.
5 - Rat chromosome.
6 - Start in rat.
7 - End in rat.
8 - Rat strand. If this is '-', the rat start/end fields
are relative to the reverse-complemented rat chromosome.
9 - Blastz score. The scoring matrix blastz uses is:
A C G T
A 91 -114 -31 -123
C -114 100 -125 -31
G -31 -125 100 -114
T -123 -31 -114 91
with a gap open penalty of 400 and a gap extension
penalty of 30. The minimum score for an alignment
to be kept was 3000 for the first pass, and then
2200 for the second pass, which just restricts
the search space to the regions between two alignments
found in the first pass.
The alignments were done with blastz, which is available
from Webb Miller's group at Pennsylvania State University (PSU).
Each chromosome was divided into 10010000 base chunks with 10000
bases of overlap. The .lav format blastz output, which does not
include the sequence, was converted to .axt with PSU's lavToAxt.
The axtTight alignments were processed with subsetAxt from
Jim Kent using the matrix:
A C G T
A 100 -200 -100 -200
C -200 100 -200 -100
G -100 -200 100 -200
T -200 -100 -200 100
with a gap open penalty of 2000 and a gap extension
penalty of 50. The minimum score was 3400. The axtTight
subset covers 6% of the human genome while axtBest covers
40%.
Name Last modified Size Description
Parent Directory -
chrY.axt.gz 2003-08-26 10:22 268K
chrX_random.axt.gz 2003-08-26 10:22 88K
chrX.axt.gz 2003-08-26 10:22 6.1M
chrUn_random.axt.gz 2003-08-26 10:22 83K
chrM.axt.gz 2003-08-26 10:22 127
chr9_random.axt.gz 2003-08-26 10:22 71K
chr9.axt.gz 2003-08-26 10:22 4.4M
chr8_random.axt.gz 2003-08-26 10:22 22K
chr8.axt.gz 2003-08-26 10:22 4.6M
chr7_random.axt.gz 2003-08-26 10:22 8.0K
chr7.axt.gz 2003-08-26 10:22 5.5M
chr6_random.axt.gz 2003-08-26 10:22 32K
chr6.axt.gz 2003-08-26 10:22 5.6M
chr5_random.axt.gz 2003-08-26 10:21 3.3K
chr5.axt.gz 2003-08-26 10:21 6.7M
chr4_random.axt.gz 2003-08-26 10:21 12K
chr4.axt.gz 2003-08-26 10:21 5.7M
chr3_random.axt.gz 2003-08-26 10:21 32K
chr3.axt.gz 2003-08-26 10:21 7.4M
chr22.axt.gz 2003-08-26 10:21 1.1M
chr2_random.axt.gz 2003-08-26 10:21 12K
chr21.axt.gz 2003-08-26 10:21 912K
chr20.axt.gz 2003-08-26 10:21 2.2M
chr2.axt.gz 2003-08-26 10:21 9.2M
chr19.axt.gz 2003-08-26 10:21 2.0M
chr1_random.axt.gz 2003-08-26 10:21 158K
chr18_random.axt.gz 2003-08-26 10:21 37
chr18.axt.gz 2003-08-26 10:21 2.4M
chr17_random.axt.gz 2003-08-26 10:21 42K
chr17.axt.gz 2003-08-26 10:21 3.8M
chr16.axt.gz 2003-08-26 10:21 3.4M
chr15_random.axt.gz 2003-08-26 10:20 32K
chr15.axt.gz 2003-08-26 10:20 3.6M
chr14.axt.gz 2003-08-26 10:20 3.5M
chr13_random.axt.gz 2003-08-26 10:20 8.2K
chr13.axt.gz 2003-08-26 10:20 3.0M
chr12.axt.gz 2003-08-26 10:20 4.5M
chr11.axt.gz 2003-08-26 10:20 5.5M
chr10_random.axt.gz 2003-08-26 10:20 26K
chr10.axt.gz 2003-08-26 10:20 4.7M
chr1.axt.gz 2003-08-26 10:20 8.9M