[Bioclusters] blastunsplit Error
zzar
bioclusters@bioinformatics.org
27 Nov 2002 18:14:53 +0900
hi
I'm using the seqsplit & blastunsplit
ftp://ftp.cgr.ki.se/pub/prog/MSPcrunch+Blixem/seqsplit.tar.Z
my work is
1. $seqsplit -c 12000 test.seq
2. $ls
test.0 test.1 test.2.test.3 test.4 test.5
3. blastwork
$ blastall -a2 -d nr -p blastx -i test.0(1,2,3,4,5) -o
test.0(1,2,3,4,5).bl
4. after I run the blastunslplit
$ blastunsplit test.seq
Reconstructing header ...
Chunk 0: 12465 residues [Total 12465]
Chunk 1: 12459 residues [Total 24924]
Chunk 2: 12456 residues [Total 37380]
Chunk 3: 12453 residues [Total 49833]
Chunk 4: 12445 residues [Total 62278]
Chunk 5: 1759 residues [Total 64037]
Total length of query with overlaps (from .bl files): 64037
FATAL ERROR: Wrong chunk and/or overlap sizes specified!
They give a total length of 254259, whereas the .bl files give 64037
--- or -----
$blastunsplit -c 12000 test.seq
Reconstructing header ...
Chunk 0: 12465 residues [Total 12465]
Chunk 1: 12459 residues [Total 24924]
Chunk 2: 12456 residues [Total 37380]
Chunk 3: 12453 residues [Total 49833]
Chunk 4: 12445 residues [Total 62278]
Chunk 5: 1759 residues [Total 64037]
Total length of query with overlaps (from .bl files): 64037
FATAL ERROR: Wrong chunk and/or overlap sizes specified!
They give a total length of 64259, whereas the .bl files give 64037
5. how can I do