GalaxyNgs

From Cheaha
Revision as of 17:09, 18 August 2011 by Curtish@uab.edu (talk | contribs) (Quick Notes on NGS sequencing analysis using our Galaxy app)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search


Attention: Research Computing Documentation has Moved
https://docs.rc.uab.edu/


Please use the new documentation url https://docs.rc.uab.edu/ for all Research Computing documentation needs.


As a result of this move, we have deprecated use of this wiki for documentation. We are providing read-only access to the content to facilitate migration of bookmarks and to serve as an historical record. All content updates should be made at the new documentation site. The original wiki will not receive further updates.

Thank you,

The Research Computing Team

Doing NGS in Galaxy

TOC

active users

* David Crossman / Genetics & CCTS / dkcrossm@uab.edu
* Curtis Hendrickson / CCTS / curtish@uab.edu
* John Osborne / CCTS / ozborn@uab.edu

Transfering large datasets to the cluster

Files > 2G can't be uploaded through a browser. See Docs - UploadLargeData

RNAseq

* Useful RNA-seq tutorial pages by Jeremy on Galaxy Pages 
* Tophat
  * formula for "mean inner pair distance":
    * !HudsonAlpha 50bp paired ends: 150, normally. 
* Cufflinks
  * To run cufflinks so that transcripts have usable gene names, use one of the patched GTF files in the Shared Data Library "Patched GTF annotation files for Cufflinks" (created by David C & Curtis H). These have the gene symbol/name patched into the gene_id attribute. 
  * NOTE: if you supply a GTF of existing annotation, it will NOT discover novel splice variants, as galaxy uses the -G flag instead of the -g flag (Cufflinks docs)
  * BUG: Runs but produces all 0 ouput if there are "_" (underscores) in the chromosome names (TreeShrew problem - names were GeneScaffold_####!) - still under investigation - may alternately have been caused by sheer number of chromosomes (scaffolds)


whole genome DNA

Exome

Command Line Tools

sam_chr_coverage

Compute read coverage by chromosome (By John ozborn@uab.edu) First argument is a SAM file to process Second argument is the reference genome fasta used to create that SAM alignment (Available Genome Abbreves: hg19, mm9, mycoplasma_fh, mycoplasma_m129, vaccinia_wr_genome)

{{{ /share/apps/galaxy/src/read_coverage/sam_chr_coverage sam_reads mm9 }}}


q_gunzip

Take a list of files (or directories) and qsub a gunzip job for every .gz file listed or findable in the directory tree. Author: Curtis curtish@uab.edu Note: Symbolic links are NOT followed!! dir [additional file

patch_GTF_with_gene_map.pl

[not yet available] Update the gene_id field in the GTF to contain the gene name, for use with Cufflinks, so it reports in terms of gene names. Gene names can be pulled from an text map file, pulled from the gene_name attribute w/in the GTF itself. Author: Curtis curtish@uab.edu {{{ /share/apps/galaxy/galaxy-tools/bin/patch_gtf_gene_id -src mm9_head10.gtf [-map UCSC_mm9_fake.txt] -prog 3 }}}