GBS722 2013

From UABgrid Documentation
(Difference between revisions)
Jump to: navigation, search
(Ref-based Protocol)
(Ref-based Protocol)
Line 30: Line 30:
  
 
== Ref-based Protocol ==
 
== Ref-based Protocol ==
# '''QC''': Demo Tools > FastQC (4x)
+
# For each sample
# '''Mapping (2x)''':  
+
## '''QC''': Demo Tools > FastQC (2x: forward & reverse reads)
## Demo Tools > BWA for Illumina (ref from history, paired end)
+
## '''Mapping (2x)''':  
## SAM Tools > SAM-to-BAM (ref from history)
+
### Demo Tools > BWA for Illumina (ref from history, paired end)
## rename & download .bam & .bai for IGV visualization
+
### SAM Tools > SAM-to-BAM (ref from history)
## '''QC''': SAM Tools > Flag Stat  
+
### rename & download .bam & .bai for IGV visualization
# '''SNP calling (2x)''': NGS: SAM Tools > Generate VCF with mpileup piped through bcftools view
+
### '''QC''': SAM Tools > Flag Stat  
# '''SNP assessment (2x)''': Demo Tools > SnpEff (Vaccinia_Western_Reserve (vacwr.NC_006998.1))
+
## '''SNP calling (2x)''': NGS: SAM Tools > Generate VCF with mpileup piped through bcftools view
 +
## '''SNP assessment (2x)''': Demo Tools > SnpEff (Vaccinia_Western_Reserve (vacwr.NC_006998.1))
 
# '''Compare Genomes''': Join, Subtract and Group > Compare two Datasets (columns 1-5)
 
# '''Compare Genomes''': Join, Subtract and Group > Compare two Datasets (columns 1-5)
  

Revision as of 11:55, 14 March 2013

Contents

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. For each sample
    1. QC: Demo Tools > FastQC (2x: forward & reverse reads)
    2. Mapping (2x):
      1. Demo Tools > BWA for Illumina (ref from history, paired end)
      2. SAM Tools > SAM-to-BAM (ref from history)
      3. rename & download .bam & .bai for IGV visualization
      4. QC: SAM Tools > Flag Stat
    3. SNP calling (2x): NGS: SAM Tools > Generate VCF with mpileup piped through bcftools view
    4. SNP assessment (2x): Demo Tools > SnpEff (Vaccinia_Western_Reserve (vacwr.NC_006998.1))
  2. 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
Personal tools
Namespaces

Variants
Actions
Navigation
Toolbox