GZIP ERROR:-2 using featureCounts (Subread)
2
0
Entering edit mode
gbreard ▴ 10
@gbreard-18447
Last seen 5.5 years ago

Hello, 

I am running STAR analysis of RNAseq samples on a Mac Terminal and I can go through the mapping and .bam file generation but when I launch my command featureCounts, I get in jam with a message "GZIP ERROR:-2" scrolling down and I can't stop the process with usual "Ctrl + C" . It seems to be a problem related to zip initial format but I am a bit lost (I begin). 

What is weird is that I could run it with previous samples using the same command. 

Here below is the command I typed: 

"featureCounts -T 4 -g GeneID -F 'SAF' -s 1 -M -O -a GB_SAF_20180531.txt -o Sample1_featureCounts.out Sample1_Aligned.sortedByCoord.out.bam & Sample1_counts.log "

Does anybody have an idea how to solve the problem?

Thank you in advance

Gwenael

 

software error rnaseq • 1.9k views
ADD COMMENT
0
Entering edit mode

Subread version number?

You might also interested in: https://www.biorxiv.org/content/early/2018/07/26/377762

ADD REPLY
2
Entering edit mode
Yang Liao ▴ 380
@yang-liao-6075
Last seen 19 days ago
Australia

Hi Gwenael,

The "GZIP ERROR:-2" message sounds like an old bug in featureCounts. Which version of featureCounts did you use? If it was not the latest version (v1.6.3), you can try our latest version, which is available on 

 https://sourceforge.net/projects/subread/files/subread-1.6.3/

since this bug had been fixed in v1.6.3 of featureCounts.

Cheers,

Yang

ADD COMMENT
1
Entering edit mode
gbreard ▴ 10
@gbreard-18447
Last seen 5.5 years ago

Yes!!!

It worked. Thanks a lot.

Gwenael

ADD COMMENT

Login before adding your answer.

Traffic: 463 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