[Genome] hgWiggle Phastcons

[Pai, Deepa]

Hi,

 

     Given a set of genomic coordinates, I want to calculate the
conservation score for the region.

      Similar question has been asked before, but I havent been
successful using the hgWiggle.

 

       For hg18,

       I have downloaded and unzipped wgEncodeUcsdNgTaf1Signal.wigVar.

       Also, downloaded phastCons17way.wib to the same folder.

        $ ln -s wgEncodeUcsdNgTaf1Signal.wigVar
wgEncodeUcsdNgTaf1Signal.wig

 

        Now,  when I run hgWiggle, on the flat files, I get the error
message,

         [paid at moria wig]$ ~/bin/i386/hgWiggle
-position=chr6:26195715-26195740 wgEncodeUcsdNgTaf1Signal.wig

        #       no values found with these constraints

        [paid at moria wig]$ ~/bin/i386/hgWiggle -chr=chr12
wgEncodeUcsdNgTaf1Signal.wig

        #       no values found with these constraints

 

        Is it because the .wib and .wig file names are different? How
can I link them?

 

Thanks

Deepa.