samtools view extract from .bed list to .bam file

Hi -

I am attempting to extract certain regions from a large .bam file into a smaller subsetted .bam file using samtools view. I'm doing so in order to have a smaller file to down/upload for viewing in IGV.

I have a .bed file with the regions of interest, tab-delimited, no extra white spaces (my "chromosomes" are actually specific genes in this case, so they each start with 0 and are comparatively short). This is a snippet of the first 3 rows. Total number of regions/rows/genes is between 60-90.

52283_0_001313  0   895
52283_0_001c02  0   3102
52283_0_001309  0   2203

I have used samtools view to attempt to extract these regions of interest:

samtools view -h -b -L regions.bed input.bam > output.bam

However, the output.bam file only consists of the first region from the .bed file (e.g., "52283_0_001313" only). I have double-checked that the "chromosome" names match in my .bam file and .bed file using samtools view -H input.bam

I am using samtools version 1.10.

Any suggestions on why this is happening or how to remedy this issue? Thank you.

samtools bam view bed

• 2,982 views • link

updated 4.8 years ago by GenoMax

• written 4.8 years ago by matthews.alix •

check the output of samtools view -c input.bam "52283_0_001c02:1-3102"

• 0 views • link

written 4.8 years ago by Pierre Lindenbaum

Thank you for your reply! The output of the above was 3378, but GenoMax's solution fixed the issue.

• 0 views • link

written 4.8 years ago by matthews.alix •

1 answer

Can you check your BED file with cat -vet your.bed to make sure your line endings are unix-compatible? Did you make this file on a PC/mac and then transfer it to your server? If so doing dos2unix your.bed may fix the problem.

• 0 views • link

written 4.8 years ago by GenoMax

Thank you for your quick reply! That was indeed the problem. I made the file on a Mac and transfered it to the server. I used mac2unix myfile.bed and that fixed it. Thank you very much!

Just in case others will find use in seeing the output:

> cat -vet myfile.bed

52283_0_001313^I0^I895^M52283_0_001c02^I0^I3102^M52283_0_001309^I0^I2203^M

Then

> mac2unix myfile.bed 
> cat -vet myfile.bed

52283_0_001313^I0^I895$
52283_0_001c02^I0^I3102$
52283_0_001309^I0^I2203$

• 0 views • link

written 4.8 years ago by matthews.alix •

More posts like this

Filter out Bam not overlapping with Bed File (keeping a read and its mate)?
written 3.8 years ago by Eliveri •

I have a Bam file for which I want to REMOVE reads overlapping with a given Bed file (so keep all other reads/alignments). I also …
Export sequence from specific region from bam file.
written 3.9 years ago by malkomalktho •

I have a bam file and I wanted to extract the sequence of a specific region. for example ch1:1-337, the sequence starts EXACTLY at position …
Extract reads overlapping a specific region in bam file
written 5.2 years ago by User000

Hello, ref --------------start-----------------------------------------stop------------------------------- r1 ------------------------------------------------------- r2 --------------------------------------------------------------- r3 ----------------- -------------------------------------- r4 ----------------------------------------------------------------------------- a. I would like to extract reads from bam file that overlap entirely …
How to extract and cut sequence from a given region from a BAM file
written 6.5 years ago by nbargues •

Hi, I have an aligned BAM file and I want to extract sequence from a given region and just that, not what's come after or …
Replace @RG group by the one of another BAM
written 7.4 years ago by da.de •

Hi, I want to replace the @RG info from one bam file (output.bam) by the @RG info from another (input.bam) creating a new.bam I tried …
Does samtools view -b -h <input.bam> -L <input.bed> keep the original header?
written 8.3 years ago by rubic

Simple question Seems like the command: samtools view -b -h <input.bam> -L <input.bed> for subsetting the `bam` file to keep only those reads intersecting the …
Extracting reads from given region
written 8.7 years ago by KVC_bioinfo

Hello all, I am extracting the reads from bam file which are falling in particular region using the following command. samtools view -h input.bam "Chr10:18000-45500" …
Samtools gives "truncated file" when trying to view a "reheaded" file
written 8.9 years ago by olavur

I have a BAM file that works fine, but after using `samtools reheader`, I cannot view it in samtools. The new file has the same …
Extract reads where both pairs maps properly and none of the pairs are split from a bam file
written 9.0 years ago by rimjhim.roy.ch

Hi, I want to extract only those reads from bam file in a given region where both pairs are mapped properly and are not split. …
Truncating BAM Files
written 10.7 years ago by JameB11 •

Hello, I am trying to create a smaller BAM file from the original BAM file. The smaller BAM file should have info on only about …