segmentation fault error
1
0
Entering edit mode
19 months ago
kalyani ▴ 10

enter image description here

Hi there,

I have a lot of BAM files and I tried counting them using featureCounts. All the files works great, but these few files throwing these error. 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 who had the same problem?

My command is:

featureCounts -a Homo_sapiens.GRCh37.75.gtf -o count.out -T 8 bams/*.bam

Thank you!

featureCounts • 852 views
ADD COMMENT
0
Entering edit mode

Run samtools quickcheck -qvvv your.bam on the problematic bams. Maybe some got corrupted during mapping.

ADD REPLY
2
Entering edit mode
19 months ago
Meisam ▴ 250

Your code is fine and since only few of the .bam files give you this error it clearly shows your pipeline is not the issue. The files that gives you error are corrupt. You can either go back and redo QC/preprocessing/alignment for those sample OR if you want to find out the problematic lines in the current bam files you can use samtools debuging as already described here by Devon Ryan

ADD COMMENT

Login before adding your answer.

Traffic: 1835 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6