Question: bam-readcount error bam_header_read
0
gravatar for fangqiong.ling
6 months ago by
fangqiong.ling0 wrote:

Hey!

I'm running bam-headcount to get allele frequency from mapped reads. I got the following error messages:

Minimum mapping quality is set to 0 [bam_header_read] EOF marker is absent. The input is probably truncated. [bam_header_read] invalid BAM binary header (this is not a BAM file). Segmentation fault: 11

My bam file was created with

samtools view -hbS sample1.sam > sample1.bam
samtools sort sample1.bam -o sample1_sorted

thus output has been specified as BAM and with header.

Any idea what might have gone wrong?

Thanks! Fangqiong

ADD COMMENTlink modified 6 months ago by Chris Miller18k • written 6 months ago by fangqiong.ling0

Have you verified that your bam looks legit with samtools view sample1_sorted.bam | less

ADD REPLYlink written 6 months ago by Chris Miller18k
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: 918 users visited in the last hour