----------------------- Need shift 1 command (have more than 9 arguments): ----------------------- ----------------------- Start of variable list: ----------------------- Sample_ID: LZ_M15 Dataset_DIR: /restricted/projectnb/waxmanlab/kkarri/liver_zonation/infusion Sample_Labels_DIR: /restricted/projectnb/waxmanlab/kkarri/liver_zonation/infusion/Scripts/00_Setup_Pipeline GTF_Files_DIR: /unprotected/projects/waxmanlab/routines/GTF_Files MODE: intersection-nonempty STRANDEDNESS_HTSeq: reverse FEATURE_ID: gene_id SCRIPT_DIR: /restricted/projectnb/waxmanlab/kkarri/liver_zonation/infusion/Scripts/07a_Extract_Counts_HTSeq ANNOTATION_FILE: genes.gtf FEATURE_TYPE: exon ----------------------- End of variable list ----------------------- ========================================================== Starting on : Fri May 11 15:11:23 EDT 2018 Running on node : scc-we2 Current directory : /restricted/projectnb/waxmanlab/kkarri/liver_zonation/infusion/Scripts/07a_Extract_Counts_HTSeq Current job ID : 6178465 Current job name : Step_07a_LZ_M15 Task index number : undefined Parameter for multiple cores : 1 ========================================================== Change dir to scratch directory Print scratch directory location: /scratch/6178465.1.w Loading required modules... List files in scratch directory: total 462M drwx------ 2 kkarri waxmanlab 4.0K May 11 15:11 . drwxrwxrwt. 27 root root 68K May 11 15:11 .. -rw-r--r-- 1 kkarri waxmanlab 0 May 11 15:11 LZ_M15_HTSeq.out -rw-r--r-- 1 kkarri waxmanlab 382M May 11 15:11 LZ_M15_primary_unique.bam -rwxr-xr-x 1 kkarri waxmanlab 81M May 11 15:11 genes.gtf Starting to run my commands Starting samtools and htseq-count samtools view LZ_M15'_sorted'.bam | htseq-count --stranded=reverse --type=exon --idattr=gene_id --order=name --mode=intersection-nonempty - genes.gtf > LZ_M15_HTSeq.out Ending samtools and htseq-count List files in scratch total 1.1G drwx------ 2 kkarri waxmanlab 4.0K May 11 15:14 . drwxrwxrwt. 26 root root 68K May 11 15:24 .. -rw-r--r-- 1 kkarri waxmanlab 236K May 11 15:26 LZ_M15_HTSeq.out -rw-r--r-- 1 kkarri waxmanlab 382M May 11 15:11 LZ_M15_primary_unique.bam -rw-r--r-- 1 kkarri waxmanlab 611M May 11 15:14 LZ_M15_sorted.bam -rwxr-xr-x 1 kkarri waxmanlab 81M May 11 15:11 genes.gtf ========================================================== Finished on : Fri May 11 15:26:01 EDT 2018 0 hours, 14 minutes and 38 seconds elapsed. ==========================================================