本文整理汇总了Python中paleomix.common.formats.msa.MSA.validate方法的典型用法代码示例。如果您正苦于以下问题:Python MSA.validate方法的具体用法?Python MSA.validate怎么用?Python MSA.validate使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类paleomix.common.formats.msa.MSA
的用法示例。
在下文中一共展示了MSA.validate方法的4个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: _read_sequences
# 需要导入模块: from paleomix.common.formats.msa import MSA [as 别名]
# 或者: from paleomix.common.formats.msa.MSA import validate [as 别名]
def _read_sequences(filenames):
results = {}
for filename in filenames:
results[filename] = MSA.from_file(filename)
MSA.validate(*results.values())
return results.iteritems()
示例2: interleaved_phy
# 需要导入模块: from paleomix.common.formats.msa import MSA [as 别名]
# 或者: from paleomix.common.formats.msa.MSA import validate [as 别名]
def interleaved_phy(msa, add_flag = False, max_name_length = _MAX_NAME_LENGTH):
MSA.validate(msa)
header = "%i %i" % (len(msa), msa.seqlen())
if add_flag:
header += " I"
result = [header, ""]
padded_len = min(max_name_length, max(len(name) for name in msa.names())) + 2
padded_len -= padded_len % -(_BLOCK_SIZE + _BLOCK_SPACING) + _BLOCK_SPACING
streams = []
spacing = " " * _BLOCK_SPACING
for record in sorted(msa):
name = record.name[:max_name_length]
padding = (padded_len - len(name)) * " "
lines = []
line = [name, padding]
for block in grouper(_BLOCK_SIZE, record.sequence, fillvalue = ""):
block = "".join(block)
if sum(len(segment) for segment in line) >= _LINE_SIZE:
lines.append("".join(line))
line = [block]
else:
line.extend((spacing, block))
lines.append("".join(line))
streams.append(lines)
for rows in zip(*streams):
result.extend(row for row in rows)
result.append("")
result.pop()
return "\n".join(result)
示例3: _run
# 需要导入模块: from paleomix.common.formats.msa import MSA [as 别名]
# 或者: from paleomix.common.formats.msa.MSA import validate [as 别名]
def _run(self, _config, temp):
# Read and check that MSAs share groups
msas = [MSA.from_file(filename) for filename in sorted(self.input_files)]
MSA.validate(*msas)
blocks = []
for msa in msas:
blocks.append(sequential_phy(msa, add_flag = self._add_flag))
with open(reroot_path(temp, self._out_phy), "w") as output:
output.write("\n\n".join(blocks))
示例4: sequential_phy
# 需要导入模块: from paleomix.common.formats.msa import MSA [as 别名]
# 或者: from paleomix.common.formats.msa.MSA import validate [as 别名]
def sequential_phy(msa, add_flag = False, max_name_length = _MAX_NAME_LENGTH):
MSA.validate(msa)
header = "%i %i" % (len(msa), msa.seqlen())
if add_flag:
header += " S"
spacing = " " * _BLOCK_SPACING
result = [header, ""]
for record in sorted(msa):
result.append(record.name[:max_name_length])
blocks = grouper(_BLOCK_SIZE, record.sequence, fillvalue = "")
lines = grouper(_NUM_BLOCKS, blocks)
for line in lines:
result.append(spacing.join("".join(block) for block in line if block))
return "\n".join(result)