本文整理汇总了Python中Bio.Emboss.Applications.WaterCommandline.asequence方法的典型用法代码示例。如果您正苦于以下问题：Python WaterCommandline.asequence方法的具体用法？Python WaterCommandline.asequence怎么用？Python WaterCommandline.asequence使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Bio.Emboss.Applications.WaterCommandline的用法示例。
在下文中一共展示了WaterCommandline.asequence方法的4个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: generate_water_cmd
# 需要导入模块: from Bio.Emboss.Applications import WaterCommandline [as 别名]
# 或者: from Bio.Emboss.Applications.WaterCommandline import asequence [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示例2: doWater
# 需要导入模块: from Bio.Emboss.Applications import WaterCommandline [as 别名]
# 或者: from Bio.Emboss.Applications.WaterCommandline import asequence [as 别名]
def doWater(contig, seq):
with open("contig.faa", "w") as stuff1:
stuff1.write(">contig\n")
stuff1.write(contig)
with open("seq.faa", "w") as stuff2:
stuff2.write(">seq\n")
stuff2.write(str(seq))
water_cline = WaterCommandline()
water_cline.asequence="contig.faa"
water_cline.bsequence="seq.faa"
water_cline.gapopen=10
water_cline.gapextend=0.5
water_cline.outfile="water.txt"
stdout, stderr = water_cline()
print(stdout + stderr)
values = getStartEnd()
return values示例3: emboss_local_pairwise_alignment
# 需要导入模块: from Bio.Emboss.Applications import WaterCommandline [as 别名]
# 或者: from Bio.Emboss.Applications.WaterCommandline import asequence [as 别名]
def emboss_local_pairwise_alignment(query_dir, seq_type):
if seq_type == 'fg':
print '\n ...pairwise comparison of functional gene sequences...\n'
elif seq_type == 'ssu':
print '\n ...pairwise comparison of SSU rRNA sequences...\n'
water_cline = WaterCommandline()
water_cline.gapopen=10
water_cline.gapextend=0.5
query_list = [query for query in sorted(glob.glob(query_dir+"/*.fa"))]
for i, a_seq in enumerate(query_list):
water_cline.asequence=str(a_seq)
for j, b_seq in enumerate(query_list[i:]):
water_cline.bsequence=str(b_seq)
align_out = query_dir+"/pairwise_"+str(i+1)+"_"+str(i+j+1)+".aln"
water_cline.outfile=str(align_out)
water_cline()
print 'Done\n'
return query_dir+"/*.aln"示例4: open
# 需要导入模块: from Bio.Emboss.Applications import WaterCommandline [as 别名]
# 或者: from Bio.Emboss.Applications.WaterCommandline import asequence [as 别名]
# http://rosalind.info/problems/swat/
from Bio.Emboss.Applications import WaterCommandline
from Bio import ExPASy, SeqIO
if __name__ == "__main__":
ids = open('rosalind_swat.txt').read().split(' ')
for i in ids:
handle = ExPASy.get_sprot_raw(i)
r = SeqIO.read(handle, "swiss")
handle.close()
with open(i, 'w') as f:
SeqIO.write(r, f, 'fasta')
water_cline = WaterCommandline()
water_cline.asequence = ids[0]
water_cline.bsequence = ids[1]
water_cline.outfile = "rosalind_swat_output.txt"
water_cline.gapopen = 10
water_cline.gapextend = 1
water_cline()
for line in open('rosalind_swat_output.txt').readlines():
if 'Score:' in line:
print(int(float(line[:-1].split(':')[-1].strip())))