Shell Programming and Scripting

Dear ALL,
I would really appreciate if you could help me in reformatting a file in this way:

The file refers to a list of genetic coordinates, each lines has a score value and the associated chromosome is listed in the line starting with chrom .

If more coordinates are found, the start refers to the start plus n lines:

original file

Code:

chrom=chr1 start=3000306
0.006
0.010
0.014
chrom=chrX start=40000306
0.014
chrom=chr9 start=80000306 
0.1
0.2

I would like to obtain a file like this

Code:

#chr #start #end #score

chr1 3000306 3000307 0.006
chr1 3000307 3000308 0.010
chr1 3000308 3000309 0.019
chrX 40000306 40000307 0.014
chr9 80000306 80000307 0.1
chr9 80000307 80000308 0.2

Let me know if you can help me, the original file in up to 1 GB and is impossible editing my hand,
Thanks,
Paolo

Try this,

Code:

awk '/^chrom/{split($1,a,"=");split($2,b,"=");next} { printf "%s\t%10d\t%10d\t%f\n",a[2],b[2],b[2]+1,$1;b[2]++}' filename

Amazing thanks again for all your support!