Question: Samtools mpileup -g and -u flags throwing errors for a specific bam file
0
gravatar for top1214
2.8 years ago by
top121420
top121420 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 COMMENTlink modified 2.8 years ago by Pierre Lindenbaum126k • written 2.8 years ago by top121420
0
gravatar for Pierre Lindenbaum
2.8 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum126k 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 COMMENTlink written 2.8 years ago by Pierre Lindenbaum126k

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

ADD REPLYlink written 2.8 years ago by top121420

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

ADD REPLYlink written 2.8 years ago by top121420

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 REPLYlink written 2.8 years ago by Pierre Lindenbaum126k

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 REPLYlink written 2.8 years ago by top121420
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: 1898 users visited in the last hour