Question: Error with samtools - indexing alignments: Writing to standard output failed: broken pipe - python
0
gravatar for coyot001
16 months ago by
coyot0010
coyot0010 wrote:

I am trying to use a pipeline for aligning data from a technique developed by another lab: https://github.com/SavageLab/dipseq/blob/master/README.md It was written a couple years ago and I think that it's not compatible with the most up to date samtools distribution.

I am trying to debug it and fix it so I can use it for my data analysis. The current error I am having is:

[E::bgzf_flush] hwrite error (wrong size) samtools view: writing to standard output failed: Broken pipe [E::bgzf_close] file write error samtools view: error closing standard output: -1

Does anyone know why it might be incopatible with samtools and how to approach solving this problem?

Thanks Willow

alignment software error • 1.3k views
ADD COMMENTlink modified 15 months ago by Biostar ♦♦ 20 • written 16 months ago by coyot0010

Not sure if this helps but I usually get a broken pipe if my computer sleeps during an analysis. If that's the case you can run it with nohup or use a screen.

ADD REPLYlink written 16 months ago by Jake Warner670

Sounds like your bam files might be too big. See here:

ANGSD hwrite error (wrong size)

ADD REPLYlink written 15 months ago by dr_bantz80

dipseq comes with a bunch of tests and concise instructions how to run them. Did you try the tests? That would be my first suggestion.

If you suspect a newer SAMtools is the culprit, install an old one, they are still available for download.

ADD REPLYlink modified 5 months ago • written 15 months ago by h.mon19k

I had the same error with samtools=1.8 so downgraded to 0.1.19 and this fixed the problem.

ADD REPLYlink written 5 months ago by jnkoberstein0
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: 1615 users visited in the last hour