Help with mapping total RNA reads to human reference but want to keep unmapped reads (with bbmap)
I am interested in removing host (human) contamination from my total-RNA seq data to retain the unmapped reads for microbial analysis. I am using BBMAP which I have used previously for DNA seq data, but I haven't used it for RNA-seq before.
Looking at the guides it does have a suggestion for doing this, and I have run the program as follows:
bbmap.sh in1=forward_reads.fq in2=reverse_reads.fq outm=mapped.fq outu=unmapped.fq maxindel=200k ambig=random intronlen=20 xstag=us
Is that correct considering what I am aiming to do? I am just interested in the unmapped.fq file in the end.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Do threads matter for host DNA removal with bbmap?
written by HPC123 •
-
Help understanding BBMAP output
written by m.radz •
-
Feature count is very low using htseq-count
written by salmon •
-
Can you assemble with merged paired end reads and unmatched reads as "single ended" reads?
written by O.rka
-
bbmap split paired-end reads back into separated fastq files?
written by c.e.chong •
-
RNA-seq : Quantification problem after mapped with BBMap
written by k.kathirvel93
-
finding variants in the low to no coverage areas of bacterial genomes
written by pdm37 •
-
BBMap: Number of reads out does not add up to number of reads in
written by Sinji
-
bbmap: Order of reference genomes seems to matter
written by mschmid
-
BBDuk trimming concern: large discarded, too aggressive (?)
written by umn_bist •
Yes that is correct.
You could have also used the following. This looks for reads that share a kmer with a reference.
Would I use that instead of the code I used?