Shell Programming and Scripting

Dear All,
I would like some help to process the following file

Code:

ID	n_exons	start	end
Zfp85-rs1	4	67747799,67751609,67752387,67755062,	67749722,67751705,67752514,67755134,
Taf9b	7	106206873,106209706,106216126,106218096,106218309,106218895,106219578,	106208657,106209817,106216196,106218231,106218446,106218977,106219842,
Zfp458	4	67254917,67259560,67260116,67269030,	67258144,67259656,67260243,67269068,
Zfp941	3	140809676,140818587,140822034,	140813314,140818714,140822178,

-first column is the gene name
-second is the number of exons
-third is the exons start cordinates comma separeted
-fourth is the exon end cordinates comma separated

I would like the following output
-first column gene name
-second column number of exons
-third column length

length should be the sum of the difference between the n cordinate in the fourth column with the n cordinate in the third column, example for gene name Zfp941

Code:

1st exon end	140813314
1st exon start	140809676
diff	3638
	
2nd exon end	140818714
2nd exon start	140818587
diff	127
	
3th exon end	140822178
3th exon start	140822034
diff	144
	
length	3909

Any help? The file is really big, almost 40000 row, like the number of genes in our genome,
Many Thanks,
Paolo

That doesn't seem too difficult a problem. Any attempts/ideas/thoughts from your side?

Thank you for your fine help,
Cheers!