Changes between Version 4 and Version 5 of GoNL_Immunochip_Data_Preparation


Ignore:
Timestamp:
Apr 21, 2011 4:39:41 PM (14 years ago)
Author:
laurent
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • GoNL_Immunochip_Data_Preparation

    v4 v5  
    4141
    4242= Liftover file =
    43 The last step in preparing the immunochip data for comparison with the sequence data is to liftover the VCF file to the same Human Genome Reference as the Sequence data so that comparisons can be made. Here is how: 
     43The last step in preparing the immunochip data for comparison with the sequence data is to liftover the VCF file to the same Human Genome Reference as the Sequence data so that comparisons can be made. Here is how:
    4444
    45  1. # Get the appropriate chain files
     45 1. Get the appropriate chain files
    4646   * From the [ftp://gsapubftp-anonymous@ftp.broadinstitute.org Broad GSA ftp] (password is blank)
    47    * From the [http://hgdownload.cse.ucsc.edu/downloads.html#human USCS Genome Browser] # 
    48  1. Get the appropriate fasta files (you'll need both from-and-to builds fasta files) 
    49    * From [ftp://ftp.ncbi.nlm.nih.gov/genomes/H_sapiens/ NCBI] 
    50    * From [http://hgdownload.cse.ucsc.edu/downloads.html#human USCS Genome Browser] 
    51  1. Index the fasta files appropriately to get .fai (samtools) and .dict files (Picard) as described on the [http://www.broadinstitute.org/gsa/wiki/index.php/Preparing_the_essential_GATK_input_files:_the_reference_genome GATK wiki] 
     47   * From the [http://hgdownload.cse.ucsc.edu/downloads.html#human USCS Genome Browser] #
     48 1. Get the appropriate fasta files (you'll need both from-and-to builds fasta files)
     49   * From [ftp://ftp.ncbi.nlm.nih.gov/genomes/H_sapiens/ NCBI]
     50   * From [http://hgdownload.cse.ucsc.edu/downloads.html#human USCS Genome Browser]
     51 1. Index the fasta files appropriately to get .fai (samtools) and .dict files (Picard) as described on the [http://www.broadinstitute.org/gsa/wiki/index.php/Preparing_the_essential_GATK_input_files:_the_reference_genome GATK wiki]
    5252 1. Get and run the [http://www.broadinstitute.org/gsa/wiki/index.php/LiftOverVCF.pl GATK LiftOver tool]
    5353
    54 Note that once the liftover VCF has successfully been created, it can be used to liftover the PLINK files. To do so: 
     54Note that once the liftover VCF has successfully been created, it can be used to liftover the PLINK files. To do so:
    5555
    56  1. Remove all SNPs that are not present in the new reference VCF file (using plink --extract) 
     56 1. Remove all SNPs that are not present in the new reference VCF file (using plink --extract)
    5757 1. Use the liftover VCF as an input to the ''liftover-bim.pl'' tool .