Viewing Multiple Regions In Samtools

Hi all

I'm wondering if there is a way to exclude regions in samtools view.

Now I awk and sed lines away in post processing via piping.

Wouldn't it save time and space to just pileup more than one region at once into a file?

Maybe this is a dumb question but I don't have the software now while I'm wondering so I can't experiment.

For example:

samtools view chrY:2709520-3709520, chrY:17443438-27443438 ǀ samtools mpileup...
samtools view chrY:2709520-3709520 AND chrY:17443438-27443438 ǀ samtools mpileup.
samtools view chrY:2709520-3709520: chrY:17443438-27443438 ǀ samtools mpileup...

...

samtools

• 6,196 views • link

updated 5.3 years ago by Pierre Lindenbaum

• written 14.5 years ago by axelwilhelm

1 answer

samtools view takes an argument to only output the alignments overlapping the regions of a BED file:

-L FILE  output alignments overlapping the input BED FILE [null]

• 0 views • link

written 14.5 years ago by Pierre Lindenbaum

update 2021: use option -M to go faster with an indexed bam.

• 0 views • link

written 5.3 years ago by Pierre Lindenbaum

More posts like this

Why Do We Need To Set Maximal Read Depth In Samtools
written 12.7 years ago by a.a.kechin •

I analyse reads from Ion Torrent platform with <a href='http://bio-bwa.sourceforge.net/'>bwa</a>-samtools pipeline. I used them in such way: <pre><code>bwa mem samtools view (convert sam to bam) …
Mpileup Problem
written 13.0 years ago by stolarek.ir

Hi. I have a problem with my variant calling using <a href='http://samtools.sourceforge.net/'>samtools</a>. I'm interested in vcf file with indels in my <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. …
Snp Calling With Bcftools Using Pileup File Format Instead Of Bcf
written 13.1 years ago by Justin

Usually when I do SNP calling, I do this: <pre><code>samtools mpileup -g -f Ref.fa File.bam | bcftools view -cv - </code></pre> where <a href='http://samtools.sourceforge.net/'>samtools</a> outputs …
Will Samtools Mpileup Produce Different Variant Calls On A Subset Of The Same Data?
written 13.5 years ago by mark.dunning

Hi al, I'm looking to parallelise my variant-calling using mpileup on a per-chromosome basis as to run in serial for a whole genome takes many, …
Varscan Error On Mpileup Input
written 13.7 years ago by GPR

Hello, I am trying to filter low quality data off my mpileup files using Varscan but keep getting stuck with the error message: Invalid format …
Samtools: Why Have I Been Using "View" Command For Mpileup
written 13.8 years ago by axelwilhelm

This is a basic question so maybe someone can make the best of it and explain some basic details so that newbies can benefit from …
How To Switch From Samtools Pileup To Mpileup
written 14.3 years ago by groundcover •

Specifically I want to run this command: <pre><code>samtools pileup -f parent1.fa.msg update_reads-aligned-parent1.bam.sorted.bam -c > update_reads-aligned-parent1.pileup </code></pre> But everything I'm reading, I can't find the -c …
Why Does Mpileup Skip My Mutation ?
written 15.1 years ago by Pierre Lindenbaum

(cross posted on seqanswers: <a href='http://seqanswers.com/forums/showthread.php?t=12542'>http://seqanswers.com/forums/showthread.php?t=12542</a> ) Hi all, I know, by capillary sequencing, that one of my samples contains a mutation at position:<code>120458492</code>. Some …
Samtools Header From A Bacterial Genome Ba File
written 15.1 years ago by Sakti

Hello dear Bio masters, I'm analyzing some reads from a bacterial genome (E. coli), and I'm trying to get the pileup to calculate coverage for …
Why Import Command In Samtools Is Not Shown In The Short Help Of Samtools?
written 15.1 years ago by Pablo Marin-Garcia

I am converting sam to bam using 'samtools import' command. I am curious why if I do 'samtools' I can not see the 'import' command? …