GBS722 2013

From Cheaha
Revision as of 17:38, 14 March 2013 by Curtish@uab.edu (talk | contribs) (Created page with "= VACV-WR WT vs DRUG protocol in Galaxy = This is an abbreviated protocol cheat-sheet to help students through the class exercises. ==Links== UAB Galaxy - https://galaxy.ua...")
(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

VACV-WR WT vs DRUG protocol in Galaxy

This is an abbreviated protocol cheat-sheet to help students through the class exercises.

Links

UAB Galaxy - https://galaxy.uabgrid.uab.edu

Ref Genome @ NCBI

Reads in UAB Galaxy

IGV Setup

  1. JNLP launch: http://www.broadinstitute.org/software/igv/download
  2. Create VACV-WR genome
    1. Genomes > Create .genome
      • Menu: Genomes > Create .genome file
      • Unique Identifier: NC_006998.1
      • Descriptive Name: VACV-WR
      • FASTA File: NC_006998.fna
      • Gene File: NC_006998.gff

Ref-based Protocol

  1. QC: Demo Tools > FastQC (4x)
  2. Mapping (2x):
    1. Demo Tools > BWA for Illumina
    2. SAM Tools > SAM-to-BAM
    3. rename & download .bam & .bai for IGV visualization
  3. QC: SAM Tools > Flag Stat
  4. SNP calling: NGS: SAM Tools > Generate VCF with mpileup piped through bcftools view
  5. SNP assessment: Demo Tools > SnpEff (Vaccinia_Western_Reserve (vacwr.NC_006998.1))
  6. Compare Genomes: Join, Subtract and Group > Compare two Datasets (columns 1-5)

Denovo Protocol (2x)

  1. Denovo contigs :
    1. NGS: Assembly > velveth (shortPairedRead & shortPairedRead2)
    2. NGS: Assembly > velvetg (MinContigLen=200, UsePairedReads=250)
  2. Map Contigs to Ref:
    1. Demo Tools > BWA SW for Contigs
    2. SAM Tools > SAM-to-BAM
    3. rename & download .bam & .bai for IGV visualization
  3. NCBI Blast un-mapped contigs: http://ncbi.nlm.nih.gov > nucleotide blast megablast

Recommended Visualization

  1. IGV, load tracks
    1. WT BWA_BAM
    2. Drug BWA_BAM
    3. WT Velvet_BWASW_BAM
    4. Drug Velvet_BWASW_BAM
Retrieved from "https://docs.uabgrid.uab.edu/w/index.php?title=GBS722_2013&oldid=4449"