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
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
|