Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Write the subset of reads from BAM file into new SAM/BAM file, using R tools. Old Pioneer Bioinformatics 0 01-27-2016 05:41 AM
intersectBed (BEDtools) generating empty output file palc Bioinformatics 1 08-28-2012 09:36 AM
How to exclude the unmapped pair reads in bam files vincent lee Bioinformatics 2 08-10-2012 03:18 AM
Exclude certain tags from BAM file yksikaksi Bioinformatics 5 05-06-2012 11:49 PM
intersectBed gives empty file Boel Bioinformatics 2 01-04-2011 07:29 AM

Thread Tools
Old 05-09-2018, 09:32 AM   #1
Location: Boston

Join Date: Mar 2013
Posts: 11
Smile question about using “intersectBed –v –s” to exclude some reads in the bam file


When I used “intersectBed -abam aligned.bam -b genes.bed -v -s -ubam > aligned.filtered.bam” to filter the reads in the aligned,bam that overlap with the entries in genes.bed, the output file “aligned.filtered.bam” has even much bigger size than the original bam file “aligned.bam”. That is not consistent the functions of removing regions in the bam file with “-v” option.

Has anyone met the similar issue when using "intersectBed -v" bedtools?

Any suggestion?

tinkering is offline   Reply With Quote
Old 05-10-2018, 05:18 AM   #2
Senior Member
Location: USA, Midwest

Join Date: May 2008
Posts: 1,155

The option "-ubam" directs intersectBed to write the output as uncompressed BAM. Since your input was undoubtedly compressed it's not surprising that the uncompressed output is larger even if it actually contains less data.
kmcarr is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:48 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO