diff --git a/genect_rnaseq/geneabundance.sh b/genect_rnaseq/geneabundance.sh index 39aad284f18c21a6c93f1e7c08ac479e893d1c74..29efd1587195b4cd50f5c9b25edc5efa02e0be0d 100644 --- a/genect_rnaseq/geneabundance.sh +++ b/genect_rnaseq/geneabundance.sh @@ -41,7 +41,7 @@ source /etc/profile.d/modules.sh module load subread/1.6.1 export PATH=/project/shared/bicf_workflow_ref/seqprg/bin:$PATH -featureCounts -s $stranded -M --fraction -J --ignoreDup -T $SLURM_CPUS_ON_NODE -p -g gene_name -a ${gtf} -o ${pair_id}.cts ${sbam} +featureCounts -s $stranded -M --fraction -J --ignoreDup -T $SLURM_CPUS_ON_NODE -t exon -g gene_name -a ${gtf} -o ${pair_id}.cts ${sbam} mkdir ${pair_id}_stringtie cd ${pair_id}_stringtie stringtie ../${sbam} -p $SLURM_CPUS_ON_NODE -G ${gtf} -B -e -o denovo.gtf -A ../${pair_id}.fpkm.txt