Viewing BLAST Results
In the previous section, we saved result file results/blastp.out
to an S3
bucket at destination s3://my-bucket/my-dir/results/blastp.out
.
One method to download the results file to your workstation and view it is to use the AWS CLI to interface with the S3 bucket where your results are stored. If you do not have the AWS CLI installed, please see the AWS CLI installation instructions for more information.
Using the AWS CLI command aws s3 cp
, the result file can be downloaded to your workstation. The
command below copies the result file at S3 URI s3://my-bucket/my-dir/results/blastp.out
to your
working directory.
$ aws s3 cp s3://my-bucket/my-dir/results/blastp.out .
download: s3://my-bucket/my-dir/results/blastp.out to ./blastp.out
To view the results, we can simply use the cat
command to output the contents of the file to our
terminal.
$ cat blastp.out
BLASTP 2.12.0+
Reference: Stephen F. Altschul, 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 composition-based statistics: Alejandro A. Schaffer,
L. Aravind, Thomas L. Madden, Sergei Shavirin, John L. Spouge, Yuri
I. Wolf, Eugene V. Koonin, and Stephen F. Altschul (2001),
"Improving the accuracy of PSI-BLAST protein database searches with
composition-based statistics and other refinements", Nucleic Acids
Res. 29:2994-3005.
Database: Nurse shark proteins
7 sequences; 922 total letters
Query= sp|P01349.2|RELX_CARTA RecName: Full=Relaxin; Contains: RecName:
Full=Relaxin B chain; Contains: RecName: Full=Relaxin A chain
Length=44
Score E
Sequences producing significant alignments: (Bits) Value
P80049.1 RecName: Full=Fatty acid-binding protein, liver; AltName... 14.2 0.96
>P80049.1 RecName: Full=Fatty acid-binding protein, liver; AltName: Full=Liver-type
fatty acid-binding protein; Short=L-FABP
Length=132
Score = 14.2 bits (25), Expect = 0.96, Method: Compositional matrix adjust.
Identities = 3/9 (33%), Positives = 6/9 (67%), Gaps = 0/9 (0%)
Query 2 LCGRGFIRA 10
+C R ++R
Sbjct 123 VCTREYVRE 131
Lambda K H a alpha
0.334 0.143 0.520 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: 22680
Database: Nurse shark proteins
Posted date: Jun 18, 2024 4:42 PM
Number of letters in database: 922
Number of sequences in database: 7
Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Neighboring words threshold: 11
Window for multiple hits: 40