Yesterday I ran some tests to compare bowtie2 and bbTools decontamination. I used the same index provided in your link for both. bbTools takes longer to compute (15-30% longer) and deletes fewer reads. It also seems like bbTools has lower sensitivity, even in heavily contaminated runs. I think I will stick with my previous workflow.
Host Decontamination: bbTools or Samtools?
Hi,
I have always used Samtools to remove any human contaminants from my gut flora data. However, bbTools also allows you to do that. All the post I've seen mention using bbTools for small contaminants removal. Is it practical to use it for removing hg19 contaminants or should I stick with Bowtie + Samtools?
1 answer
Yes it is totally practical. @Brian Bushnell wrote a pipeline (which uses BBTools) to do that. You can find the post here. You could also use bbsplit.sh from BBMap tools and bin your reads in multiple ways.
Log in to answer this question.
More posts like this
-
Query Regarding Outlier Removal in RRBS Data Analysis Using MethylKit
written by Srinka •
-
Input and output of picard markduplicate for deepvariant
written by hello_bioinf •
-
What are recommended parameters for removing rRNA from meta[transcript][gen]omics using bbduk.sh?
written by O.rka
-
Is It a Necessity to Remove Human Data from Samples for Metagenomic Analysis?
written by gtasource
-
Reheadering a subsetted bam gives truncated file
written by Luciferine •
-
Find human contaminants in metagenomics microbiome studies
written by David
-
Questions : Data management of whole exome sequencing
written by mangfu100
-
What Is The Best Tool For Mouse (Mm9 Or Mm10) Variant Annotations?
written by Charles Warden
-
Finding All Possible Exact Matches With Bowtie
written by jensenkristian89 •
-
Bioconductor/R For Ngs Analysis
written by Travis
Well, samtools has multiple subcommands, and you can do a lot with it. It is probably most helpful if you describe more in detail what you did (which commands you used).
Basic samtools stuff: align against hg19 with bowtie2, extract reads with both ends unmatched, sort and convert filtered bams to fastq.
I was just wondering if I could remove one docker container from my pipeline by using bbTools for both quality trimming and decontamination
Yes you can do that.