Hello.
I am new to differential gene expression analysis.
I was trying the tutorial of griffith test data provided in BIOSTAR HANDBOOK.
I have done the alignment using HISAT2 and when I am doing quantification using the 

> featureCounts

I am getting the GZIP ERROR:-2 as follows:

> GZIP ERROR:-2

the data is paired ended so I am using the same command given in the manual of subhead package.

>featureCounts -p -a annotation.gtf -t exon -g gene id -o counts.txt mapping results PE.bam

Also, the path and working directory all are fine, still, I don't know why this error is coming up.
Please suggest as I am stuck at this step and not able to proceed.

Thanks in advance!
It seems the error encounter due to `BAM` file input with some bad blocks. Did you check if the validity of the `BAM` file that is generated correctly and not damaged/corrupted? I'd suggest trying a couple of things:

1. Could you try re-running the command with `SAM` file as input in `featureCounts`? 

2. Could you try running the command on a different BAM file (maybe from publically available data?)

Have a look at a discussion forum regarding similar error: https://groups.google.com/forum/#!topic/subread/S4smWRfBNPM

Hope this helps!