featureCounts segmentation fault

[shi]

There seems to be a problem with accessing reads in the BAM file. featureCounts tried to parse the bam file and read in first few reads, but do not know what happened there that caused seg fault. Could you please send us the link to the bam file you downloaded so we can take a close look?

[santhilalsubhash]

FeatureCounts problem

Same problem here, I have tried with Subread 1.4.5-p1 and 1.4.6-p5 Linux-x86_64 versions. I tried downloading the BAM file 7 times and I do not think it has corrupted during downloading because everything goes well while downloading from CGhub genetorrent (no errors). Please some body help me in resolving this issue. I have also crosschecked the BAM files with samtools, it works fine.



featureCounts -Q 10 -F GTF -a MY_ANNOTATION.gtf -t exon -g gene_id -o mypath/out_counts.txt mypath/celline1.bam

========== _____ _ _ ____ _____ ______ _____
===== / ____| | | | _ \| __ \| ____| /\ | __ \
===== | (___ | | | | |_) | |__) | |__ / \ | | | |
==== \___ \| | | | _ <| _ /| __| / /\ \ | | | |
==== ____) | |__| | |_) | | \ \| |____ / ____ \| |__| |
========== |_____/ \____/|____/|_| \_\______/_/ \_\_____/
v1.4.5-p1

//========================== featureCounts setting ===========================\\
|| ||
|| Input files : 1 BAM file ||
Segmentation fault (core dumped)

[shi]

The latest version is 1.4.6-p1. But I don't think that will make a difference.

Have you tried @GenoMax's suggestion to check if your bam file is corrupted? You may also try to subset your bam file to try to find offending reads as suggested by @dpryan. Or alternatively, you may convert your bam file to a sam file to see if you will work.

Wei

[batel]

Hi all,
Thanks for yo9ur answers!
I'm using version v1.4.6, Is it the latest?
I also tried using Rsubread instead of directly, but it also crashes.
The command is:
featureCounts -p -T 5 -g gene_id -a combinedDB5.lincs.f1.pure.gtf -o fc.txt G28059.KMBC-2.1.bam

The full output is:
========== _____ _ _ ____ _____ ______ _____
===== / ____| | | | _ \| __ \| ____| /\ | __ \
===== | (___ | | | | |_) | |__) | |__ / \ | | | |
==== \___ \| | | | _ <| _ /| __| / /\ \ | | | |
==== ____) | |__| | |_) | | \ \| |____ / ____ \| |__| |
========== |_____/ \____/|____/|_| \_\______/_/ \_\_____/
v1.4.6

//========================== featureCounts setting ===========================\\
|| ||
|| Input files : 1 BAM file ||
Segmentation fault (core dumped)

When I used featureCounts via R The command was:
output <- featureCounts("data/BLCA/0aefd782-d636-4308-b940-63054b37e7b0/G28059.KMBC-2.1.bam", annot.ext="combinedDB5.lincs.f1.pure.gtf", isPairedEnd=TRUE, nthreads=3)

========== _____ _ _ ____ _____ ______ _____
===== / ____| | | | _ \| __ \| ____| /\ | __ \
===== | (___ | | | | |_) | |__) | |__ / \ | | | |
==== \___ \| | | | _ <| _ /| __| / /\ \ | | | |
==== ____) | |__| | |_) | | \ \| |____ / ____ \| |__| |
========== |_____/ \____/|____/|_| \_\______/_/ \_\_____/
Rsubread 1.16.1

//========================== featureCounts setting ===========================\\
|| ||
|| Input files : 1 BAM file ||

*** caught segfault ***
address (nil), cause 'unknown'

Traceback:
1: .C("R_readSummary_wrapper", as.integer(n), as.character(cmd), PACKAGE = "Rsubread")
2: featureCounts("/home/batel/CCLEdata/BLCA/0aefd782-d636-4308-b940-63054b37e7b0/G28059.KMBC-2.1.bam", annot.ext = "/home/batel/lncGTF/combinedDB5.lincs.f1.pure.gtf", isPairedEnd = TRUE, nthreads = 3)

Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
Selection: 1
aborting ...
Segmentation fault (core dumped)

Thank you!

[shi]

that's why I don't think that the problem is at the BAM files..
Do you think it's the BAM file?

Hi @batel,

You do not need to use '-b' option when you provide BAM format input. featureCounts automatically detects input format for you.

Please make sure you are using the latest version (1.4.6-p1). '-b' was an option used in old versions.

If the problem persists after upgrading to the latest version, please provide the complete featureCounts output. This will be helpful for diagnosing the problem.

Wei

[dpryan]

If cufflinks worked with the file then presumably it's a bug in featureCounts. It would be good if you could notify the author.

[batel]

I tried cufflinks and it works great..

that's why I don't think that the problem is at the BAM files..
Do you think it's the BAM file?

[dpryan]

If what GenoMax suggested doesn't indicate where the problem is, then try htseq-count. If that runs into problems too then you know with certainty that you have issues with the BAM files. If it doesn't then maybe it's a featureCounts bug (in which case, try subsetting the file until you can find an alignment that causes the problem).

[GenoMax]

You could have somehow corrupted those BAM files.

Have you tried to look at the problem BAM files (with samtools view) to see if there are any obvious problems? Do the file sizes look ok?

[batel]

I have a lot of BAM files (from CCLE database), and I tried counting them using featureCounts.
Some of the files works great, but many(!) of them featureCounts throws segmentation fault right after the 'Input files : 1 BAM file' line.
I'm using the same annotation file all the time so I guess that's not the problem, I also tried to decrease the number of threads- didn't change anything..
Can someone help me? had the same problem?
my command is:
featureCounts -p -T 5 -g gene_id -a /home/batel/GTF/combinedDB5.f1.pure.gtf -o fc.txt -b G30584.UM-UC-3.1.bam
Thank you!