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

I've seen a lot of references to removing rRNA with bbduk.sh 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. bbduk.sh to bin out rRNA (might do something w/ this later)
  3. rnaspades.py to assemble the transcripts

I'm very familiar with steps 1 and 3 but I haven't experimented much with bbduk.sh. 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: http://seqanswers.com/forums/showthread.php?t=42776&page=17

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.

Help
Access

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