我正在尝试根据Python中的全局对齐算法(或Needleman-Wunsch算法)绘制矩阵.
我不知道matplotlib是否是这种情况下最好的工具.我试图使用Bokeh,但结构很难适合我想要的矩阵.
我正在使用Bio.SeqIO(BioPython的标准序列输入/输出接口)来存储两个序列.
我得到的结果与此图像类似:
这可能在Matplotlib?我怎样才能做到这一点?
UPDATE
最后,我能够根据ImportanceOfBeingErnest给出的答案构建算法.结果如下:
以下是此实现的要点:plot_needleman_wunsch.py
这是整个项目(正在进行中):bmc-sequence-alignment