jj_zyl48的博客

原创 Linux shell 提取wig文件中指定染色体的数据

wig数据通常如下图所示： UCSC对wig格式的解释链接：https://genome.ucsc.edu/goldenPath/help/wiggle.html 以19号染色体为例，提取wig文件中19号染色体的数据： cat brain_1_methy.wig | sed -n '/^variableStep chrom=chr19/,/^variableStep chrom=chr2/p'>brain119.wig sed -i '$d' brain119.wig sed -i '