just wondering if anyone else in Galaxy land was having this problem, I emailed tophat.cufflinks@ to let them know
We have been experiencing segmentation faults with cufflinks 1.0.3. Our platform is CentOS 5.5, and we compiled from source using gcc 4.1.2 (Red Hat 4.1.2-48).
The segmentation fault occured when we passed in a GTF annotation file with the -G option. Using gdb, I tracked it down to line 218 in bundles.cpp:
string rs = (rna_seq) ? rna_seq:"";
It appeared that in our case rna_seq was uninitialized and did not point to a valid C string.
char* rna_seq;
if (loadSeqs && faseq) {
rna_seq = rna.getSpliced(faseq, false, &seqlen);
}
...
if (loadSeqs) {
string rs = (rna_seq) ? rna_seq:"";
...
}
As you can see, it would be possible for the first if-block to not execute, but the second one to still execute (loadSeqs = true, faseq = false). I tried making the following change to the declaration of rna_seq to initialize it to NULL:
char* rna_seq = 0;
this stopped the segmentation fault we were seeing
--
Glen L. Beane
Senior Software Engineer
The Jackson Laboratory
(207) 288-6153