本文整理汇总了Python中Bio.Emboss.Applications.WaterCommandline.stdout方法的典型用法代码示例。如果您正苦于以下问题:Python WaterCommandline.stdout方法的具体用法?Python WaterCommandline.stdout怎么用?Python WaterCommandline.stdout使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Bio.Emboss.Applications.WaterCommandline
的用法示例。
在下文中一共展示了WaterCommandline.stdout方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: generate_water_cmd
# 需要导入模块: from Bio.Emboss.Applications import WaterCommandline [as 别名]
# 或者: from Bio.Emboss.Applications.WaterCommandline import stdout [as 别名]
def generate_water_cmd(macierz, pliki_fasta_rodzina):
"""
Generuje polecenia wywolania programu water EMBOSS dla wszystkich sekwencji podanych jako nazwy plikow je zawierajacych
:param macierz: lokalizacja/nazwa pliku z macierza substytucji PAM/BLOSUM
:param pliki_fasta_fodzina: lista lokalizacji/nazw plikow z sekwencjami bialkowymi fasta nalezacymi do danej rodziny
:return: polecenie wywolania programu water
"""
records = []
for file in pliki_fasta_rodzina:
handle = open(file, "rU")
records.extend(list(SeqIO.parse(handle, "fasta")))
handle.close()
from Bio.Emboss.Applications import WaterCommandline
all_water_cmd = []
for i in range(len(records)):
for j in range(len(records)):
if i < j:
water_cmd = WaterCommandline(gapopen=100, gapextend=10)#maksymalne wartosci aby uzyskac uliniowienia bezspacjowe
water_cmd.asequence = "asis:" + str(records[i].seq)
water_cmd.bsequence = "asis:" + str(records[j].seq)
water_cmd.stdout = True
water_cmd.sprotein=True
water_cmd.datafile=macierz
all_water_cmd.append(str(water_cmd))
return all_water_cmd