[Bioperl-l] Fwd: blast.pm patch
Robson Francisco de Souza
robfsouza at gmail.com
Fri Dec 18 18:07:34 UTC 2009
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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: blast_patched.pm
Type: application/octet-stream
Size: 91820 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20091218/3771d91c/attachment-0004.obj>
More information about the Bioperl-l
mailing list