This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Getting 'Samtools View' To Work On Multiple Bamfiles Simultaneously

Hello,

I want to extract regions of multiple bamfiles kept as seperate files using the following samtools command:

samtools view -b a.bam refsequence :113111720-113112765 > aoutput.bam

How can I modify the above command so it works on multiple bamfiles example:a.bam b.bam c.bam d.bam simultaneously? I am currently doing this one bamfile at a time and it is taking forever. Thanks

samtools bam
• 8,122 views • link
updated by lh3 lh3 • written by Kwame • Kwame

2 answers

To output one BAM:

samtools merge -r chr1:100-1000 - *.bam > out.bam

To output multiple BAMs:

ls *.bam|sed s,.bam,,|xargs -i echo samtools view -b {}.bam chr1:100-1000 > part-{}.bam|sh
• 0 views • link
written by lh3 lh3

Or far less complicated:

for f in *.bam ; do samtools view -b $f chr1:100-1000 > part-$f ; done

• 0 views • link
written by Marvin • Marvin

I don't think you can. If you made the .bam files with read group information, you could merge them togehter, then run the view command on that.

• 1 views • link
written by Swbarnes2 Swbarnes2

Log in to answer this question.

More posts like this