Passing named pipe to samtools mpileup

[kjpr]

Hi all
This is my first post here although I have been reading discussions for a while and it has helped me no end. I am trying to pass a named pipe to the samtool mpileup command but it complains of a broken pipe.

mkfifo input.fifo

/samtools mpileup -EDSgu -d 5000 -f /Homo_sapiens.GRCh37.73.dna.primary_assembly.fa -L 5000 -F 0.05 input.fifo -> /Variants.vcf &

/samtools view -bq 20 -F 1796 /bamfile.bam | tee input.fifo | /samtools mpileup -EDS -d 5000 -f /Homo_sapiens.GRCh37.73.dna.primary_assembly.fa -L 5000 -F 0.05 -> /testmpileup.txt


The Output shows the following error
Running tool: Run command



tee: input.fifo: Broken pipe
[bam_header_read] EOF marker is absent. The input is probably truncated.
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000

I get the data being piped in as expected after the 'tee' command and the testmpilup.txt file is populated but when I pass the named pipe ie 'input.fifo' to the mpileup command, I get a Broken pipe error message and the Variants.vcf file is not populated. I know I can pass the command twice for samtools view and I probably won't lose that much processing time but I successfully implemented this approach on another part of my pipeline but wondered why I can't do the same here.

Any ideas?