Is it possible to output bwa output directly into samtools view to convert into bam format without using intermediate files? I am running low on hard disk space and I am attempting to automate my script. I would like to convert my bwa output to bam, sort it, and index it. I know the sam-bam conversion can be piped into the sort command, but is it possible for the samtools view to take its input from STDIN?
-
-
figured out on linux at least, if I give it the filename /dev/stdin that will satisfy it since the OS does the work, even though samtools does not appear to be expressly programmed to switch to stdin.
It appears to be working too, which is going to save me some headaches.
--edit
apparently I could also have done this with named pipesLast edited by Kotoro; 07-18-2011, 04:06 PM. -
bwa + samtools have been developed with pipes in mind:Originally posted by Kotoro View Post
HTHCode:$ bwa aln [OPTIONS] [DB] [FASTQ] | bwa samse [OPTIONS] [DB] - [FASTQ] | samtools view -Su - | samtools sort - [PREFIX]
dComment
-
-
Almost, this should work (essentially it is a double process substitution).Originally posted by gavin.oliver View Post
Code:$ bwa sampe [OPTIONS] [DB] <(bwa aln [OPTIONS] [DB] [FASTQ1]) <(bwa aln [OPTIONS] [DB] [FASTQ2]) [FASTQ1] [FASTQ2] | samtools view -Su - | samtools sort - [PREFIX]
Comment
-
Thanks!
The wonders of the command line, eh?
Comment
-
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment