Question: What are recommended parameters for removing rRNA from meta[transcript][gen]omics using
gravatar for O.rka
8 weeks ago by
O.rka170 wrote:

I've seen a lot of references to removing rRNA with but I haven't seen specific commands that should be used. Is there a recommended parameter setting that could be used? Would this change with respect to read length?

My plan was to do the following:

  1. kneaddata to run trimmomatic and remove human contamination
  2. to bin out rRNA (might do something w/ this later)
  3. to assemble the transcripts

I'm very familiar with steps 1 and 3 but I haven't experimented much with I'm very impressed with the rest of the bbsuite so I'm looking forward to adding this into the pipeline.

rna-seq metagenomics • 90 views
ADD COMMENTlink written 8 weeks ago by O.rka170

perhaps this would help:

Alternatively, you could use sortmerna

ADD REPLYlink written 8 weeks ago by andres.firrincieli590

@O.rka: Take a look at C: How to identify 16s sequences from binning data(contigs)? in addition (specifically 2nd comment from @Brian in answer there).

If you are working with a eukaryotic sample then you may want to use NCBI's ITS, LSU and SSU rRNA sequences they make available in pre-made blast indexes directory. You will need to extract fasta files and make k-mers from them.

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by genomax83k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1301 users visited in the last hour