解析使用biopython

Question

我需要幫助解析以下BLASTP輸出BLASTP輸出對準兩個序列：

BLASTP 2.2.28+ 


Query= 
Length=237 

Subject= 
Length=268 


Score = 429 bits (1104), Expect = 2e-157, Method: Compositional matrix adjust. 
Identities = 237/268 (88%), Positives = 237/268 (88%), Gaps = 31/268 (12%) 

Query 1 MERYENLFAQLNDRREGAFVPFVTLGDPGIEQSLKIIDTLIDAGADALELGVP------- 53 
      MERYENLFAQLNDRREGAFVPFVTLGDPGIEQSLKIIDTLIDAGADALELGVP  
Sbjct 1 MERYENLFAQLNDRREGAFVPFVTLGDPGIEQSLKIIDTLIDAGADALELGVPFSDPLAD 60 

Query 54 --TIQNANLRAFAAGVTPAQCFEMLALIREKHPTIPIGLLMYANLVFNNGIDAFYARCEQ 111 
       TIQNANLRAFAAGVTPAQCFEMLALIREKHPTIPIGLLMYANLVFNNGIDAFYARCEQ 
Sbjct 61 GPTIQNANLRAFAAGVTPAQCFEMLALIREKHPTIPIGLLMYANLVFNNGIDAFYARCEQ 120 

Query 112 VGVDSVLVADVPVEESAPFRQAALRHNIAPIFICPPNADDDLLRQVASYGRGYTYL---- 167 
      VGVDSVLVADVPVEESAPFRQAALRHNIAPIFICPPNADDDLLRQVASYGRGYTYL  
Sbjct 121 VGVDSVLVADVPVEESAPFRQAALRHNIAPIFICPPNADDDLLRQVASYGRGYTYLLSRS 180 

Query 168 ---------------LIEKLKEYHAAPALQG-GISSPEQVSAAVRAGAAGAISGSAIVKI 211 
          LIEKLKEYHAAPALQG GISSPEQVSAAVRAGAAGAISGSAIVKI 
Sbjct 181 GVTGAENRGALPLHHLIEKLKEYHAAPALQGFGISSPEQVSAAVRAGAAGAISGSAIVKI 240 

Query 212 IEKNLASP--MLAELRSFVSAMKAASRA 237 
      IEKNLASP MLAELRSFVSAMKAASRA 
Sbjct 241 IEKNLASPKQMLAELRSFVSAMKAASRA 268 



Lambda  K  H  a   alpha 
    0.320 0.136 0.386 0.792  4.96 

Gapped 
Lambda  K  H  a   alpha sigma 
    0.267 0.0410 0.140  1.90  42.6  43.6 

Effective search space used: 51972 




Matrix: BLOSUM62 
Gap Penalties: Existence: 11, Extension: 1 
Neighboring words threshold: 11 
Window for multiple hits: 40 

Answer 1

你必須閱讀教程此項：http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec100

您將學習如何啓動本地BLAST並獲取XML輸出。一旦你有你的XML文件（而不是你粘貼在你的問題的文本），你可以解析它：

NCBIXML.parse(xml_results)