作业帮 > 综合 > 作业

编写一段关于blast提取序列的Perl语言

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/04/28 13:26:22
编写一段关于blast提取序列的Perl语言
目的:根据信息,提取名称
举例:
.
大概就是这样的,“*”是其他信息,主要就是在“20140101 seqs1”和“20140102 seqs1”这里,blast结果是seqs1的20140101和20140102都提取出来,而blast结果是seqs2的不提取.
$  grep -n “seqs1”  |less如果想输出到某个文件 |less 改为 >youwant.txt
再问: ��д��ϸ�����ҵ�perl��Ϊ0.0000000001. �鷳�ˣ�����
再答: ����˵�ɣ� grep ��liunx�²��ҵ�һ�����grep 'seqs1'����˵�ҵ�����seqs1����,Ȼ��test.txt������㽨�ģ�����������ļ����������ij�����ļ�·���Ϳ��� �����ǵ�Ŀ��������������ֻ���뿴��| less �ǽ�֮ǰ�õ��Ľ���ҵ�����seqs1���У���ʾ�������Ļ������һ�� >youwant.txt��˵ָ��һ����Ϊyouwant.txt���ļ���Ϊ������������ij�����Ҫ�����·������,-n��ʾѡ���������Ƕ��٣�����Ҫ����ȥ����ֱ�������liunx���������� grep 'seqs1' test.txt |less ����˵���뵥��Ҫһ��peri�Ĵ��룿
再问: ��������Ҫһ��perl�Ĵ��룬�ҵ�liunx=0.���ԡ�����������
再答: #!usr/bin/perluse warnings;use strict;use Getopt::Long;my $testFile;GetOptions("testFile:s" => \$testFile,     );print "not exist testFile\n" if ( ! -e $testFile);open (IN,"<$testFile") or die ;while (my $line=<IN>) {     if ($line=~/(^.+seqs1)/){         my $seqs1=$1;         print "$seqs1\n";     } }����Ӧ�ú�����ˣ���windowϵͳ�� ��windowȻ������cmd��֮��ͼperl �������Ŀ�ִ���ļ������������д�Ĵ��룬·�ߺ����ָij����Լ��ģ�-testFile ������ļ���·�����ij����Լ��Ϳ��ԡ�