[Bioperl-l] Slow restriction digestion

Shawn Hoon shawnh at fugu-sg.org
Mon May 12 01:11:16 EDT 2003


On Tuesday, May 13, 2003, at 12:02  AM, anthony.underwood at hpa.org.uk  
wrote:

> Hi All,
>
>
> I'm trying to list all the MseI sites (recognition site TTAA)  in a 4
> megabase genome, using the Bio::Tools::RestrictionEnzyme module.
> This works but takes a VERY long time. Is there a way of speeding this  
> up?.
> I realise it maybe a perl string searching problem- since searching  
> long
> strings in perl tends to be slow and when I've written perl code to do  
> the
> same without Bio::Perl it is also very slow. Just wondered if there  
> were any
> neat tricks?
>
>


U might want to avoid using perl check out EMBOSS apps

http://www.hgmp.mrc.ac.uk/Software/EMBOSS/Apps/ 
nucleic_restriction_group.html


hope that helps..

shawn


> Thanks,
>
> Anthony
>
> Dr Anthony Underwood
> Bioinformatics Unit
> Central Public Health Laboratory
> Health Protection Agency
> 61 Colindale Avenue
> London
> NW9 5HT
> t:    0208 2004400 ext. 3618
> f:    0208 3583138
> e: anthony.underwood at hpa.org.uk
>
> *********************************************************************** 
> ***
> The information contained in the EMail and any attachments is
> confidential
> and intended solely and for the attention and use of the named
> addressee(s).
> It may not be disclosed to any other person without the express
> authority of
> the HPA, or the intended recipient, or both. If you are not the
> intended
> recipient, you must not disclose, copy, distribute or retain this
> message or
> any part of it.
>
> For information on how to send data to the HPA in encrypted form via
> E.Mail, visit www.HPA.org.uk.
>
> This footnote also confirms that this EMail has been swept for computer
> viruses, but please re-sweep any attachments before opening or saving.
>
> HTTP://www.HPA.org.uk
> *********************************************************************** 
> ***
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at bioperl.org
> http://pw600a.bioperl.org/mailman/listinfo/bioperl-l
>



More information about the Bioperl-l mailing list