欧美三区_成人在线免费观看视频_欧美极品少妇xxxxⅹ免费视频_a级毛片免费播放_鲁一鲁中文字幕久久_亚洲一级特黄

使用python實現BLAST

系統(tǒng) 1651 0

最近在自學python,又用python實現了一下BLAST。

這次更新了打分函數如下,空位罰分改為-5,但不區(qū)分gap open 和 gap extend。

使用python實現BLAST_第1張圖片

            
''''' 
@author: JiuYu 
''' 
 
def score(a,b):#scoring function 
  score=0 
  lst=['AC','GT','CA','TG'] 
  if a==b: 
    score +=2 
  elif a+b in lst: 
    score += -5 
  else: 
    score += -7 
  return score 
 
def BLAST(seq1,seq2):#Basic Local Alignment Search Tool 
  l1 = len(seq1) 
  l2 = len(seq2) 
  GAP =-5   #-5 for any gap 
  scores =[] 
  point =[] 
   
  for j in range(l2+1): 
    if j == 0: 
      line1=[0] 
      line2=[0] 
      for i in range(1,l1+1): 
        line1.append(GAP*i) 
        line2.append(2) 
    else: 
      line1=[] 
      line2=[] 
      line1.append(GAP*j) 
      line2.append(3) 
    scores.append(line1) 
    point.append(line2) 
   
  #fill the blank of scores and point 
  for j in range(1,l2+1): 
    letter2 = seq2[j-1] 
    for i in range(1,l1+1): 
      letter1 = seq1[i-1] 
      diagonal_score = score(letter1, letter2) + scores[j-1][i-1] 
      left_score = GAP + scores[j][i-1] 
      up_score = GAP + scores[j-1][i] 
      max_score = max(diagonal_score, left_score, up_score) 
      scores[j].append(max_score) 
       
      if scores[j][i] == diagonal_score: 
        point[j].append(1) 
      elif scores[j][i] == left_score: 
        point[j].append(2) 
      else: 
        point[j].append(3) 
         
  #trace back 
  alignment1='' 
  alignment2='' 
  i = l2 
  j = l1 
  print 'scores =',scores[i][j] 
  while True: 
    if point[i][j] == 0: 
      break 
    elif point[i][j] == 1: 
      alignment1 += seq1[j-1] 
      alignment2 += seq2[i-1] 
      i -= 1 
      j -= 1 
    elif point[i][j] == 2: 
      alignment1 += seq1[j-1] 
      alignment2 += '-' 
      j -= 1 
    else: 
      alignment1 += '-' 
      alignment2 += seq2[i-1] 
      i -= 1 
       
  #reverse alignment 
  alignment1 = alignment1[::-1] 
  alignment2 = alignment2[::-1] 
  print 'The best alignment:' 
  print alignment1 
  print alignment2 
 
seq1=raw_input('Please input your first sequences:\n') 
seq2=raw_input('input second sequences:\n') 
BLAST(seq1, seq2) 
          

運行結果:

使用python實現BLAST_第2張圖片

無疑python對字符串的處理更加強大,語言也更加簡單,優(yōu)雅。比如最后逆序輸出alignment,java我是單獨寫了一個逆序函數,而python只用一個語句就可以完成相同任務。

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦?。。?/p>

發(fā)表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 久久久久久久久成人 | 精品a在线观看 | 秦浩最新电视剧 | 亚洲欧美日韩综合二区三区 | 精品久久久久久久久久久 | www.亚洲在线 | 国产精品综合色区小说 | 国产福利不卡一区二区三区 | 中文字幕av在线播放 | 99福利视频 | 国产偷国产偷亚洲高清在线 | 免费国产一级特黄久久 | 午夜性啪啪A片免费播放 | 欧美日韩一区视频 | 免费观看欧美一级片 | 久久se精品一区精品二区 | 婷婷人人爽人人做人人添 | 午夜福利国产在线观看1 | 欧美福利 | 黄色亚洲视频 | 九一传媒在线观看 | 小视频网址 | 久久国产精品久久 | 久久秋霞理论电影 | 久久永久免费中文字幕 | 最近免费日本视频在线 | 九九在线视频 | 亚洲精品久久久久久中文字幕小说 | 免费电影av | 99久久精品国产亚洲 | 久久亚洲国产精品日日av夜夜 | 久久香蕉国产线看观看网站 | 亚洲精品一区二区网址 | 午夜日韩精品 | 男女视频在线 | 奇米影视555 | 亚洲第一男人天堂 | 亚洲一区影院 | 97久久精品午夜一区二区 | 精品久久久久久久人人人人传媒 | 久久精品人人做人人 |