Question: bam-readcount error bam_header_read
0
gravatar for fangqiong.ling
7 weeks 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 weeks ago by Chris Miller18k • written 7 weeks ago by fangqiong.ling0

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

ADD REPLYlink written 6 weeks 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: 493 users visited in the last hour