-5

2 つのシーケンスを整列させるために、次の 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
4

1 に答える 1

0

チュートリアルでこのエントリを読む必要があります: http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec100

BLAST をローカルで起動して XML 出力を取得する方法を学習します。XML ファイル (質問に貼り付けたテキストではない) を取得したら、次のように解析できます。

NCBIXML.parse(xml_results)
于 2014-10-06T12:16:47.530 に答える