1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
BLASTP 2.2.26+ [I/O Module of SeqAn-2.0.0, https://www.seqan.de]
Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer,
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997),
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs", Nucleic Acids Res. 25:3389-3402.
Reference for SeqAn: Doering, A., D. Weese, T. Rausch, K. Reinert (2008): SeqAn --
An efficient, generic C++ library for sequence analysis. BMC Bioinformatics,
9(1), 11. BioMed Central Ltd. doi:10.1186/1471-2105-9-11
Database: The Foo Database
2 sequences; 400 total letters
Query= Query_Numero_Uno with args
Length=12
Score E
Sequences producing significant alignments: (Bits) Value
Subject_Numero_Uno 23 0.0005
Subject_Numero_Dos 22 0.0009
ALIGNMENTS
> Subject_Numero_Uno
Length=200
Score = 23.1 bits (48), Expect = 0.0005
Identities = 10/14 (71%), Positives = 12/14 (86%), Gaps = 2/14 (14%)
Query 1 VAYAQ--PRKLCYP 12
VAYAQ PR+LC+P
Sbjct 79 VAYAQTKPRRLCFP 92
> Subject_Numero_Dos
Length=200
Score = 22.3 bits (46), Expect = 0.0009
Identities = 8/8 (100%), Positives = 8/8 (100%), Gaps = 0/8 (0%)
Query 3 YAQPRKLC 10
YAQPRKLC
Sbjct 157 YAQPRKLC 164
Lambda K H
0.267 0.0410 0.1400
Gapped
Lambda K H
0.267 0.0410 0.1400
Effective search space used: 4800
Query= Query_Numero_Dos with args
Length=10
***** No hits found *****
Lambda K H
0.267 0.0410 0.1400
Gapped
Lambda K H
0.267 0.0410 0.1400
Effective search space used: 4000
Query= Query_Numero_Tres with args
Length=8
Score E
Sequences producing significant alignments: (Bits) Value
Subject_Numero_Dos 18 0.007
Subject_Numero_Uno 16 0.03
ALIGNMENTS
> Subject_Numero_Dos
Length=200
Score = 18.9 bits (37), Expect = 0.007
Identities = 8/8 (100%), Positives = 8/8 (100%), Gaps = 0/8 (0%)
Query 1 AVITSFTQ 8
AVITSFTQ
Sbjct 10 AVITSFTQ 17
> Subject_Numero_Uno
Length=200
Score = 16.9 bits (32), Expect = 0.03
Identities = 7/8 (88%), Positives = 8/8 (100%), Gaps = 0/8 (0%)
Query 1 AVITSFTQ 8
AVITS+TQ
Sbjct 184 AVITSWTQ 191
Lambda K H
0.267 0.0410 0.1400
Gapped
Lambda K H
0.267 0.0410 0.1400
Effective search space used: 3200
Database: The Foo Database
Number of letters in database: 400
Number of sequences in database: 2
Matrix:BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
|
