Question: Does anyone know the genome alignment format like this?
0
gravatar for wang.yiguan
2.6 years ago by
wang.yiguan30
wang.yiguan30 wrote:

I have the whole genome alignment results from two species, but I am confused about the results format. Does anyone know the format, and can you make an explanation? Thank you very much!

##matrix=axtChain 16 91,-114,-31,-123,-114,100,-125,-31,-31,-125,100,-114,-123,-31,-114,91
##gapPenalties=axtChain O=400 E=30
# lastz.v1.04.00 --masking=50 --hspthresh=2200 --ydrop=3400 --gappedthresh=4000 --inner=2000 --format=axt --output=xxxxxx.axt
#
# hsp_threshold      = 2200
# gapped_threshold   = 4000
# x_drop             = 910
# y_drop             = 3400
# gap_open_penalty   = 400
# gap_extend_penalty = 30
#        A    C    G    T
#   A   91 -114  -31 -123
#   C -114  100 -125  -31
#   G  -31 -125  100 -114
#   T -123  -31 -114   91
chain 589271 MTTC0100016x.1 36216 + 2 13582 scf718000060xxxx 211718 + 199665 211441 1250
16      4       0
59      0       8
5       0       12
9       0       11
54      19      0
9       140     118
34      1       0
29      2       0
8       20      0
69      0       7
133     9       0
42      2       0
25      10      0
88      1       0
78      6       0
39      0       6
36      0       5
55      1       0
20      8       0
28      0       1
37      66      0
68      8       0
38      0       2
21      1       0
73      1       0
61      0       21
77      0       1
13      0       22
26      1       0
44      0       1
ADD COMMENTlink modified 2.6 years ago by h.mon31k • written 2.6 years ago by wang.yiguan30
0
gravatar for h.mon
2.6 years ago by
h.mon31k
Brazil
h.mon31k wrote:

Your post would be a good candidate for the google your own title™, had you used a more informative title. This is an axt alignment file.

But wait! No, this is an axtChain file:

The axtNet and axtChain alignments are produced by processing the alignment files with additional utilities written by Jim Kent at UCSC.

Probably used to transfer coordinates between genome versions.

P. S.: how you obtained the file? Did you ran the alignment yourself or downloaded the file? Which software you used to align the genomes? And so on... All this information is useful for answering questions, but alas, you didn't provide any of it.

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by h.mon31k

Thanks for providing me these info. I found some documents about the format: https://groups.google.com/a/soe.ucsc.edu/forum/#!topic/genome/m5eTOXbkiec http://genome.ucsc.edu/goldenPath/help/chain.html

ADD REPLYlink written 2.6 years ago by wang.yiguan30
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: 1932 users visited in the last hour