[EMBOSS] prophet BLOSUM question
Daniel Barker
db60 at st-andrews.ac.uk
Thu Apr 10 07:10:11 EDT 2008
Dear EMBOSS users,
I am new to protein profile searches and am using prophecy and prophet
for a distant similarity search. I'm trying to use BLOSUM45, but prophet
is reporting BLOSUM62.
I'm generating profiles using prophecy and the BLOSUM45 matrix, thus:
prophecy -auto -warning -type H -datafile EBLOSUM45 -name
homologene_1825.aln -outfile homologene_1825.prophecy
This seems fine.
I'm using the profiles to search a protein database, thus:
prophet -auto -warning -infile homologene_1825.prophecy -outfile
prophet_out.txt -sequence
NC_000909_Methanocaldococcus_jannaschii_DSM_2661_prots.fasta
I am puzzled by this line I get in my prophet output (prophet_out.txt in
this case):
# Matrix: EBLOSUM62
I can't seem to make it go away. The profile file itself
(homologene_1825.prophecy) does correctly report BLOSUM45:
Matrix EBLOSUM45
Is there any way to make prophet also use BLOSUM45? (Does it in fact
make any difference to prophet?)
Thank you very much for advice.
Best regards,
Daniel
--
Daniel Barker
http://bio.st-andrews.ac.uk/staff/db60.htm
The University of St Andrews is a charity registered in Scotland :
No SC013532
More information about the EMBOSS
mailing list