This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Trinity : average transcript size check + quality assessment

Hello all

I have performed trinity assembly , now I want assess the quality of the the assembled transcript .. with the most closely related species of my plant

I have created database of annotated reference transcriptome

tuttu@kau:~/ngs/trinityrnaseq-2.1.1$ makeblastdb -in test.fasta -dbtype nucl

Building a new DB, current time: 12/09/2015 17:59:24
New DB name:   test.fasta
New DB title:  test.fasta
Sequence type: Nucleotide
Keep Linkouts: T
Keep MBits: T
Maximum file size: 1000000000B
Adding sequences from FASTA; added 61538 sequences in 4.21683 seconds.

After that I tried to run megablast to align the known transcripts to the Trinity assembly:

This is the command used :

blastn -query ~/ngs/trinityrnaseq-2.1.1/trinity_output_data2/Trinity.fasta -db /home/jinu/ngs/trinityrnaseq-2.1.1/test.fasta -out Trinity_vs_test.blastn -evalue 1e-20 -dust no -task megablast -num_threads 2 -max_target_seqs 1 -outfmt 6

But my attempt was unsuccessful

it results in this error

command line argument error: Argument "out". File is not accessible:  `test.blastn'

Thank you

ngs assembly megablast rna-seq

I already gone through that post ... but still its not working for me

For other hand I guess that the file which is not accessible is 'Trinity_vs_test.blastn', isn't it? Because 'test.blastn' is not specified in your command line..

Also the error could be related to permissions problem that prevents writing to a file in the given folder. Have you check it?

ls -lh command to see permissions, Then you can use chmod to change them. If you are not able to do it, google has a lot of examples related to permissions issues :).

0 answers

No answers yet.

Log in to answer this question.