Chromosome Position In The Ucsc Phylop46Way.Placental.Wigfix File
3
0
Entering edit mode
11.2 years ago
Abdel ▴ 150

Hi everyone
I downloaded the chr1.phyloP46way.placental.wigFix file from ucsc , i need it to get a phyloP score for a given position of the chr1 . My problem is that there is no position field in this file , just the scores...
0.056
0.056
0.056
0.064
0.056
0.064
0.056
. . .

Did someone know how am i supposed to figure out a chromosomic position ?

wiggle format ucsc • 4.9k views
2
Entering edit mode
11.2 years ago

You have a fixed step wiggle file. Check the Wiggle specification:

http://genome.ucsc.edu/goldenPath/help/wiggle.html

0
Entering edit mode

if i'm not mistaken , the phyloP score for the chr1 10918 is 0.064 it is intuitive , the confusing thing is that when querying ucsc table browser for this position ,i found a different value 0.0400079 (clade :mammal ,Group :comparative genomics ,track :conservation ,table :mammals Cons phyloP46wayPlacental) is that normal ?

2
Entering edit mode
11.2 years ago

Istvan is right and you missed the very first line of the wig file in your input: the wig starts at 10918 on chr1 and the distance between each value is 1. See the wig spec.

curl -s "ftp://biomirror.auckland.ac.nz/goldenpath/currentGenomes/Homo_sapiens/phyloP46way/placentalMammals/chr1.phyloP46way.placental.wigFix.gz" |\
gunzip -c |\

fixedStep chrom=chr1 start=10918 step=1
0.064
0.056
0.064
0.056
0.064
0.064
0.064
0.064
0.064

0
Entering edit mode

if i'm not mistaken , the phyloP score for the chr1 10918 is 0.064 it is intuitive , the confusing thing is that when querying ucsc table browser for this position ,i found a different value 0.0400079 (clade :mammal ,Group :comparative genomics ,track :conservation ,table :mammals Cons phyloP46wayPlacental) is that normal ?

2
Entering edit mode
10.9 years ago
Emma ▴ 140

I have been having the same problem, i.e. I am getting different scores between the wigfix files (chr21.phastCons46way.placental.wigFix) and the galaxy platform (Placental Mammal Conservation by PhastCons), but for the phastCons score (hg19 in both databases):

fixedStep.chrom.chr21.start.9411194.step.1
0.053
0.044
0.033
0.017
0.010
0.008


and from galaxy:

variableStep chrom=chr21 span=1
9411194 0.0493701
9411195 0.0394961
9411196 0.029622
9411197 0.014811
9411198 0.00987402
9411199 0.00493701
9411200 0.00987402


I would be grateful if you mentioned how you solved your problem and what I could be missing.

0
Entering edit mode

Please don't post a new question in the answer section of an old question. Post a new standalone question if you have one that way you can get it answered.