如何从BioAlignments.jl比对中检索比对区域的序列索引？

问题描述

我想访问从BioAlignments.jl中的Julia对对函数返回的比对数据的索引，以了解在原始序列的上下文中比对发生的位置。

using BioAlignments
using BioSequences
scoremodel = AffineGapScoreModel(EDNAFULL,gap_open=-5,gap_extend=-1);
my_alignment = pairalign(LocalAlignment(),dna"ATATTAGGTATTGATTATTGTACGCGGCCCGGC",dna"TTGATTATTGT",scoremodel)
alignment(my_alignment)

例如，这样的脚本将输出一个对齐对象，我可以通过score（）函数从中访问得分。但是，我想知道在作为输入提供的原始序列中哪里发生了对齐，并且知道如何调用存储此索引的变量。在文档的任何地方似乎都找不到。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

bioinformatics julia