js DNA动态序列比对代码

<title>动态序列比对</title> 
<img src="/_New/_cj/myeditor/sysimage/space.gif" class="ewebeditor__script" _ewebeditor_fake_tag="script" _ewebeditor_fake_value="%3Cscript%3E%20%0D%0Afunction%20pairaln(seq1%2C%20seq2)%7B%20%0D%0A%2F%2F%20scoring%20scheme%20%0D%0A%0D%0Avar%20MATCH%20%3D%201%3B%20%2F%2F%20%2B1%20for%20letters%20that%20match%20%0D%0Avar%20MISMATCH%20%3D%20-1%3B%20%2F%2F%20-1%20for%20letters%20that%20mismatch%20%0D%0Avar%20GAP%20%3D%20-1%3B%20%2F%2F%20-1%20for%20any%20gap%20%0D%0A%2F%2F%20initialization%20%0D%0A%0D%0Avar%20score%3D%5B%5B%5D%5D%3B%20%0D%0Avar%20pointer%3D%5B%5B%5D%5D%3B%20%0D%0Ascore%5B0%5D%5B0%5D%20%3D%200%3B%20%0D%0Apointer%5B0%5D%5B0%5D%20%3D%200%3B%20%0D%0Avar%20l1%3Dseq1.length%20%0D%0Avar%20l2%3Dseq2.length%20%0D%0A%0D%0Afor(var%20j%20%3D%201%3B%20j%20%3C%3D%20l1%3B%20j%2B%2B)%20%7B%20%0D%0A%0D%0Ascore%5B0%5D%5Bj%5D%20%3D%20GAP%20*%20j%3B%20%0D%0Apointer%5B0%5D%5Bj%5D%20%3D%202%3B%20%0D%0A%7D%20%0D%0Afor%20(var%20i%20%3D%201%3B%20i%20%3C%3D%20l2%3B%20i%2B%2B)%20%7B%20%0D%0A%0D%0Ascore%5Bi%5D%3D%5B%5D%3B%20%0D%0Apointer%5Bi%5D%3D%5B%5D%3B%20%0D%0Ascore%5Bi%5D%5B0%5D%20%3D%20GAP%20*%20i%3B%20%0D%0Apointer%5Bi%5D%5B0%5D%20%3D%203%3B%20%0D%0A%7D%20%0D%0A%2F%2F%20fill%20%0D%0A%0D%0Afor(var%20i%20%3D%201%3B%20i%20%3C%3D%20l2%3B%20i%2B%2B)%20%7B%20%0D%0Avar%20letter2%20%3D%20seq2.charAt(i-1)%3B%20%0D%0A%0D%0Afor(var%20j%20%3D%201%3B%20j%20%3C%3D%20l1%3B%20j%2B%2B)%20%7B%20%0D%0A%0D%0Avar%20diagonal_score%2C%20left_score%2C%20up_score%3B%20%0D%0A%2F%2F%20calculate%20match%20score%20%0D%0A%0D%0Avar%20letter1%20%3D%20seq1.charAt(j-1)%3B%20%0D%0A%0D%0Aif%20(letter1%20%3D%3D%20letter2)%20%7B%20%0D%0A%0D%0Adiagonal_score%20%3D%20score%5Bi-1%5D%5Bj-1%5D%20%2B%20MATCH%3B%20%0D%0A%7D%20%0D%0A%0D%0Aelse%20%7B%20%0D%0Adiagonal_score%20%3D%20score%5Bi-1%5D%5Bj-1%5D%20%2B%20MISMATCH%3B%20%0D%0A%7D%20%0D%0A%2F%2F%20calculate%20gap%20scores%20%0D%0A%0D%0Aup_score%20%3D%20score%5Bi-1%5D%5Bj%5D%20%2B%20GAP%3B%20%0D%0Aleft_score%20%3D%20score%5Bi%5D%5Bj-1%5D%20%2B%20GAP%3B%20%0D%0A%2F%2F%20choose%20best%20score%20%0D%0A%0D%0Aif%20(diagonal_score%20%3E%3D%20up_score)%20%7B%20%0D%0A%0D%0Aif%20(diagonal_score%20%3E%3D%20left_score)%20%7B%20%0D%0A%0D%0Ascore%5Bi%5D%5Bj%5D%20%3D%20diagonal_score%3B%20%0D%0Apointer%5Bi%5D%5Bj%5D%20%3D%201%3B%20%0D%0A%7D%20%0D%0A%0D%0Aelse%20%7B%20%0D%0Ascore%5Bi%5D%5Bj%5D%20%3D%20left_score%3B%20%0D%0Apointer%5Bi%5D%5Bj%5D%20%3D%202%3B%20%0D%0A%7D%20%0D%0A%0D%0A%7D%20else%20%7B%20%0D%0A%0D%0Aif%20(up_score%20%3E%3D%20left_score)%20%7B%20%0D%0A%0D%0Ascore%5Bi%5D%5Bj%5D%20%3D%20up_score%3B%20%0D%0Apointer%5Bi%5D%5Bj%5D%20%3D%203%3B%20%0D%0A%7D%20%0D%0A%0D%0Aelse%20%7B%20%0D%0Ascore%5Bi%5D%5Bj%5D%20%3D%20left_score%3B%20%0D%0Apointer%5Bi%5D%5Bj%5D%20%3D%202%3B%20%0D%0A%7D%20%0D%0A%7D%20%0D%0A%7D%20%0D%0A%7D%20%0D%0A%2F%2F%20trace-back%20%0D%0A%0D%0Avar%20align1%20%3D%20%5B%5D%3B%20%0D%0Avar%20align2%20%3D%20%5B%5D%3B%20%0D%0A%2F%2F%20start%20at%20last%20cell%20of%20matrix%20%0D%0Avar%20j%20%3D%20l1%3B%20%0D%0Avar%20i%20%3D%20l2%3B%20%0D%0Awhile%20(1)%20%7B%20%0D%0A%0D%0Aif(pointer%5Bi%5D%5Bj%5D%20%3D%3D%200)break%20%2F%2F%20ends%20at%20first%20cell%20of%20matrix%20%0D%0Aif%20(pointer%5Bi%5D%5Bj%5D%20%3D%3D%201)%20%7B%20%0D%0A%0D%0Aalign1%5Balign1.length%5D%20%3D%20seq1.charAt(j-1)%3B%20%0D%0Aalign2%5Balign2.length%5D%20%3D%20seq2.charAt(i-1)%3B%20%0D%0Ai--%3B%20%0D%0Aj--%3B%20%0D%0A%7D%20%0D%0A%0D%0Aif%20(pointer%5Bi%5D%5Bj%5D%20%3D%3D%202)%20%7B%20%0D%0A%0D%0Aalign1%5Balign1.length%5D%20%3D%20seq1.charAt(j-1)%3B%20%0D%0Aalign2%5Balign2.length%5D%20%3D%20%22-%22%3B%20%0D%0Aj--%3B%20%0D%0A%7D%20%0D%0A%0D%0Aif%20(pointer%5Bi%5D%5Bj%5D%20%3D%3D%203)%20%7B%20%0D%0A%0D%0Aalign1%5Balign1.length%5D%20%3D%20%22-%22%3B%20%0D%0Aalign2%5Balign2.length%5D%20%3D%20seq2.charAt(i-1)%3B%20%0D%0Ai--%3B%20%0D%0A%7D%20%0D%0A%7D%20%0D%0Aalign1%20%3D%20align1.reverse()%3B%20%0D%0Aalign2%20%3D%20align2.reverse()%3B%20%0D%0Areturn%20%22%E5%BA%8F%E5%88%97%EF%BC%9A%3Cp%3E%22%2Bseq1%2B%22%3Cbr%3E%22%2Bseq2%2B%22%3Cp%3E%E6%AF%94%E5%AF%B9%E7%BB%93%E6%9E%9C%EF%BC%9A%3Ctable%3E%3Ctr%3E%3Ctd%3E%22%2Balign1.join(%22%3C%2Ftd%3E%3Ctd%3E%22)%2B%22%3C%2Ftd%3E%3C%2Ftr%3E%3Ctr%3E%3Ctd%3E%22%2Balign2.join(%22%3C%2Ftd%3E%3Ctd%3E%22)%2B%22%3C%2Ftr%3E%3C%2Ftable%3E%22%20%0D%0A%7D%20%0D%0Adocument.write(pairaln('CTGGGCTGACTGA'%2C%20'GACTAGCTAGACTGA'))%20%0D%0A%3C%2Fscript%3E">

　　 [Ctrl+A 全选注:如需引入外部Js需刷新才能执行]

教程首页更多教程