如何将提取domain的脚本修改为提取所有domain？

基因家族分析课程脚本中

#截取得到序列上的保守结构域序列，注意脚本提取的是第一个domain，如要提取其他domain，请修改脚本27行$a[9]==1为第一个，$a[9]==2为第二个，依次类推

请问如果想把脚本（script/domain_xulie.pl ）修改为提取所有domain应该怎么修改呢？

默认排序时间排序

1 个回答

安生水 2019-07-12 09:35

擅长：perl,基因家族,linux,chip-seq

你可以提取出所有domain，然后合并到一块，注意基因ID需要修改一下。