Question: samtools mpileup's options -b (more than one bam file) and -r cannot be used together
0
gravatar for jing.mengrabbit
5 months ago by
jing.mengrabbit10 wrote:

Hi, I would like to generate a mpileup file from 2 bam files in some regions instead of the whole chromosomes:

samtools -b bam_list_file -r regions_to_be_parsed.bed > mpileup_from_two_bam_files.pileup

The content in the bam_list_file is

~/1.bam

~/2.bam

but an error was given:

[E::mpileup] fail to parse region 'regions_to_be_parsed.bed' with 1.bam.

I do not think there is something wrong with the bam files or the bed file, because it works when I ran with the following command:

samtools -r regions_to_be_parsed.bed ~/1.bam > mpile1.pileup

samtools -r regions_to_be_parsed.bed ~/2.bam > mpile2.pileup

samtools -b bam_list_file > mpileup_from_two_bam_files.pileup

Is there some way to generate a mpileup file from two bam files with regions in bed file? Thanks for your time!

sequence software error • 237 views
ADD COMMENTlink modified 5 months ago by h.mon20k • written 5 months ago by jing.mengrabbit10

Hello jing.mengrabbit,

Please use the formatting bar (especially the code option) to present your post better.

Formatting bar

For example:

samtools -b bam_list_file -r regions_to_be_parsed.bed > mpileup_from_two_bam_files.pileup

Besides, I believe your command is:

samtools mpileup -b bam_list_file -r regions_to_be_parsed.bed > mpileup_from_two_bam_files.pileup
ADD REPLYlink written 5 months ago by h.mon20k

Thanks for your advice.

ADD REPLYlink written 5 months ago by jing.mengrabbit10
2
gravatar for h.mon
5 months ago by
h.mon20k
Brazil
h.mon20k wrote:

Are you sure about -r regions_to_be_parsed.bed? The help says:

-r, --region REG region in which pileup is generated

I think -r should be used as -r chr1:200-210. For using bed files as intervals, you should use -l:

-l, --positions FILE skip unlisted positions (chr pos) or regions (BED)

ADD COMMENTlink modified 5 months ago • written 5 months ago by h.mon20k

Thanks for your reply. I misunderstood the helper documentation. It should be -l instead of -r as you said.

ADD REPLYlink written 5 months ago by jing.mengrabbit10
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: 1892 users visited in the last hour