Changes between Version 4 and Version 5 of GoNL_Immunochip_Data_Preparation
- Timestamp:
- Apr 21, 2011 4:39:41 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
GoNL_Immunochip_Data_Preparation
v4 v5 41 41 42 42 = 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: 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: 44 44 45 1. # Get the appropriate chain files45 1. Get the appropriate chain files 46 46 * 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] 52 52 1. Get and run the [http://www.broadinstitute.org/gsa/wiki/index.php/LiftOverVCF.pl GATK LiftOver tool] 53 53 54 Note that once the liftover VCF has successfully been created, it can be used to liftover the PLINK files. To do so: 54 Note that once the liftover VCF has successfully been created, it can be used to liftover the PLINK files. To do so: 55 55 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) 57 57 1. Use the liftover VCF as an input to the ''liftover-bim.pl'' tool .