Sortmerna: index fna file from ncbi
1
0
Entering edit mode
7 months ago
Stella • 0

Hello,

after running the rrna removal on my rnaseq samples I saw that I have some overpresented sequences in one of them. By blasting, I saw that they are rrna but from subunits 45s,40s,48s,5.8s. As these are not included in the sortmerna database and my sequences were not compared against them, they were left over and that sample is a major outlier.

Is it possible to index fna files that contain the sequences for these other subunits downloaded from ncbi and rerun the removal including them? I haven't succeded with the normal function (.stats file is few bytes big).

Many thanks in advance!

sortmerna index rrna • 553 views
ADD COMMENT
0
Entering edit mode
7 months ago
GenoMax 107k

Take a look at indexdb_rna -h command help to see how you can do that.

ADD COMMENT
0
Entering edit mode

I don't see some other option on indexing... I have used the standard command of --ref .fasta,db.idx

ADD REPLY
0
Entering edit mode

You should be able to replace the input fasta with one you want to use. .fna is NCBI's fasta file designation.

ADD REPLY
0
Entering edit mode

Ok, thank you. I will try running it even if the file seems empty.

ADD REPLY
1
Entering edit mode

No that would not help. You should figure out why the index files are empty. Perhaps sortmeRNA requires the extension to be .fasta etc.

ADD REPLY

Login before adding your answer.

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