Versions Compared

Key

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

...

The Tufts CNR Genomics Core supplies links to bioinformatics resources related to their operation. See Tufts CNR Genomics Core Resources for more information.

c. Genome Indexes on Cluster

Several mammalian genomes, indexes, and annotations are located on the Tufts HPC cluster.  Currently the genomes are listed below in the indicated directory tree are UCSC genome builds, except for canFam3 which is a NCBI build.

/cluster/tufts/genomes
  /HomoSapiens
    /hg18
    /hg19
  /MusMusculus
    /mm9
    /mm10
  /RattusNorvegicus
    /rn4
  /CanisFamiliaris
    /canFam2
    /canFam3

  Within each build subdirectory, there are two subdirectories.

  /Annotation
  /Sequence

In the Annotation directory there are subdirectories for gene annotations ( Gene), and depending upon the degree of annotation, directories for smallRNA and Variation. 

Under the Sequence directory, there are subdirectories containing indexes for popular short read sequence mapping programs.

  /AbundantSequences -- data files with over-represented sequences 
  /BlastDB  -- blast formatted genomic indexes: use genome.fa as reference name
  /Bowtie2Index  -- Bowtie2 formatted indexes: use genome as reference name
  /BowtieIndex  -- Bowtie formatted indexes: use genome as reference name
  /BWAIndex  -- BWA formatted indexes: use genome.fa as reference name
  /Chromosomes  -- individual chromosomes as fasta files
  /Transcriptome  -- Bowtie2 formatted index of transcriptome sites: use transcript as ref name
  /WholeGenomeFasta -- Genome as one file with accessory files

Please read the documentation for a mapping program to understand the way in which the reference indexes are referred.

Example: BWA

It helps to set up environmental variables to avoid having to type long paths. Here a set of short reads ( myreads.fq) are mapped to the mouse genome (mm10) with a SAM formatted file as output. Note that bwa uses genome.fa as a reference index name.

  module load bwa/0.7.9a
  
  export MM10=/cluster/tufts/genomes/MusMusculus/mm10/Sequence/BWAIndex/genome.fa
  export MYDATADIR=/cluster/shared/myutln/mmdata
  bwa mem $MM10 $MYDATADIR/myreads.fq >$MYDATADIR/myreads.sam

Example: Bowtie2

Similarly, environmental variables can be set up, and in the case of bowtie2 BOWTIE2_INDEXES must be set also. Here we have an example of a paired end analysis, with minimal options. Note Bowtie2 uses genome as reference index name.

  module load bowtie2
  export BOWTIE2_INDEXES=/cluster/tufts/genomes/MusMusculus/mm10/Sequence/Bowtie2Index
  export MYDATADIR=/cluster/shared/myutln/mmdata
  
  bowtie2 -q -x genome -1 $MYDATADIR/myreads_1.fq -2 $MYDATADIR/myreads_2.fq -S myreads.sam