Hi Jim
I am trying to prepare a bed file to load as accustom track on the UCSC genome
browser.
I have a data frame that looks like the one below.
> x
V1 V2 V3
1 chr1 11255 55
2 chr1 11320 29
3 chr1 11400 45
4 chr2 21680 35
5 chr2 21750 84
6 chr2 21820 29
7 chr2 31890 46
8 chr3 32100 29
9 chr3 52380 29
10 chr3 66450 46
I would like to insert the following 4 lines at the beginning:
browser position chr1:1-10000
browser hide all
track type=wiggle_0 name=sample description=chr1_sample visibility=full
variableStep chrom=chr1 span=1
and then insert 2 lines before each chromosome:
track type=wiggle_0 name=sample description=chr2_sample visibility=full
vriableStep chrom=chr2 span=1
The final result should be tab delimited file that looks like this:
browser position chr1:1-10000
browser hide all
track type=wiggle_0 name=sample description=chr1_sample visibility=full
variableStep chrom=chr1 span=1
chr1 11255 55
chr1 11320 29
chr1 11400 45
track type=wiggle_0 name=sample description=chr2_sample visibility=full
variableStep chrom=chr2 span=1
chr2 21680 35
chr2 21750 84
chr2 21820 29
track type=wiggle_0 name=sample description=chr3_sample visibility=full
variableStep chrom=chr3 span=1
chr3 32100 29
chr3 32170 45
chr3 32240 45
Any kind of help or guidance will be much appreciated.
Joseph
____________________________________________________________________________________
Be a better friend, newshound, and
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
