Changes between Version 14 and Version 15 of GoNL_Immunochip_Data_Preparation


Ignore:
Timestamp:
Feb 13, 2012 2:30:43 PM (13 years ago)
Author:
laurent
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • GoNL_Immunochip_Data_Preparation

    v14 v15  
    2525
    2626 * From VCF:  grep -v '!^#' in.vcf | awk '{OFS="\t";print $1,$2,$2+1}' > out.bed
     27 * OR (depending on the ref seq)  From VCF:  grep -v '!^#' in.vcf | awk '{OFS="\t";print $1,$2-1,$2}' > out.bed
    2728
    2829Once you have the input file, simply run ''fastaFromBed'' on it giving the Human Reference corresponding to the chip data as the other input. For more information on ''fastaFromBed'', see the [http://code.google.com/p/bedtools/ BEDTools] Manual.