Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Code Block
languagebash
titleStart an idev session
idev -m 180120 -p normal -A UT-2015-05-18OTH21164 -N 1 -nr 68CoreNGSday4

Go ahead and load the bowtie2 module so we can examine some help pages and options.

Code Block
languagebash
module biocontainers
module load bowtiebowtie2

 Now that it's loaded, check out the options. There are a lot of them! In fact for the full range of options and their meaning, Google "Bowtie2 manual" and bring up that page (http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml). The Table of Contents is several pages long! Ouch!

...

Code Block
languagebash
# Get the FASTQ to align
mkdir -p $SCRATCH/core_ngs/alignment/fastq
cp $CORENGS/alignment/*fastq.gz $SCRATCH/core_ngs/alignment/fastq/

# Set up the bowtie2 index
mkdir -p $SCRATCH/core_ngs/references/bt2/vibCho
cp $CORENGS/idx/bt2/vibCho/*.*  $SCRATCH/core_ngs/references/bt2/vibCho

Make sure you're in an idev session with the bowtie2 BioContainers module loaded:

Code Block
languagebash
idev -p normal -m 120 -A UT-2015-05-18OTH21164  -N 1 -nr 68CoreNGSday4
module load biocontainers
module load bowtie

Now set up a directlry directory to do this alignment, with symbolic links to the bowtie2 index directory and the directory containing the FASTQ to align:

...

We'll be aligning the V. cholerae reads now in ./fq/cholera_rnaseq.fastq.gz (how many sequences does it contain?)

...