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).