sortmerna command not working
0
0
Entering edit mode
4.9 years ago
Shahzad ▴ 30

Dear All I need some help. I am getting an error in running sortmerna tool. I have successfully indexed the reference rRNA database. command used for indexing db

./indexdb_rna --ref ./rRNA_databases/silva-euk-28s-id98.fasta,./index/silva-euk-28s-db -v it made 4 files in the index folder named silva-euk-28s-db.bursttrie_0.dat
silva-euk-28s-db.kmer_0.dat
silva-euk-28s-db.pos_0.dat
silva-euk-28s-db.stats

Next step is to filter RNA-seq reads from rRNA indexed database command used for this was

./sortmerna --ref '/home/laptop/Desktop/NGS_tools/sortmerna-2.1/rRNA_databases','/home/laptop/Desktop/NGS_tools/sortmerna-2.1/index' --reads '/home/laptop/Desktop/RNA-seq/subdir1/output.fastq' --aligned sortmeout --fastx --other RNAsortme --log

this is resulting in following error

ERROR: The index '/home/laptop/Desktop/NGS_tools/sortmerna-2.1/index.stats' does not exist. Make sure you have constructed your index using the command indexdb'. Seeindexdb -h' for help.

Please help if you notice that I did something wrong here.

RNA-Seq sortmerna NGS Bioinformatics • 3.9k views
ADD COMMENT
1
Entering edit mode

I don't know the tool, but having glimpsed at the manual, should you not define a fasta file as a reference to the sortmerna tool (instead of just some folder like in your command)? Likewise, it looks to me like you're not defining the location of the index file correctly (also in error message). Shouldn't it be:

./sortmerna --ref /home/laptop/Desktop/NGS_tools/sortmerna-2.1/rRNA_databases/silva-euk-28s-id98.fasta,/home/laptop/Desktop/NGS_tools/sortmerna-2.1/index/silva-euk-28s-db --reads ..etc
ADD REPLY
0
Entering edit mode

Thank you @5heikki It solved my problem. Regards.

ADD REPLY

Login before adding your answer.

Traffic: 2155 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