This post is useful for removing specific chromosomes from a BAM: Remove mitochondrial reads from BAM files -Remove mitochondrial reads from BAM files
Removing Human Reads From Metagenome BAM File
So I have metagenome bam files which contain reads for both human and bacteria. I would like to remove all the human reads so I can focus on the bacterial reads. Has anyone been able to remove human reads from their RNA-seq data using samtools or any other software?
• 3,996 views
•
link
updated
by
WouterDeCoster
•
written
by
Kevin.Y •
•
written
by
Kevin.Y •
2 answers
Align your reads to the human genome first, then extract all the unmapped reads to work with further. The unmapped reads will have all your bacteria (as well as other non-human things). Extracting unmapped paired end from an alignment - C: Extracting unmapped paired end from an alignment
This sounds like a job for BBsplit
• 0 views
•
link
written
by
WouterDeCoster
Agree with @WouterDeCoster, BBSplit can be used, although it is very slow. Here is an example of it Tool to separate human and mouse rna seq reads from one of my posts.
You can probably use seal.sh from BBmap itself which is very fast. The example of Seal.sh is in the above link.
Log in to answer this question.
More posts like this
-
RNA-seq alignment in UNIX
written by Hien •
-
What is the MAPQ value of flag 512 for bam file?
written by hyjforesight •
-
RNA-seq of human cells and bacteria
written by bart •
-
Bacterial reads in ATAC-seq
written by C4 •
-
How to detect and remove bacterial sequences in assembled eukaryote genome
written by Seq225
-
Remove contamination before assembly
written by bruseq •
-
Using Taxonomizr to Call Bacterial Taxonomy From CSV File
written by Kevin.Y •
-
RNA-seq Bacterial Alignments Using NCBI Database
written by Kevin.Y •
-
Remove specific mapped reads from SAM file
written by Jennifer •
-
Using bowtie2 to create BAM files to view strand-specific coverage in Artemis
written by Chris S.
You have it in a bam file, meaning it's already aligned? To what?
It was aligned to a metagenome database for human and bacteria.