Segfault from bedtools in prepare_results.R about leafcutter HOT 7 CLOSED

this is pretty odd. Every time I've run prepare_results.R with the hg19 bed files I get two warnings about inconsistent records but I've never had a seg fault. Have you updated to the latest version of bedtools?

If these weird chromosome names are the case then maybe we could get rid of them when creating the bed files?

from leafcutter.

That warning may just be a coincidence, not sure. I'll try updating and see if it helps!

from leafcutter.

Hmm so upgrading bedtools from v2.25.0 to v2.26.0 didn't fix this.

from leafcutter.

I removed the "weird" chromosome names from gencode_hg19_all_introns.bed.gz. I don't get the warning anymore but I still get the segfault, so I don't think this was actually the root of the problem. @jackhump can you try running on your machine and see if you the same problem?

from leafcutter.

I tried with a new annotation_code generated with wrangle_annotation.sh and still have the same issue.

from leafcutter.

So I'm starting to think this is a bug in bedtools, I've raised an issue here: arq5x/bedtools2#576

from leafcutter.

I've reworked prepare_results.R to use dplyr joins instead of bedtools now so we avoid these issues. I've also changed the CLI to be a bit simpler to document.

from leafcutter.