What is index in sortmerna?
0
0
Entering edit mode
7 months ago
taufiqbge • 0

Hi, I am a beginner of NGS data analysis. I am using SortMeRNA to remove rRNA from my mapping sequence. I read all the available instruction SortMeRNA but still confused. My command is-

 ./indexdb_rna --ref ./silva-bac-16s-id90.fasta, ./folder containing input fasta file --reads ./folder containing input fasta file/input.fastq

ERROR: -m [INT] must be a positive integer value (in Mbyte).

Any help would be much appreciate.

sortmerna RNA-Seq next-gen sequencing • 514 views
ADD COMMENT
0
Entering edit mode

Multiple things you need to pay attention to.

You can't use random , in your command lines. ./folder containing input fasta file need to be replaced with a real entry` (did you copy that from the manual example command line).

./indexdb_rna --ref ./rRNA_databases/silva-bac-16s-id90.fasta,./index/silva-bac-16s-db -v
ADD REPLY
0
Entering edit mode

Thanks for your reply. I checked the (,) and also went through the manual. I just mentioned "folder containing input fasta file" to ensure that the path is correct before and after the --reads. However, still I am getting error. I input fastq file which is already mapped.

ADD REPLY
1
Entering edit mode

Ah I see. You used folder containing input fasta file as an indicator that you have a real value there.

I posted an example command for the indexing step from sortmerna manual. This is done independently of running sortmerna. You don't provide any reads/alignments in this step.

After the indexing is complete you provide aligned reads when running the program.

./sortmerna --ref db.fasta,db.idx --reads file.fq --aligned base_name_output
ADD REPLY
0
Entering edit mode

Sorry for not understanding the trick. Just confused that should I make a folder name "index" in the same path. I am confused this portion only-

> *./index/silva-bac-16s-db -v*

Thank you for your assistance.

ADD REPLY

Login before adding your answer.

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