如何用perl处理bed文件?

发布网友发布时间：2022-05-25 01:27

共2个回答

热心网友时间：2023-11-22 06:26

#!/usr/bin/perl -w
use strict;
die "perl $0 <bedfilename>\n" unless(@ARGV == 1);
open IN,$ARGV[0];
while(<IN>){
chomp;
my @tmp=split/[\t ]+/,$_;
print "$tmp[1]\t$tmp[5]\n";
}
close IN;

将上面代码保存为run.pl,假设你的蛋白质序列文件为a.txt 运行的时候，在命令行输入 perl run.pl a.txt 回车即可追问太感谢了，结果现在出来了，怎么把结果放到一个新的txt文件中啊？

热心网友时间：2023-11-22 06:26

awk '{print $2,$6}' infile > outfile

声明：本网页内容为用户发布，旨在传播知识，不代表本网认同其观点，若有侵权等问题请及时与本网联系，我们将在第一时间删除处理。E-MAIL:11247931@qq.com