Mm9 multiple alignment
From genomewiki
Revision as of 19:36, 16 August 2007 by Hiram (talk | contribs) (→blastz alignment parameters and results)
To avoid artifacts in downstream processing of the UCSC multiple alignments, it is important to be careful on the use of the parameters used in the blastz processing pipeline. There are a number of steps in the pipeline and a variety of tunable parameters involved. This page will track the various parameters used in the alignments as they proceed toward the completion of a multiple alignment conservation track on the mm9 mouse (NCBI build 37) assembly
axtChain parameters and end results
| sequence | distance | axtChain minScore |
axtChain linearGap |
% of mm9 matched |
% of other matched by mm9 |
done |
|---|---|---|---|---|---|---|
| rat rn4 | 0.1587 | 3000 | medium | 68.357 | 69.541 | 16 August |
| human hg18 | 0.4667 | 3000 | medium | 38.499 | 35.201 | 16 August |
blastz alignment parameters details
| target | query | target size (overlap> |
query size (overlap> |
|---|---|---|---|
| mm9 | rat rn4 | 10M (10K) | 10M (0) |
| human hg18 | mm9 | 10M (0) | 10M (10K) |
matrix parameters
The "medium" gap score matrix, tuned for the mouse-human distance is:
tableSize 11 smallSize 111 position 1 2 3 11 111 2111 12111 32111 72111 152111 252111 qGap 350 425 450 600 900 2900 22900 57900 117900 217900 317900 tGap 350 425 450 600 900 2900 22900 57900 117900 217900 317900 bothGap 750 825 850 1000 1300 3300 23300 58300 118300 218300 318300
The "loose" gap score matrix, tuned for the chicken-human distance is:
tablesize 11 smallSize 111 position 1 2 3 11 111 2111 12111 32111 72111 152111 252111 qGap 325 360 400 450 600 1100 3600 7600 15600 31600 56600 tGap 325 360 400 450 600 1100 3600 7600 15600 31600 56600 bothGap 625 660 700 750 900 1400 4000 8000 16000 32000 57000
