1-场景: 有100个序列A1、A2…A100,需要比较序列俩俩之间的identity,以保持序列的不冗余程度
2-用jalview进行比对
一对一对地比较则有 100*(100+1)/2=5050对,能比较到山无棱天地合的时候
可以用软件jalview(免费软件、官网download即可)进行比对。
1-打开fasta/aln文件
2-选中所有序列,然后点击caculate-pairwase alignment进行比对
3-这将生成所有选中序列的排列组合的pairwise alignment文件
当序列不多的时候,之间读取Percentage ID即可,但是当序列多的时候,挨个读取选择非常麻烦。
用python读取jalview的pairwise alignment文件
path=