Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Samtools mpileup -g and -u flags throwing errors for a specific bam file
0
gravatar for top1214
2.8 years ago by
top1214 • 20
top1214 • 20 wrote:

I'm trying to run to get variants called from a (merged) bam file (presumably using bfctools). 

$ samtools mpileup merged.bamfile.bam

works fine, spitting out in pileup format.

However, 

$ samtools mpileup -g merged.bamfile.bam

I get the following error message

   [mpileup] 1 samples in 1 input files
   <mpileup> Set max per-file depth to 8000
   Assertion failed: (q), function group_smpl, file bam_plcmd.c, line 287.
   Abort trap: 6

The same thing happens with -u or -v flags.

This is not a general problem, as I can use the flags on a different bam file.
samtools • 1.1k views
ADD COMMENT • link •
modified 2.8 years ago by Pierre Lindenbaum ♦ 126k • written 2.8 years ago by top1214 • 20
0
gravatar for Pierre Lindenbaum
2.8 years ago by
Pierre Lindenbaum ♦ 126k
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum ♦ 126k wrote:

there is no such 'assertion' in bam_plcmd.c, line 287. https://github.com/samtools/samtools/blob/master/bam_plcmd.c#L287 , you should update your version of samtools.

otherwise, there is a problem with your read-groups. Did you add some Read-group in your bam ? Adding Read Groups To Bam Files
ADD COMMENT • link written 2.8 years ago by Pierre Lindenbaum ♦ 126k

I didn't realize there was an update 2 months ago. I was using 1.3.1

ADD REPLY • link written 2.8 years ago by top1214 • 20

I did have read-groups during alignment. How do I diagnose what the read-group problems are?

ADD REPLY • link written 2.8 years ago by top1214 • 20

How do I diagnose what the read-group problems are?

there is an assertion about RG a line 299: https://github.com/samtools/samtools/blob/master/bam_plcmd.c#L299

ADD REPLY • link written 2.8 years ago by Pierre Lindenbaum ♦ 126k

Treating the samples as different was not actually important to my project, so just using the -R flag to ignore the read group lines solved my particular problem.

ADD REPLY • link written 2.8 years ago by top1214 • 20
Please log in to add an answer.

Similar posts • Search »

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