博文
用cut命令从n多列中提取特定的几列
|||
例如：cut -f 1,2,3,4,5,6,7,8 Cul_Wild_1529lines.chr2>1529snp_chr2
这个文件有1529+8=1537列，前8列如下：
Chromosome | Position | RefBase | SnpBase | TotalHitNum | MAF | RefNum | SnpNum |
|---|
后1529列为不同品种在相应位置的SNP。文件Cul_Wild_1529lines.chr2的大小为2.3G，后面产生的1529snp_chr2的文件大小为33.2M。这个就小多了，可以把所有数据导入一个SQL或者MySQL数据库。就这样一个看似简单的命令，完成的工作，可不那么简单。这可能是Linux的魅力。
https://m.sciencenet.cn/blog-479743-745965.html
上一篇：blast in Linux
下一篇：变异效应预测工具Variant Effect Predictor
