Thread: RDXplorer
View Single Post
Old 08-14-2012, 01:01 AM   #3
eulbra
Junior Member
 
Location: Beijing

Join Date: Mar 2012
Posts: 8
Default

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.

Quote:
Originally Posted by songsy View Post
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
eulbra is offline   Reply With Quote