Question: Pipeline BWA aln+samse
0
gravatar for Mick
7 months ago by
Mick10
Mick10 wrote:

Hi I'm trying to pipeline the bwa aln and samse commands.

This works:

bwa samse reference.fa <(bwa aln reference.fa test.fq) test.fq > test.sam

But this doesn't:

bwa aln reference.fa test.fq | bwa samse reference.fa test.fq > test.sam

Why is that, what is the difference? Thank you.

bwa • 227 views
ADD COMMENTlink modified 7 months ago by ATpoint38k • written 7 months ago by Mick10
1
gravatar for ATpoint
7 months ago by
ATpoint38k
Germany
ATpoint38k wrote:

Wrong order of arguments or rather the sai argument missing:

$ bwa samse
Usage: bwa samse [-n max_occ] [-f out.sam] [-r RG_line] <prefix> <in.sai> <in.fq>

bwa aln reference.fa test.fq | bwa samse reference.fa - test.fq > test.sam

By the way, doesn't work is no error message.

ADD COMMENTlink modified 7 months ago • written 7 months ago by ATpoint38k

Ah ok, I had no idea you could use the dash that way. Now it makes sense. Thank you!

ADD REPLYlink written 7 months ago by Mick10

You can also use /dev/stdin instead in this case but most tools where Heng Li had his hands on, e.g. bwa, samtools, seqtk interpret - as read from stdin.

ADD REPLYlink written 7 months ago by ATpoint38k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1585 users visited in the last hour