DECONSEQ error: cannot find all database files
0
1
Entering edit mode
8.8 years ago

Dear all, I created a database of human chromosomes using the ftp.ncbi.nlm.nih.gov/genomes repository. The database was created in the folder refChr and contained the files:

hs_ref_GRCh38_p2.amb  hs_ref_GRCh38_p2.pac   hs_ref_GRCh38_p2.sa
hs_ref_GRCh38_p2.ann  hs_ref_GRCh38_p2.rbwt  hs_ref_GRCh38_p2_split.fa
hs_ref_GRCh38_p2.bwt  hs_ref_GRCh38_p2.rpac  hs_ref_GRCh38_p2_split.fa.log
hs_ref_GRCh38_p2.fa   hs_ref_GRCh38_p2.rsa   hs_ref_GRCh38_p2_split_PS.fa.fasta

I then modified the DECONSEQ config file DeconSeqConfig.pm in order to point to the right database with the lines:

use constant DBS => {hs_ref_GRCh38_p2 => {name => 'Human Reference GRCh38',  
db => 'hs_ref_GRCh38_p2'}, ...
use constant DB_DEFAULT => 'refChr';

and then I ran DECONSEQ as follows:

perl /usr/local/lib/Deconseq/deconseq.pl -f fu_1.fq -dbs ./refChr/hs_ref_GRCh38_p2 -i 90 -c 90 -out_dir DECONSEQ

where fu_1.fq is the file to be analyzed but I obtained:

ERROR: database "./refChr/hs_ref_GRCh38_p2" does not exist in config file.

When running:

perl /usr/local/lib/Deconseq/deconseq.pl -f fu_1.fq -dbs hs_ref_GRCh38_p2 -i 90 -c 90 -out_dir DECONSEQ

the result was:

ERROR: cannot find all database files for database "hs_ref_GRCh38_p2" in dir "db/".

Could you please tell what I am getting wrong?
Thanks
Luigi

sequence • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 1935 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6