This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filtering for redundancy in transcriptome de novo assembly using CD-HIT EST_fatal error: file opening failed Program halted

Hi,

I'm performing transcriptome de novo assembly from plant RNA-seq data using two strategies (the multiple k-mer approach and the single k-mer approach uising different assemblers) and then merging the several pre-assemblies. Now I'm trying to filter for redundancy using CD-HIT EST (Program: CD-HIT, V4.6), but the program stopped with an error.

Error looks like this:

228085  finished     175242  clusters

Apprixmated maximum memory consumption: 455M
writing new database

Fatal Error:
file opening failed
Program halted !!

Can someone help me to understand what happens?

Thanks

rna-seq assembly
• 5,115 views • link
updated by Ram Ram • written by c.evang • c.evang

Is this the first write operation? Could be a permission issue..

• 0 views • link
written by 5heikki 5heikki

Could you write the command you have used, please?
And as @5heikki has pointed out.. maybe is a permission issue. Check out folder permissions, files permissions...

• 0 views • link
written by iraun iraun

This is the command:

cd-hit-est -i /media/universita/Archivio/WATBIO/CHIARA/ARUNDO_RNA-Seq_092015/TransAbySS_03-04122015/transabyss_merge/MergeAssembly_K21_25/MergeAssembly_k21_25 \
           -o /media/universita/Archivio/WATBIO/CHIARA/ARUNDO_RNA-Seq_092015/cd-hit-est_output \
           -c 0.95 -n 8 -M 0 -T 8
• 0 views • link
updated by Ram Ram • written by c.evang • c.evang

Please edit your question instead of posting an 'answer', this will increase the readability of the post for future readers.

For other hand, I guess that this /media/universita/Archivio/WATBIO/CHIARA/ARUNDO_RNA-Seq_092015/TransAbySS_03-04122015/transabyss_merge/MergeAssembly_K21_25/MergeAssembly_k21_25 file is a fasta file, isn't it? Because you did not specify the extension and it is always a good practice.

• 0 views • link
updated by Ram Ram • written by iraun iraun

Sorry!

It is a fasta file; I tried adding the extension but I obtained always the same error.

• 0 views • link
updated by Ram Ram • written by c.evang • c.evang

Now it works properly; I have to add also the file extension for the output file.

Thanks for all your suggestions!

• 0 views • link
written by c.evang • c.evang

0 answers

No answers yet.

Log in to answer this question.

More posts like this