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

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

ADD REPLYlink written 19 days ago by Chris Miller17k
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: 516 users visited in the last hour