本文整理汇总了Python中Bio.SeqFeature.SeqFeature.qualifiers['ref']方法的典型用法代码示例。如果您正苦于以下问题：Python SeqFeature.qualifiers['ref']方法的具体用法？Python SeqFeature.qualifiers['ref']怎么用？Python SeqFeature.qualifiers['ref']使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Bio.SeqFeature.SeqFeature的用法示例。
在下文中一共展示了SeqFeature.qualifiers['ref']方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: main
# 需要导入模块: from Bio.SeqFeature import SeqFeature [as 别名]
# 或者: from Bio.SeqFeature.SeqFeature import qualifiers['ref'] [as 别名]
def main():
genome_record = SeqIO.read(INPUT_GENOME, 'genbank')
with open(SNP_CSV_DATA_FILE) as input_fh:
reader = csv.DictReader(input_fh, SNP_FIELD_NAMES)
reader.next() # Ignore header.
for row in reader:
feature_ref = row['ref']
feature_alt = row['alt'].replace('[', '').replace(']', '')
feature_start = int(row['position']) - 1 # pythonic
feature_end = feature_start + len(feature_alt)
feature_location = FeatureLocation(feature_start, feature_end)
feature = SeqFeature(location=feature_location,
type=REALIGNED_SNP_TYPE,
strand=1)
feature.qualifiers['ref'] = row['ref']
feature.qualifiers['alt'] = row['alt']
genome_record.features.append(feature)
with open(OUTPUT_GENOME, 'w') as output_fh:
SeqIO.write(genome_record, output_fh, 'genbank')