RDXplorer

[songsy]

I'm using SVmerge. One of the tools is RDXplorer. I run the run.sh. It first picked q0 items and write into a sam file. And then it uses samtools mpileup to call variants. But it doesn't work, showing that :
[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).

The script truncate the bam file itself and it doesn't move on. What can I do?
Thank you very much

[gprakhar]

Hello,

Something wrong with your bam!
BTW Is it a bam?

--
pg

[eulbra]

I add a com in

if fu.isExist(q_mapped_sam):
com = 'wc -l ' + q_mapped_sam + ' > ' + q_mapped_sam + '.log'
v6.execute(com, debug)
com = glob.SAMTOOLS_PATH + '/samtools view -T ' + reference +' -bS '+ q_mapped_sam + ' > ' + q_mapped_bam
v6.execute(com, debug)
com = glob.SAMTOOLS_PATH + '/samtools mpileup -Sf ' + reference + ' ' + q_mapped_bam + ' > ' + q_pileup
v6.execute(com, debug)

then it work, but I still can not run through the whole programe.

Originally Posted by songsy

I'm using SVmerge. One of the tools is RDXplorer. I run the run.sh. It first picked q0 items and write into a sam file. And then it uses samtools mpileup to call variants. But it doesn't work, showing that :
[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).

The script truncate the bam file itself and it doesn't move on. What can I do?
Thank you very much