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).

sortmerna index rrna • 553 views
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.

0
Entering edit mode

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

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.

0
Entering edit mode

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

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.