retrieve sequence alignment score produced by emboss in biopython -


i'm trying retrieve alignment score of 2 sequences compared using emboss in biopython. way know retrieve output text file produced emboss. problem there hundreds of these files iterate over. there easier/cleaner method retrieve alignment score, without resorting that? main part of code i'm using.

from bio.emboss.applications import stretchercommandline needle_cline = stretchercommandline(asequence=,bsequence=,gapopen=,gapextend=,outfile=) stdout, stderr = needle_cline() 

i had same problem , after time spent on searching neat solution popped white flag.

however, speed processing of output files did following things:

1) used re python module handling regular expressions extract data needed.

2) created ramdisk space output files. use of ramdisk here allowed processing , exchanging data in ram memory (much faster writing , reading output files hard drive, not mention saves hdd in case of processing massive number of alignments).


Popular posts from this blog

php - How should I create my API for mobile applications (Needs Authentication) -

5 Reasons to Blog Anonymously (and 5 Reasons Not To)

Google AdWords and AdSense - A Dynamic Small Business Marketing Duo