This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Coverage file conversion to fixed step wig

I'm having some file conversion problem in unix. I would like to create a fixed step wig file from a coverageBed output.

coverageBed output example (tab-separated):

chr1    0 250     2       63      250     0.2520000
chr1    250 500     2       62      250     0.2480000
chr1    500 750     1       51      250     0.2040000
n
chr2 0   250  1   51 250 0.2040000
m

Desired output for example:

track type=wiggle_0 name="fixedStep" description="fixedStep format" visibility=full autoScale=off viewLimits=0:1000 color=0,200,100 maxHeightPixels=100:50:20 graphType=points priority=20
fixedStep chrom=chr1 start=0 step=250 span=250
2
2
1
n
fixedStep chrom=chr2 start=0 step=250 span=250
1
m

Maybe rtracklayer can do this? Thank you for your help!

format fixed wig coverage

maybe you could try to write a script to do it yourself

0 answers

No answers yet.

Log in to answer this question.