本文整理汇总了Python中Bio.SeqFeature.SeqFeature.qualifiers["sequence_length"]方法的典型用法代码示例。如果您正苦于以下问题:Python SeqFeature.qualifiers["sequence_length"]方法的具体用法?Python SeqFeature.qualifiers["sequence_length"]怎么用?Python SeqFeature.qualifiers["sequence_length"]使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Bio.SeqFeature.SeqFeature
的用法示例。
在下文中一共展示了SeqFeature.qualifiers["sequence_length"]方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: SeqFeature
# 需要导入模块: from Bio.SeqFeature import SeqFeature [as 别名]
# 或者: from Bio.SeqFeature.SeqFeature import qualifiers["sequence_length"] [as 别名]
# original row as of nov23 read
#feature = SeqFeature(FeatureLocation(st-1,en), strand=framepart2,type="repeathit")
# start location should not be modified by -1 in this case
# changed Featurelocation start, older line:
#feature = SeqFeature(FeatureLocation(st,en), strand=framepart2,type="repeathit")
# Genbank magic adds one to start position?
feature = SeqFeature(FeatureLocation(st-1,en), strand=framepart2,type="repeathit")
#feature.id=rename
#feature.qualifiers["hit"]=rename
feature.id=isname
feature.qualifiers["hit"]=isname
feature.qualifiers["sequence_length"]=str(len(contigseq[myhsp.sbjct_start-1:myhsp.sbjct_end]))
#
feature.qualifiers["score"]=str(myhsp.score)
if myhsp.score>bestscorehsp.score:
bestscorehsp=myhsp
if myhsp.expect<bestexphsp.expect:
bestexphsp=myhsp
# hit length is disance from hits start plus hit end
feature.qualifiers["hitlength"]=str(1+max(st,en)-min(st,en))
#
feature.qualifiers["expect"]=str(myhsp.expect)
feature.qualifiers["query_start"]=str(myhsp.query_start)
feature.qualifiers["query_end"]=str(myhsp.query_end)
feature.qualifiers["sbjct_start"]=str(myhsp.sbjct_start)