This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Running RepeatMasker with RepeatModeler database

Hello,

I have a draft assembly from Pacbio sequencing. I would like to evaluate the TEs abundance from it.

That's why I would like to use RepeatModeler and RepeatMasker. It does not exist a database corresponding to my species in RepeatMasker.

Should I run RepeatModeler on my draft assembly? And then, use the created database for RepeatMasker? In that way :

 BuildDatabase -name my_db -engine ncbi my_draft_assembly.fasta

 RepeatModeler -engine ncbi -pa 16 -database my_db

 RepeatMasker -s -xsmall -a -gff -pa 50 -u -lib  my_db.fa -dir final_RepeatMasker_out my_draft_assembly.fasta

I will also probably use REPET and/or Repbase to enrich my database. I will need to add the new sequences to the my_db.fa file, right?

Best

repeatmasker repeatmodeler
• 3,194 views • link
updated by Juke34 Juke34 • written by pablo pablo

1 answer

Should I run RepeatModeler on my draft assembly? And then, use the created database for RepeatMasker?

Yes, see here for more details on how to create de novo repeat ibrary using RepeatModeler: Create de novo repeat library

I will also probably use REPET and/or Repbase to enrich my database. I will need to add the new sequences to the my_db.fa file, right?

Yes you can do like that.

You can also consider using EDTA to create your repeat library.

• 0 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this