Run BLAST
Align the file p450s.fasta to the database we created in previous steps.
When doing the search, use an E-value cutoff of 1e-6, keep the top one target sequences, and produce an output file called p450s_blastp.blast in output format 11 (-outfmt 11).
p450s_blastp.blast existsp450s_blastp.blast contains the output of aligning sequences from p450s.fasta to the orf_trans database.