This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
wget http://supportres.illumina.com/documents/downloads/productfiles/humanmethylation450/humanmethylation450_15017482_v1-2.csv |
Next, look at the manifest file and determine what columns are for the chromosome and the base pair coordinates (usually some name like Chr and Coordinate). Put these column numbers into the below line of UNIX code and you will have a .bed file from the manifest file. Here I am making a file only for the X chromosome, but this can be modified to fit your respective needs as well.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
cut -d "," -f15,16 humanmethylation450_15017482_v1-2.csv | sed -e 's/,/\t/' | grep "^X" | sort -n -k2 | awk 'BEGIN{OFS="\t"};{print "chr"$1,$2,$2}' > chrX_HumanMethylation450.bed |
You should be all set with a .bed file from your manifest file. Just ensure that when you are making comparisons you are using the correct genome build.
No comments:
Post a Comment