[Biopython-dev] Hmmer integration modules/package

Glen van Ginkel gvg500 at york.ac.uk
Thu Apr 7 08:10:29 EDT 2005


Hi Guys

As a project for my MRes in bioinformatics I have had to write a Python 
package that would expand Biopythons capabilities in interacting with useful 
programs. Here I submit code to interact with the Hmmer suite of programs. 

Basically, you instantiate the Hmmer object (Like a Hmmer commandline builder 
object), build a commandline, execute it and grab the results. Really simple 
for the user. This only works in UNIX since Hmmer is not available on Windows 
platforms and it obviously assumes you have Hmmer already installed on your 
PC. 
I have also tried to integrate the Bio.Align.Generic alignment object sothat 
the Hmmer object is able to handle alignment objects by writing the records 
of the alignment object to a temporary fasta file. 

Since I have to write this up as a report I would greatly appreciate any 
criticism of any kind and perhaps some suggestions as to how I might improve 
the code. Also, If you can think of any other ways to implement the Hmmer 
stuff please let me know. 

At the moment I am working on a test suite for the project. If you would like 
the code I have written to exercise some of the methods please let me know 
and I'll send it over. 
I also have a bit of documentation I would like to add to the Application 
package because I feel if failed to help me concerning certain aspects. How 
would I go about this?

I look forward to hearing your suggestions. 

Glen van Ginkel

Attatched is the code as well as an example of a test suite using the files 
from Eddy2003 and an experiment file experiment.fas
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Applications.py
Type: application/x-python
Size: 16558 bytes
Desc: not available
Url : http://portal.open-bio.org/pipermail/biopython-dev/attachments/20050407/aab5b8bd/Applications-0001.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: __init__.py
Type: application/x-python
Size: 320 bytes
Desc: not available
Url : http://portal.open-bio.org/pipermail/biopython-dev/attachments/20050407/aab5b8bd/__init__-0001.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: HmmerStandalone.py
Type: application/x-python
Size: 9139 bytes
Desc: not available
Url : http://portal.open-bio.org/pipermail/biopython-dev/attachments/20050407/aab5b8bd/HmmerStandalone-0001.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testsuite.py
Type: application/x-python
Size: 2427 bytes
Desc: not available
Url : http://portal.open-bio.org/pipermail/biopython-dev/attachments/20050407/aab5b8bd/testsuite-0001.bin
-------------- next part --------------
>S134211|S134211|GLOBIN - BRINE SHRIMP1
DKATIKRTWATVTDLPSFGRNVFLSVFAAK
>S134212|S134212|GLOBIN - BRINE SHRIMP2
PEYKNLFVEFRNIPASELASSERLLYHGGR
>S134213|S134213|GLOBIN - BRINE SHRIMP3
VLSSIDEAIAGIDTPDRAVKTLLALGERHI
>S134214|S134214|GLOBIN - BRINE SHRIMP4
SRGTVRRHFEAFSYAFIDELKQRGVESADL
>S134215|S134215|GLOBIN - BRINE SHRIMP5
AAWRRGWDNIVNVLEAGLLRRQIDLEVTGL
>S134216|S134216|GLOBIN - BRINE SHRIMP6
SCVDVANIQESWSKVSGDLKTTGSVVFQRM
>S134217|S134217|GLOBIN - BRINE SHRIMP7
INGHPEYQQLFRQFRDVDLDKLGESNSFVA


More information about the Biopython-dev mailing list