HTSeq through module rcbio/1.3.3

Start an interactive job, with a walltime of 2 hours, 2000MB of memory. 

Make alignment using: STAR through module rcbio/1.1 , bowtie2 through module rcbio/1.1 or bwa through module rcbio/1.1 

After alignment finishes, go to that folder, for example star

Copy the example htseq bash script:

Now you can modify the options as needed. For example, if your file are huge, you may need increase mem and runtime -t for sbatch command. 

To edit the htseq bash script:

To test the pipeline run the following command if you did alignment by bwa as on this page. If you did alignment using bowtie2 following this page, please change -p bwaOut to -p bowtieOut, or if you used star following this page, please change -p bwaOut to -p starOut. Jobs will not be submitted to the scheduler.

To run the pipeline:

To understand how 'runAsPipeline' works, how to check output, how to re-run the pipeline, please visit: Run Bash Script As Slurm Pipeline

Now you are ready to run an rcbio workflow

To instead run the workflow on your own data, transfer the sample sheet to your local machine following this wiki page and modify the sample sheet. Then you can transfer it back to O2 under your account, then go to the build folder structure step.