[Bioperl-l] Fwd: blast.pm patch

Chris Fields cjfields at illinois.edu
Fri Dec 18 18:33:44 UTC 2009


Robson, 

Any chance you could check this against SVN?  We haven't used the CVS tree for a few years (had a number of releases along the way as well).

Not sure about bioperl-bugs, we have bugzilla still running though:

http://bugzilla.open-bio.org/

chris


On Dec 18, 2009, at 12:07 PM, Robson Francisco de Souza wrote:

> Hi,
> 
> I've been dealing with an apparent bug in the output of NCBI's BLAST
> programs (blastall, blastpgp) which sometimes produces output like the
> one below.
> I think I've managed to produce a work around for Bioperl blast.pm
> parser and would like to contribute it to Bioperl.
> The fix is based on blast.pm from the CVS tree (downloaded some months
> ago...) and is attached to this message.
> Best,
> Robson
> 
> PS: what happened to the bioperl-bugs mailing list? It does not seem
> to be working...
> 
>> gi|156552846|ref|XP_001600053.1| PREDICTED: similar to conserved
>           hypothetical protein [Nasonia vitripennis]
>          Length = 1774
> 
>  Score = 75.9 bits (185), Expect = 1e-11,   Method: Compositional matrix adjust.
>  Identities = 85/393 (21%), Positives = 175/393 (44%), Gaps = 28/393 (7%)
> 
> Query: 0   -
> 
> Sbjct: 328 P                                                            328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 0
> 
> Sbjct: 328                                                              328
> 
> Query: 612 VPPPPGSGIPMPPGGGFFGMKTKLP-----KLPELKATKDTKKIHIAG             654
>            P PP +   + P       KTK+      K+P  K         +
> Sbjct: 329 TPEPPNNSAKLLPQQEIPTPKTKMKTINWNKIPNHKVIGKRNIWSLVA             376
> 
> Query: 655 DKINNKDIEGTGWMSILEENAEKMSKIFDKN-LFENNFQKKETRDAPSQEKENVPTLVSF 713
>           ++  N  +    W  +     +++  +   N    NN       D   +E    PT ++
> Sbjct: 377 NEHQNSPMADLDWAEMEGLFCQQVPPMIPANTTCSNNLGNGVDTDKRRRE----PTEIAL 432
> 
> Query: 714 LDSKTSYQLALLLGFLKKNEREIRKHVIDLNEKELQKQTIHSLKDLCPEEDKFKEIESFV 773
>           LD K S  + + L   + +  +I + + D    ++  + +  L  + PE D+ + ++SF
> Sbjct: 433 LDGKRSLNVNIFLKQFRSSNEDIIQLIKDGGHDDIGAEKLRGLLKILPEVDELEMLKSF- 491
> 
> Query: 774 QKGDGYLEQLEPGDKLFYAMKDIPRLKQRFTAWSSQIYFEGSVISVEPDIESLNRACKNI 833
>              DG   +L   +K F  +  +P  K R      +  F  ++  +EP I S+  A +++
> Sbjct: 492 ---DGDKLKLGNAEKFFLQLIQVPNYKLRIECMLLKEEFAANMSYLEPSINSMILAGEDL 548
> 
> Query: 834 VQCKSLQRLMTLIVLLVNFLNKAKTDKDRVYGFKLNFLTKLGDIKSSSDPNRSMMNYLCE 893
>           +  KSLQ ++ ++++  NFLN      +   G KL+ L KL +I++    N+  MN L
> Sbjct: 549 MTNKSLQEVLYMVLVAGNFLNSGGYAGN-AAGVKLSSLQKLTEIRA----NKPGMN-LIH 602
> 
> Query: 894 FLLAKDDKLIPELLKELK--DYAEVGSRIELPELKKEIGKLNESLKVIQTELEFYKKEQK 951
>           ++  + ++   +LL   +  +  +  ++  + +L  E   L+  +K I+++++    E
> Sbjct: 603 YVAMQAERKRKDLLNFARGMNALDSATKTTVEQLTNEFNALDTRIKKIRSQIQLPTTEA- 661
> 
> Query: 952 FINDKFPKQLDEFYQYAKSEMQKINKAQEKLEKILKEVAKFFGE 995
>                  +Q+ +F Q A+ EM ++ +  E+L+ + + +A+FF E
> Sbjct: 662 ----DIQEQMAQFLQMAEQEMSQLKRDMEELDGVRRTLAEFFCE 701
> <blast_patched.pm>_______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l





More information about the Bioperl-l mailing list