Error in pipe output to samblaster from bwa-mem2

Error in pipe output to samblaster from bwa-mem2

0

Hi,

I am trying to upgrade my command from bwa to bwa-mem2. This command usually works.

bwa mem -M -R “@RGtID:idtSM:sampletLB:lib” human_g1k_v37.fasta
sample.1.fq sample.2.fq | samblaster -M –excludeDups –addMateTags
–maxSplitCount 2 –minNonOverlap 20 | samtools view -S -b – > sample.bam

Once I start using bwa-mem2, the piping breaks.

bwa-mem2 mem -M -R “@RGtID:idtSM:sampletLB:lib” human_g1k_v37.fasta
sample.1.fq sample.2.fq | samblaster -M –excludeDups –addMateTags
–maxSplitCount 2 –minNonOverlap 20 | samtools view -S -b – > sample.bam

samblaster: Version 0.1.26
samblaster: Inputting from stdin
samblaster: Outputting to stdout
Executing in AVX2 mode!!
SA compression enabled with xfactor: 8
[E::bwa_set_rg] no ID at the read group line
Important parameter settings:
BATCH_SIZE: 512
MAX_SEQ_LEN_REF: 256
MAX_SEQ_LEN_QER: 128
MAX_SEQ_LEN8: 128
SEEDS_PER_READ: 500
SIMD_WIDTH8 X: 32
SIMD_WIDTH16 X: 16
AVG_SEEDS_PER_READ: 64
samblaster: Input file is empty. Exiting.
samblaster: Premature exit (return code 1).

Do you have any idea to run samblaster from bwa-mem2 ?? Thanks !


mapping


markduplicate


bwa-mem2


samblaster

• 32 views

updated 22 minutes ago by

53k

written 2 hours ago by

0

Read more here: Source link