Question: Viewing Multiple Regions In Samtools
1
gravatar for axelwilhelm
7.3 years ago by
axelwilhelm100
Sweden
axelwilhelm100 wrote:

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 • 3.4k views
ADD COMMENTlink modified 7.3 years ago by Pierre Lindenbaum119k • written 7.3 years ago by axelwilhelm100
6
gravatar for Pierre Lindenbaum
7.3 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum119k wrote:

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]
ADD COMMENTlink written 7.3 years ago by Pierre Lindenbaum119k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1786 users visited in the last hour