Dear Sir or Madam,

I'd like to aggregate wgEncodeDukeUniqueness35bp table with some gene tracks 
like Ensembl to generate whole-genome gene-level summary (like mean or median) 
of mappability. Could you tell me whether there is some simple way to do this? 

So far as I know, one solution is to transfer this track to galaxy and use its 
Datapoint aggregation tool. However, wgEncodeDukeUniqueness35bp track is too 
big to transfer. Thus, I guess what I can do is to download the raw file, 
wgEncodeDukeUniqueness35bp.wib (hgdownload.cse.ucsc.edu:/gbdb/hg18/wib/) and do 
some local coding. However, wib is binary not ascII format. I searched the mail 
list and found people uses hgWiggle to unpack wib format. But, when I ran 
"hgWiggle wgEncodeDukeUniqueness35bp.wib", it crashes with the error, "Line too 
long (more than 536870912 chars) line 1 of wgEncodeDukeUniqueness35bp.wib". 

Could you help me to handle this situation? 

Many thanks ahead.
                                
Best regards,
-------------------------------------------------------------
Yong Zhang
Ph.D, Research Scholar
Manyuan Long's Lab
University of Chicago

_______________________________________________
Genome maillist  -  [email protected]
https://lists.soe.ucsc.edu/mailman/listinfo/genome

Reply via email to