Hello,
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?
Thanks.
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?
Thanks.
Comment