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!
• 564 views
•
link
0 answers
No answers yet.
Log in to answer this question.
maybe you could try to write a script to do it yourself