This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to obtain and compare genomic coordinates from RepeatMasker output?

Hello folks!

I want to mask a genome with a particular repeat library using RepeatMasker.

Then I want to cross the coordinates of the repeats with those of gene annotations to find overlaps between them and study associations and stuff.

I'm only starting to consider feasible ways to do that so any input would be great.

Thanks!

repeats repeatmasker genome browser gbrowse
• 4,336 views • link
updated by Ram Ram • written by gbdias gbdias

2 answers

Sounds like a job for BedTools. You should be able to make a GFF from the RepeatMasker output and use bedtools intersect to find overlaps.

• 0 views • link
updated by Ram Ram • written by Jon Jon

Thank you very much! I reading BEDTools' documentation and it really seems to be the right tool for this task.

• 1 views • link
written by gbdias gbdias

The answer by Jon is a good one to get overlaps, and I can tell you how to make the GFF. In the 'util' directory of the RepeatMasker distribution there is a script called 'rmOutToGFF3.pl' that will do the conversion. The usage is pretty simple since it writes to stdout:

perl rmOutToGFF3.pl my_repeatmasker.out > my_repeatmasker.gff
• 0 views • link
written by SES SES

Thank you very much! This will help a lot.

• 1 views • link
written by gbdias gbdias

Log in to answer this question.

More posts like this