[Bioperl-l] automated stand alone blast with repeat masker

Dave Messina David.Messina at sbc.su.se
Mon May 12 17:58:40 EDT 2008


I haven't done this myself, but from a quick search on the BioPerl website,
it looks like you'll want to use the
Bio::Tools::Run::RepeatMasker<http://doc.bioperl.org/releases/bioperl-current/bioperl-run/Bio/Tools/Run/RepeatMasker.html>module
to create a repeat-masked fasta file.

If you RepeatMask your query sequence(s), then you need to specify that
sequence when you create your Bio::Seq object.

If you instead RepeatMask your database, you'll need to create a blast
database from the repeat-masked sequences and specify that db in your
@params. I don't think there's a module for running formatdb, but you can do
it through a system call.



Dave


More information about the Bioperl-l mailing list