使用 biopython 解析 blastp 输出答案

【问题标题】：parsing blastp output using biopython使用 biopython 解析 blastp 输出
【发布时间】：2014-10-05 22:15:05
【问题描述】：

我需要帮助来解析以下 blastp 输出以对齐 2 个序列：

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

【问题讨论】：

你尝试了什么？你有什么问题？你看过 biopython 教程吗？

标签： bioinformatics biopython blast

【解决方案1】：

您必须阅读教程中的此条目：http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec100

您将学习如何在本地启动 BLAST 并获得 XML 输出。一旦你有了你的 XML 文件（而不是你粘贴在你的问题中的文本），你可以解析它：

NCBIXML.parse(xml_results)

【讨论】：