[BiO BB] BLAST problems/Questions
Raj
Rajkumar.Bondugula at gmail.com
Mon Sep 27 15:24:56 EDT 2004
I am trying to search a local database with the profile of a sequence
rather than a sequence itself. For this, I first generate profile for
the sequence using 'nr' database and in the second step use the
generated profile to search the local database. I use the 24/7/2004,
windows version BLASTPGP program from NCBI website.
I was running the following two commands in the same order:
1)'blastpgp -d nr -i seq.txt -j 3 -C profile.bin -o seq.out -a 2'
2)'blastpgp -d sel25culled -i seq.txt -R profile.bin -o seq0.txt -a 2'
The first run is fine, but the second run produces the following lines:
[NULL_Caption] WARNING: posReadCheckpoint: Attempting to recover data
from previous checkpoint
[NULL_Caption] WARNING: posReadCheckpoint: Data recovered successfully
Why am I getting these messages, will these messages effect my
results in anyway?
Also, when I have ungapped alignments while search, I still see
gaps inserted into the protein sequences in the output when I run the
following command:
'blastpgp -d sel25culled -i seq.txt -R profile.bin -o seq0.txt -a 2 -g F'
Why is this happening? is there anyway to have just ungapped alignment
while searching?
Regards,
Raj
===================================
Rajkumar Bondugula
Digital Biology Laboratory
316 Engineering Bldg. North
University of Missouri-Columbia
===================================
More information about the BBB
mailing list