Merge pull request #231 from bioperl/pir-no-desc-line
[bioperl-live.git] / t / data / bug1986.blastp
blobd0400e39d017555a5d3bfbf55c378462c434ae33
1 BLASTP 2.2.11 [Jun-05-2005]
4 Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, 
5 Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
6 "Gapped BLAST and PSI-BLAST: a new generation of protein database search
7 programs",  Nucleic Acids Res. 25:3389-3402.
9 Query= NP_006065
10          (442 letters)
12 Database: Homo_sapiens.NCBI36.apr.pep.fa 
13            48,851 sequences; 23,910,368 total letters
15 Searching..................................................done
17                                                                  Score    E
18 Sequences producing significant alignments:                      (bits) Value
20 ENSP00000350182 pep:novel clone::BX322644.8:4905:15090:-1 gene:E...   120   3e-27
21 ENSP00000350182 pep:novel clone::BX322644.8:4905:15090:-1 gene:E...   120   3e-27
22 ENSP00000327738 pep:known-ccds chromosome:NCBI36:4:189297592:189...   115   8e-26
24 >ENSP00000350182 pep:novel clone::BX322644.8:4905:15090:-1 gene:ENSG00000137397
25            transcript:ENST00000357569 
26           Length = 425
28  Score =  120 bits (301), Expect = 3e-27
29  Identities = 76/261 (29%), Positives = 140/261 (53%), Gaps = 21/261 (8%)
31 Query: 9   IEKEVTCPICLELLTEPLSLDCGHSFCQACITAKIKESVIISRGESSCPVCQTRFQPGNL 68
32            +++EV CPICL++L +P+++DCGH+FC  CIT +I E+   S G   CP+C+T  +   +
33 Sbjct: 10  LQEEVICPICLDILQKPVTIDCGHNFCLKCIT-QIGET---SCGFFKCPLCKTSVRKNAI 65
35 Query: 69  RPNRHLANIVERVKEVKMSP-QEGQKRDVCEHHGKKLQIFCKEDGKVICWVCELSQEHQG 127
36            R N  L N+VE+++ ++ S  Q  +K   C  H +    FC++DGK +C+VC  S++H+ 
37 Sbjct: 66  RFNSLLRNLVEKIQALQASEVQSKRKEATCPRHQEMFHYFCEDDGKFLCFVCCESKDHKS 125
39 Query: 128 HQTFRINEVVKECQEKLQVALQRLIKEDQEAEKLED------DIRQERTAWKIERQKILK 181
40            H    I E  +  Q ++Q  +Q L ++++E  +++       D+  ++   + E+Q+IL 
41 Sbjct: 126 HNVSLIEEAAQNYQGQIQEQIQVLQQKEKETVQVKAQGVHRVDVFTDQV--EHEKQRILT 183
43 Query: 182 GFNEMRVILDNEEQRELQKL----EEGEVNVLDNLAAATDQLVQQRQDASTLISDLQRRL 237
44             F  +  +L+ E+   L ++     EG       +A+   QL     D   L+  L+ + 
45 Sbjct: 184 EFELLHQVLEEEKNFLLSRIYWLGHEGTEAGKHYVASTEPQL----NDLKKLVDSLKTKQ 239
47 Query: 238 TGSSVEMLQDVIDVMKRSESW 258
48                 ++L+     + RSE +
49 Sbjct: 240 NMPPRQLLEVTQPHLPRSEEF 260
52 >ENSP00000327738 pep:known-ccds chromosome:NCBI36:4:189297592:189305643:1
53            gene:ENSG00000184108 transcript:ENST00000332517
54            CCDS3851.1
55           Length = 468
57  Score =  115 bits (289), Expect = 8e-26
58  Identities = 101/410 (24%), Positives = 180/410 (43%), Gaps = 39/410 (9%)
60 Query: 8   DIEKEVTCPICLELLTEPLSLDCGHSFCQACITAKIKESVIISRGESSCPVCQTRFQPGN 67
61            ++ +E+TC ICL+  + P++ +CGHSFC  C+    +E         SCP C    +  +
62 Sbjct: 9   NLREELTCFICLDYFSSPVTTECGHSFCLVCLLRSWEE----HNTPLSCPECWRTLEGPH 64
64 Query: 68  LRPNRHLANIVERVKEVKMSPQEGQKRDVCEHHGK-----KLQIFCKEDGKVICWVCELS 122
65             + N  L  +    ++++   Q  Q  D    +G+     K     ++ G        ++
66 Sbjct: 65  FQSNERLGRLASIARQLR--SQVLQSEDEQGSYGRMPTTAKALSDDEQGGSAF-----VA 117
68 Query: 123 QEHQGHQTFRINEVVKECQEKLQVALQRLIKEDQEA------EKLEDDIRQERTAWKIER 176
69            Q H  ++    +E  +  +EKLQ  L  L    +EA      EK    + QE T  K  +
70 Sbjct: 118 QSHGANRVHLSSEAEEHHREKLQEILNLLRVRRKEAQAVLTHEKERVKLCQEET--KTCK 175
72 Query: 177 QKILKGFNEMRVILDNEEQRELQKLEEGEVNVLDNLAAATDQLVQQRQDASTLISDLQRR 236
73            Q ++  + +M   L  EEQ +LQ LE+ E   +  L     +L QQ +  S +I+ ++  
74 Sbjct: 176 QVVVSEYMKMHQFLKEEEQLQLQLLEQEEKENMRKLRNNEIKLTQQIRSLSKMIAQIESS 235
76 Query: 237 LTGSSVEMLQDVIDVMKRSESWTXXXXXXXXXXXXXXFRVPDLSGMLQVLKELTDVQYYW 296
77               S+ E L++V   ++RSE                   +  ++GM ++L++ +      
78 Sbjct: 236 SQSSAFESLEEVRGALERSE----PLLLQCPEATTTELSLCRITGMKEMLRKFS------ 285
80 Query: 297 VDVMLNPGSATSNVAISVDQRQVKTVRTCTFKNSNPCDF-SAFGVFGCQYFSSGKYYWEV 355
81             ++ L+P +A + + +S D + VK   +      NP  F  +  V G Q F+SG++YWEV
82 Sbjct: 286 TEITLDPATANAYLVLSEDLKSVKYGGSRQQLPDNPERFDQSATVLGTQIFTSGRHYWEV 345
84 Query: 356 DVSGKIAWILGVHSKISSLNKRKSSGFAFDPSVNYSKVYSRYRPQYGYWV 405
85            +V  K  W +G+     S    +       P   +S +  +    Y  WV
86 Sbjct: 346 EVGNKTEWEVGICKDSVS----RKGNLPKPPGDLFSLIGLKIGDDYSLWV 391
89   Database: Homo_sapiens.NCBI36.apr.pep.fa
90     Posted date:  Jun 15, 2006  8:56 PM
91   Number of letters in database: 23,910,368
92   Number of sequences in database:  48,851
93   
94 Lambda     K      H
95    0.319    0.133    0.398 
97 Gapped
98 Lambda     K      H
99    0.267   0.0410    0.140 
102 Matrix: BLOSUM62
103 Gap Penalties: Existence: 11, Extension: 1
104 Number of Hits to DB: 20,900,506
105 Number of Sequences: 48851
106 Number of extensions: 899179
107 Number of successful extensions: 6075
108 Number of sequences better than 1.0e-25: 105
109 Number of HSP's better than  0.0 without gapping: 18
110 Number of HSP's successfully gapped in prelim test: 87
111 Number of HSP's that attempted gapping in prelim test: 5632
112 Number of HSP's gapped (non-prelim): 157
113 length of query: 442
114 length of database: 23,910,368
115 effective HSP length: 107
116 effective length of query: 335
117 effective length of database: 18,683,311
118 effective search space: 6258909185
119 effective search space used: 6258909185
120 T: 11
121 A: 40
122 X1: 16 ( 7.4 bits)
123 X2: 38 (14.6 bits)
124 X3: 64 (24.7 bits)
125 S1: 41 (21.8 bits)
126 S2: 289 (115.9 bits)