[Bioperl-l] blastall and blosum62 problem
Mariano Latorre A
malatorr at genoma.ciencias.uchile.cl
Tue Jun 10 16:39:43 EDT 2003
I have a standalone blastall and a html version of blast on my server
and i am having problems with my results:
When I try the HTML version, I submit this pair of sequences
("sequences.in"):
>Test1
AGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGC
TTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAA
TATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACC
ATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAG
>Test2
CCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAA
GTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCC
AGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTG
AAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTT
I try a TBLASTX with BLOSUM62 and I get good results.
But I need to use blastall with the script version...so I try with this
command:
cat sequences.in | blastall -p tblastx -d /usr/local/blast/data/db1 -M
BLOSUM62 | more
and I get this errors and warnings:
[blastall] WARNING: [000.000] Test1: Unable to open BLOSUM62
[blastall] WARNING: [000.000] Test1: BlastScoreBlkMatFill returned
non-zero status
[blastall] WARNING: [000.000] Test1: SetUpBlastSearch failed.
[blastall] WARNING: [000.000] Test2: Unable to open BLOSUM62
[blastall] WARNING: [000.000] Test2: BlastScoreBlkMatFill returned
non-zero status
[blastall] WARNING: [000.000] Test2: SetUpBlastSearch failed.
But remember this: when I try blastn I get good results again:
cat sequences.in | blastall -p blastn -d /usr/local/blast/data/db1 -M
BLOSUM62 | more
I need this working first, then I need to code some perl scripts.
thanks,
Mariano
More information about the Bioperl-l
mailing list