percentage of repeat region every 500kb

[rashmijain]

I have gff3 file for repeat region, which shows repeat start and end position. I want to calculate what is % of repeat region in every 500 kb.
I have tried to find it on internet, but didn't find.

Any help is appreciated!

Thank you

Rashmi

[westerman]

I suspect that you will have to write some code to parse out the GFF and output what you want. This does not seem like a common request and thus I doubt if there is pre-existing code.

[rashmijain]

Thank you for the reply.....

[syfo]

I'd try with bedtools. Something like converting the gff to bed, generate a bed file with 500Kb genomic bins, run bedtools intersect or coverage to get the overlap in each bin and parse the resulting output to compute the precise values. If you have some knowledge on bioinformatics that should not be too problematic.

[rashmijain]

Hi Thank you for the reply

we have written code for that.