BASH SortMeRNA--Fasta Must be followed by an index name
0
0
Entering edit mode
3.4 years ago

Help! I'm trying to submit a script to sort out rRNA contamination from my RNA-sequencing using SortMeRNA. I tried running the script but just got this error message back:

ESC[0;31mERRORESC[0m: the FASTA reference file name /uoa/home/r01dj21/Software/sortmerna/sortmerna-2.1b/rRNA_databases/silva-euk-18s-id95.fasta must be followed by an index name.

What is an index name and how do I enter that into my script? It wasn't in the original guidelines for SortMeRNA.

Thank you!

index bash rRNA • 1.5k views
ADD COMMENT
1
Entering edit mode

From manual (copy/pasted in parts):

For Indexing:

The executable indexdb rna indexes an rRNA database.

Example 1: indexdb rna using one database
>> ./indexdb_rna --ref ./rRNA_databases/silva-bac-16s-id90.fasta,./index/silva-bac-16s-db -v

usage: ./indexdb_rna --ref db.fasta,db.idx [OPTIONS]:
--------------------------------------------------------------------------------------------------------
| parameter          value          description             default |
--------------------------------------------------------------------------------------------------------
--ref              STRING,STRING      FASTA  reference file, index file        mandatory
(ex. --ref /path/to/file1.fasta,/path/to/index1)
If passing multiple reference sequence files, separate
them by ’:’,
(ex. --ref /path/to/file1.fasta,/path/to/index1:/path/to/file2.fasta,path/to/index2)

For sortmeRNA

usage: ./sortmerna --ref db.fasta,db.idx --reads file.fa --aligned base_name_output [OPTIONS]:
-------------------------------------------------------------------------------------------------------------
| parameter value description default |
-------------------------------------------------------------------------------------------------------------
--ref STRING,STRING FASTA reference file, index file mandatory
(ex. --ref /path/to/file1.fasta,/path/to/index1)
If passing multiple reference files, separate
them using the delimiter ’:’,
(ex. --ref /path/to/file1.fasta,/path/to/index1:/path/to/file2.fasta,path/to
ADD REPLY

Login before adding your answer.

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